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Abstract 



An Isotropic Diffuse Gamma-Ray Background (IDGRB) in the spectral range 30- 
10,000 MeV was first reported in the early 1970's using measurements made by the 
SAS 2 instrument. Data recorded by the Energetic Gamma Ray Experiment Tele- 
scope (EGRET) on board the Compton Gamma Ray Observatory (CGRO) over the 
last 4 years are analysed in order to extract the best measurement yet made of the 
IDGRB. Extensive analysis of the EGRET instrumental background is presented in 
order to demonstrate that an uncontaminated data set can be extracted from the 
EGRET data. 

A model of the high latitude galactic diffuse foreground emission is presented 
and the existence of an IDGRB is confirmed. Spatial and spectral analysis of this 
background is presented. In addition, point source analysis at high galactic latitudes 
is performed to reveal the existence of a population of extragalactic sources. The 

characteristics of this population are examined and models of its flux distribution 
are reported. The question of whether the IDGRB is composed of unresolved point 
sources is addressed using fluctuation analysis. 

Finally, possible future directions for 7-ray astronomy arc examined through sim- 
ulations of a future 7-ray telescope: the Gamma-ray Large Area Space Telescope 
(GLAST). The GLAST baseline design is described and its scientiflc performance 
is evaluated. The ability of this telescope to detect 1,000-10,000 new extragalactic 
sources is demonstrated and the likely impact on the study of the IDGRB is consid- 
ered. 
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Chapter 1 
Introduction 

1.1 The Isotropic Diffuse Gamma- Ray Background 
(IDGRB) 

If one were to be able to look up at the high energy (E > lOOMeV) 7-ray sky, it would 
not be dominated by thousands of point sources as is the familiar night sky, but rather 
it would be seen to be dominated by a bright diffuse swath cut across the sky by the 
plane of our own Milky Way galaxy. Diffuse 7-rays from our galaxy constitute the 
single brightest source of celestial gamma radiation. But even the regions well away 
from this bright galactic feature would be seen to glow faintly. Since the 1970's and 
the results oi SAS 2 , there has been evidence that the faint 7-ray glow at high galactic 
latitudes is brighter than can be accounted for by the small amount of galactic gas 
observed looking through the plane of our galaxy. This glow appears to be isotropic 
leading to the conclusion that it is likely to be of extra-galactic origin. It is the nature 
of this Isotropic Diffuse Gamma-Ray Background (IDGRB) that is the subject of this 
discussion. 

Of course even if our eyes were sensitive to gamma radiation the spectacular fea- 
tures of the 7-ray sky would remain hidden from us because of the opacity of the 
earth's atmosphere to 7-rays . The emergence of the field of 7-ray astronomy in the 
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last half century has been driven by the development of satellite borne instrumen- 
tation. The most sensitive high energy 7-ray telescope yet developed, the Energetic 
Gamma-ray Experiment Telescope {EGRET ), was launched in 1991 aboard the 
Compton Gamma Ray Observatory (CGRO). This instrument has yielded the best 
quality data to date on high energy 7-ray phenomena. Traditionally, 7-ray astronomy 
has been limited by poor counting statistics. Even objects that release comparatively 
large energies produce relatively few 7-ray photons because of the large energy car- 
ried by each such event. EGRET's large collecting area has allowed unprecedented 
numbers of photons to be collected (over a million photons of celestial origin since 
launch) which has enabled the investigation of the high energy phenomena of the 
universe to be conducted with unprecedented power. 

EGRET is particularly well suited to the study of diffuse radiation. Its relatively 
low altitude and low inclination orbit minimize the instrument's exposure to the 
background cosmic rays and multiple systems on the instrument reject the cosmic 
ray induced events with high selectivity. As will be discussed in Chapter 3, this 
results in a low instrumental background which is critical to measurements of faint 
diffuse features. Furthermore, the existence of a good instrument calibration allows 
one to do accurate absolute photometry. 



1.2 Physical Models 

The details of the measurement of the IDGRB will be discussed in subsequent chap- 
ters. This section is intended to provide the physical motivation for the study of the 
IDGRB as well as a brief history of the physical models invoked in an attempt to 
understand this phenomenon. 



Figure |0| shows a reproduction of the diffuse spectrum reported by Fichtel et al. 
( Pichtel et al. 19781) using SAS 2 data. The background was measured to be a power 
law with a soft spectral index of 2. 7 ±0.4. The integrated intensity above 35 MeV was 
reported to be 5.7 ± 1.3 x lO^^ph cm~^s^^sr^^. Also shown in this figure are balloon 
measurements of the soft 7-ray background. These early measurements indicated 
an enhancement above a continuous power law in the energy range from 1-10 MeV 
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Figure 1.1: Figure from (Fichtel et al. 1978) showing the SAS 2 measurement of the 
diffuse 7-ray background as well as other lower energy balloon measurements. 



( [Trombka et al. 1977|) . These two measurements motivated much of the early work 
on the IDGRB. 
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1.2.1 Baryonic Halo Models 

Interactions between cosmic rays and atomic and molecular gas in the galactic plane 
are responsible for the copious 7-ray production in the plane of our galaxy. This 
process will be discussed in detail in Chapter 4. In addition to the observed matter 
in our galaxy, there is widely believed to be a dark halo which is invoked in order to 
explain the measured galactic rotation curves. If this halo were baryonic in nature it is 
possible that the same sorts of interactions that are responsible for the diffuse galactic 
7-ray emission are occurring at large galactocentric radii leading to a comparatively 
isotropic diffuse source of 7-radiation. 

The proposed halo must meet several requirements for this theory to be applicable. 
Firstly, the halo must be composed of baryonic matter in order to generate significant 
targets for the cosmic rays. Recent measurements by the EROS ( |Aubourg et al. 1993|) 
and MACHO ( |Alcock et al. 19931) collaborations have provided evidence in support 
of this idea. Monitoring of millions of stars in the LMC has resulted in the detection 
of a handful of microlensing events indicating halo objects with masses ~ O.IMq. 
Because it is unlikely that there is a mechanism that transforms gas into massive 
halo objects with perfect efficiency, it is to be expected that some gaseous clouds are 
also present. For these clouds to be invisible to terrestrial radio telescopes, they must 
be composed of molecular hydrogen, H2. 

Lastly, this halo must be exposed to a sufficient flux of cosmic rays for this mat- 
ter to be illuminated. The cosmic ray flux at large galactocentric distances can be 
estimated from galactic mass loss (see e.g. Breitschwerdt 1991 ) but is quite uncertain. 

De Paolis et al. ( PePaolis et al. 1995| ) have suggested a scenario in which all the 
above requirements are met. In this scenario, the same sorts of proto-galactic clouds 
which result in globular clusters can result in molecular clouds at larger galactocentric 
radii. The factor that determines the fate of such a proto-galactic cloud is the intensity 
of UV radiation it is exposed to in the early universe. This emission (due to a 
population of massive young stars in the galactic center) imprints a characteristic 
temperature on clouds at distances less than ~ 10 kpc resulting in their evolution 
into globular cluster. Clouds at larger distances cool and collapse simultaneously and 
consequently are able to fragment into smaller clouds. These clouds would constitute 
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the raw material from which to form massive halo objects as well as targets for cosmic 
rays interactions. 

The same authors have predicted the resultant 7-ray flux that would emerge from 
such a scenario. Assuming that the dark matter halo is entirely baryonic in nature 
and that the cosmic ray spectrum and composition is similar to that measured locally, 
they show, 

7^ ~ e X 1.7 X 10-^ph cm'^ s"^ sr"^ ,E>1 GeV (1.1) 

where e reflects the anisotropy which causes most of the resultant 7-ray fiux to be 
directed outward from the galaxy because of the net radial motion of the cosmic rays. 
Because this factor is not well known and could easily be smaller than 0.1, these 
models can be made consistent with the SAS 2 measurement. 

The best way to unambiguously test this source of the IDGRB would be to detect 
an anisotropy resulting from the sun's offset within this halo. The effect of this offset 
is to make the emission more intense in directions toward the galactic center. The 
SAS 2 measurement revealed no such anisotropy and resulted in lower limits on the 
halo size of > 45 kpc. This is somewhat larger than predicted by the above authors 
but there are large uncertainties associated with measuring the IDGRB flux toward 
the galactic center because the inverse Compton emission is so pronounced in this 
direction. 

1.2.2 Baryon Symmetric Big Bang Theories 

Baryon Symmetric Big Bang (BSBB) Cosmologies are cosmological models in which 
equal amounts of matter and anti-matter arc formed in the big bang. Subsequently, 
there was a separation of the two forms of matter into zones which then evolved into 
matter or anti-matter zones that are the size of galactic clusters. Such cosmologies 
were first introduced in order to explain how a matter symmetric big bang (which is 
aesthetically preferable) could result in the dominance of matter of anti-matter in the 
local universe. While the explanation of this phenomenon now centers on weak CP 
violation in the early universe, experimental constraints on CP violating interactions 
are not sufficient to definitively identify them as the cause of the matter asymmetry 
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in the universe. 

BSBB models are summarized in |Stecker 1978| . In these models, bubbles of matter 
and antimatter form in the early universe due to the phase separation between matter 
and antimatter ( |Cisneros 1973| ). The free energy of the system scales as the ratio of 
the surface area of such a bubble to the volume and as a result the domains tend 
to coalesce to form large matter or anti-matter domains ( pmnes 1972| ). Models to 
explain galaxy formation in such domains have been proposed ( [Stecker 19721 ). 

The main observational consequence of BSBB models is that the interfaces be- 
tween matter/antimatter domains would be the site of copious nucleon-antinucleon 
annihilations. This would generate a diffuse 7-ray flux through the subsequent decay 
of 7r°'s. Because these processes are well understood, it is possible to calculate the 
expected spectrum of such annihilation subject only to cosmological free parameters. 
Figure |1.2| shows a reproduction of a figure from Stecker ( ^tecker 1978| ) in which 
the nucleon-antinucleon annihilation spectrum is fit to the observed IDGRB. The 
main features of the IDGRB spectrum are well reproduced by this model. Both the 
'MeV-bump' and the steep high energy spectrum match the theoretical spectrum. 

While this early agreement between BSBB models and the IDGRB lends credi- 
bility to these theories there are other observational problems. The most important 
additional constraint on BSBB theories is the distortion of the microwave background 
caused by interactions between the supra-thermal electrons that would result from the 
annihilations and the background photons. Electron bremsstrahlung would distort 
the low energy (Rayleigh- Jeans) tail of the microwave background while Thompson 
scattering would distort the high energy (Wien) tail ( [Ramani 1976|) . Measurements 
by the Cosmic Background Explorer (COBE) have shown no evidence of such distor- 
tion. 

A further challenge to this theory results from big bang nucleosynthesis. Because 
the matter/antimatter domains are smaller than the neutron diffusion length at the 
time of nucleosynthesis, the neutrons are annihilated and cannot participate in nu- 
cleosynthesis ( [Combes et al. 1975| ). The He abundance in the universe must instead 
be generated in galaxies in so-called 'little bangs' ( |Wagoner 19741 ). The success of big 
bang nucleosynthesis at predicting He abundances makes this an unlikely explanation. 
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Figure 1.2: Figure from Stecker (1978) showing the theoretical matter-antimatter 
annihilation spectrum fit to the IDGRB. This model reproduces the 'MeV bump' as 
well as the steep spectrum at high energy 

Recent measurements by the COMPTEL instrument aboard CGRO ( |Kappadath| 



et al. 1995|) removed one of the strongest pieces of evidence in favor of this theory. 
The modulation of the diffuse spectrum with orbital rigidity (see chapter 3) has 
allowed the unambiguous separation of instrumental and celestial backgrounds for 
the first time in this energy range. The most recent reports of the soft 7-ray diffuse 
spectrum suggest that the 'MeV-bump' is an instrumental artifact and that the 7-ray 
background in this energy range is instead a power law consistent with the spectrum 
at higher energies. 
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1.2.3 Exotic Particle Decays 

The discovery of the IDGRB led in part to interest in the field of particle astrophysics. 
Very high energy phenomena which are out of the reach of terrestrial accelerators 
could be important processes in the early universe. Observational constraints of the 
early universe provide the best insight available into physics beyond the standard 
model. Many particle theories predict the existence of new massive particles. The 
observable end products of the eventual decay of such particles are 7-rays. These 
7-rays are produced either directly as products of the particle decay or as the end 
products of a baryonic shower with neutral pion production and decay. Attempts 
have been made to explain the IDGRB using light super- symmetric particle decay 
( [Ellis et al. 198^ ), and cosmic string interactions ( [MacGibbon et al. 1993|) . 

These theories have the common theme that that some of the IDGRB is due to 
the decay of Weakly Interacting Massive Particles (WIMPs). If these WIMPs are 
long lived compared to the age of the universe it is possible that they are significant 
contributors to the dark matter of the universe. The 7-ray spectrum expected from a 
WIMP particle X has been calculated by Kamionkowski ([Kamionkowski 1995|) . The 
calculation is outlined here. 

If we assume that in the early universe, the particle X was in thermal equilibrium 
and that T » mx, there would consequently be as many X particles as photons. In 
this state the X particle is constantly undergoing equilibrated conversions to quark 
anti-quark pairs or lepton anti lepton pairs: XX ^ qq. As the universe cools below 
T = mx the density of X is exponentially suppressed by the Boltzmann factor. 
Consequently, if the universe were to remain in thermal equilibrium as it cools there 
would be no remaining massive WIMPs. However, the non-thermal expansion of the 
universe causes the X to 'freeze out'. This happens when the annihilation rate F = 
{o'Av)nx drops below the expansion rate H. Kamionkowski shows that the parameter 
{(Jav), the thermally averaged cross-section for annihilation times the relative velocity, 
determines the relic density that survives the early universe. 

^ ,2 ,3xlO-27cmV\ , , 

= ' — (^) — > ■ 
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where fix is the X density in units of the critical density and h is the Hubble constant 
in units of 100 km s~^Mpc~^. Note that this result is independent of mx and that 
the density increases with decreasing cross section. 

If this relic distribution of WIMPs decays at a characteristic red-shift zd, the 
resultant 7-ray flux can be calculated to be, 

where, 

Eo = = ^(-)^/^ , (1.4) 

21 + ZD 2 4o 

and r is the WIMP lifetime and to is the age of the universe. 

If instead the WIMP lifetime r is O(to) and the decay happens today, 

E-^ ^ 2.4 X l04-_-iki^(^)3/2,^-2 -1 -1 ^^ 5) 

dQdE {T/to){mx/GeVyEo' ^ ' 

for E < mx/2. 

In either case the WIMP annihilation can produce observable levels of background 
with values of Qx as low as 10~^. Figure shows a figure reproduced from Kamion- 



lowski (IKamionkowski 1995|) showing these calculated spectra superimposed on the 
IDGRB spectrum. While it is evident that a single WIMP line cannot account for 
the entire IDGRB, it is possible that a spectral 'bump' in the IDGRB signals the 
presence of a WIMP decay. 

The above discussion describes the scenario in which the WIMP particles are 
unstable on the time scale of the age of the universe. It is perhaps more intriguing 
to consider the case in which the WIMP particles have lifetimes greater than the 
age of the universe. In this case WIMP particles could still be abundant today and 
thus could constitute an important fraction of the dark matter. These stable WIMPs 
would cluster along with luminous matter into galaxies. They are particularly good 
candidates for massive halos. Kamionkowski has noted that the expected order of 
magnitude annihilation cross-section for a particle with weak-scale interactions is, 
{aAv) ~ Q;^/(100GeV)^ ~ lO^^^cm^s"^ Substitution into equation 1.2 yields the 
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Figure 1.3: Figure from Kamionkowski (1995) showing the diffuse 7-ray background 
spectrum enveloped by a conservative upper hmit (dotted hne). Theoretical calcu- 
lation of the diffuse 7-ray spectrum due to primordial WIMP decays in the early 
universe (solid lines) and in today's galactic halo (dashed lines) are shown superim- 
posed. 

suggestive result that the resulting WIMP density is remarkably close to the critical 
density. 

While stable WIMPs would not decay into 7-rays as described above, there is 
a finite cross section that couples the WIMPs and anti- WIMPs to ordinary matter. 
Thus WIMPs in our galactic halo can annihilate just as they did in the early uni- 
verse to produce 7-rays. There are two distinct channels that are of interest. If the 
products of the WIMP annihilation are quarks and leptons, the subsequent shower 
from hadronization of the quarks produce a 7-ray spectrum with a broad peak at 
~ 1/10 of the WIMP rest mass ( [Silk 198^ ). This background could be a significant 



contribution to the IDGRB even though if the halo is large enough to obscure any 



CHAPTER 1. INTRODUCTION 



11 



quadrupole moment caused by the earth's position within the halo. 

While by definition there is not direct coupling between WIMPs and 7-rays, it is 
quite likely that a pair of WIMPs may couple to a pair of photons through Feynmann 



loop diagrams ([Srednicki 1986| ). If this is the case, the low mean velocity of a halo 
WIMP (~ 300km s~^) leads to the production of essentially mono-energetic 7-rays. 
The observation of a 7-ray line at several tens of GeV would be a powerful signature 
of the presence of WIMPs in the galactic halo. 

1.2.4 Active Galaxies 

Prior to the launch of EGRET, COS B had detected one active galaxy in the Virgo 
region: 3C 273 ( [Bignami et al. 1981 ). While this detection did not allow much more 



than speculation as to the mechanisms of particle acceleration in AGN, it was im- 
mediately suggested that active galaxies could contribute significantly to the IDGRB 
( pchonfelder et al. 1978| , [Fichtel et al. 198l| ). Predictions as to the angular distribu- 



tion that would result from such a collection of point sources were also put forward 
( Pao et al. 1990| ). As will be discussed at some length in subsequent chapters, the 



EGRET detection of a large population of these sources have brought these theories 
to the forefront. A large part of this discussion will be devoted to the extraction of 
the amount of background that can be attributed to these sources. 



1.3 Summary 

The relative transparency of the universe to 7-rays as well as their ubiquity in high 
energy processes in the universe give 7-ray astronomers a unique insight into otherwise 
unmeasurable phenomena. The study of the IDGRB which is most likely cosmological 
in nature, can inform astrophysical models as well as fundamental particle physics 
theories. In the following chapters, the EGRET data will be presented and analyzed 
to extract the extragalactic signal. A detailed discussion of the contribution from 
AGN will be presented and prospects for the future study of this phenomenon will 
be discussed. 
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Throughout this discussion, unless otherwise indicated, "flux" will refer to the inte- 
grated photon flux from a point source above 100 MeV in units of photons cm~^ s~^. 
"Intensity" will refer to to the integrated diffuse flux above 100 MeV in units of 
photons cm~^ s~^ sr~^. 



Chapter 2 

The EGRET Instrument 



2.1 Introduction 

High-energy 7-ray observations provide one of the most direct views of energetic pro- 
cesses occurring in the universe. Unhke charged particles, 7-rays will pass undeflected 
through the magnetic fields of interstellar space, retaining their directional informa- 
tion. In addition, the small interaction cross-section of high-energy 7-rays allows 
them to travel essentially unattenuated through space. A high-energy 7-ray can pass 
through the central plane of the galactic disk with only a one percent chance of be- 
ing absorbed. In contrast, an optical photon can only penetrate about one-tenth of 
the distance from the galactic center to the Earth as it travels through the galactic 
disk. Unfortunately, 7-rays cannot penetrate very deeply into the Earth's atmosphere. 
Even using high-altitude balloon observations, the interaction of cosmic rays with the 
Earth's atmosphere produces such a large number of background 7-rays that only 
the strongest point sources can be detected without leaving the Earth's atmosphere. 
Thus, the progress of observational high-energy 7-ray astrophysics was slow until the 
1970's when detectors were first placed on orbiting sateUites. 

Since the fiux of celestial 7-rays is relatively low, detectors must have a large 
sensitive area, a high detection efficiency, and they must be able to distinguish a 
7-ray from charged particle cosmic radiation, whose flux can be a factor 10^ more 
than that of 7-rays. In addition, the detector must measure the incident direction. 
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arrival time, and energy of the incoming photon. Using a spark chamber assembly 
surrounded by an anti-coincidence dome, SAS 2 was able to efficiently reject charged 
particle radiation and achieve an angular resolution of ~ 2° ( [Fichtel et al. 1975| ). 
COS B was similar in area to SAS 2, but it also incorporated a crystal scintillator to 
aid in the measurement of the 7-ray energy ( [Bignami et al. 1975|) . The results from 
these two experiments moved high-energy 7-ray astronomy past the discovery phase 
and into an exploration phase. 

In the mid 1970's, NASA began the Creat Observatories for Space Astrophysics 
program for the express purpose of mapping the electromagnetic spectrum from in- 
frared to 7-rays with unparalleled detail. As part of this program, the Compton 
Gamma Ray Observatory (CCRO) was approved in 1977 to explore the 7-ray win- 
dow from less than 0.1 MeV to more than 10 GeV. Since substantially different 
detection methods are required to detect photons in different parts of the 7-ray spec- 
trum, four instruments were selected to be placed on board CCRO. The Burst and 
Transient Source Experiment (BATSE) is designed to detect 7-ray bursts over the en- 
ergy range 25 keV-2 MeV and consists of eight uncollimated detector modules placed 
on the eight corners of CCRO, providing nearly uniform coverage of the sky ([Fish 



man et al. 1989|) . The Oriented Scintillation Spectrometer Experiment (OSSE) uses 
four shielded NaI(Tl)-CsI(Na) phoswich detectors to study astrophysical sources in 
the 0.05-10 MeV energy range, with an effective field-of-view of 3?8 x 11?4 full width 
at half- maximum (FWHM) ([Johnson et al. 1993|) . The imaging Compton Telescope 
(COMPTEL) explores the 1-30 MeV energy range with a field of view of ~ 1 sr, 
relying on a Compton scattering in one detector array, and a second interaction in 
a lower detector array ( [Schonfelder et al. 1993 ). The ECRET instrument is a pair- 
conversion telescope, sensitive to photons in the energy range 20 MeV to 30 GeV, 
with a field of view of ~ 1 sr. 

CCRO was launched aboard the Space Shuttle Atlantis (STS-37) from the Ken- 
nedy Space Center in Florida on 1991 April 5 and was deployed two days later. 
ECRET was activated on 1991 April 15 and began taking data on April 20. The 
ECRET instrument and its calibration have been described extensively in Hughes 
et al. (1980), Kanbach et al. (1988, 1989), Nolan et al. (1992), and Thompson et al. 
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(1993). The EGRET instrument has been designed, cahbrated, and maintained in a 
collaborative effort by scientists from the following institutions: 

• NASA/Goddard Space Fhght Center, Greenbelt MD, U.S.A. 

• Stanford University, Stanford CA, U.S.A. 

• Max-Planck Institut fiir Extraterrestrische Physik, Garching, FRG. 

• Grumman Aerospace Corporation, Bethpage NY, U.S.A. 

2.2 Design of the Detector 

EGRET is multilevel spark chamber triggered by a scintillator coincidence system, 
with a large Nal(Tl) crystal spectrometer to measure the energy. The instrument 



is shown schematically in Figure |2.1| . It is similar in design to SAS 2 and GOS B, 
but it employs a much larger spark chamber with improved resolution. Above the 
critical energy of a given material, generally on the order of a few tens of MeV, pair 
production becomes the dominant process by which 7-rays interact with matter. Pair 
production is an attractive interaction to observers because the resulting electron- 
positron pair retains almost all of the information about the incident direction and 
energy of the original photon. The EGRET telescope is designed so that the electron- 
positron pairs produced by high-energy 7-rays can be clearly separated from other 
events that might trigger the detector. 

A 7-ray which enters the top of the EGRET instrument will pass undetected 
through the large anticoincidence scintillator surrounding the spark chamber and has 
a probability ~ 33% of converting into an electron-positron pair in one of the thin 
tantalum (Ta) sheets interleaved between the 28 closely spaced spark chambers in the 
upper portion of the instrument. The effective critical energy of Ta is 10.7 MeV. The 
high atomic number Z of Ta increases the probability of pair production relative to 
the Compton process and reduces the energy loss per unit length by the converted e^, 
while the thickness of the sheets is kept to a minimum to reduce multiple scattering. If 
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Figure 2.1: Schematic of the EGRET instrument. 



pair production occurs, the secondary charged particles will ionize the spark chamber 
gas along their trajectories. 

Below the conversion stack are two 4x4 arrays of plastic scintillation detector tiles 
spaced 60 cm apart which register the passage of charged particles. The general di- 
rection of a charged particle passing through both scintillator arrays is determined by 
measuring the time-of-flight delay with a resolution of ~ 1.5 ns, and by the combina- 
tion of individual upper and lower scintillator tiles that detected the charged particle. 
In between the upper and lower arrays of scintillator tiles are six widely spaced spark 
chamber modules with very little scattering material. These provide additional infor- 
mation on the path of a charged particle as it travels between the upper and lower 
scintillator arrays. If the time-of-flight delay indicates a downward moving particle 
which passed through a valid combination of upper and lower scintillator tiles, and 
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the anticoincidence system has not been triggered by a charged particle, a high volt- 
age pulse is applied to orthogonal wire arrays in the spark chamber modules, and the 
track information is recorded digitally in ferrite cores, which are then read out and 
included in the instrument telemetry. In this manner, a three-dimensional picture of 
the path of the electron-positron pair is measured, allowing it to be identified with 
certainty and its basic properties to be determined. The high- voltage trigger and path 
readout will occur even if only one charged particle is detected by the time-of-flight 
coincidence system. 

The energy deposition in the Nal(Tl) Total Absorption Shower Counter (TASC) 
located directly below the lower array of plastic scintillators is used to estimate the 
photon energy. The TASC has a thickness of about eight radiation lengths and 
can determine photon energy from a few tens of MeV up to several tens of GeV 
with an energy resolution of ~ 20% FWHM over most of that range. TASC energy 
measurements must be corrected for the energy lost by the while traveling through 
the spark chambers and scintillators. This correction is estimated from the path 
length in each material and is typically in the range 15-30 MeV. Another energy 
correction is made to those photons for which one of the electron-positron pair was 
not intercepted by the TASC. Once an event has been registered, it is tagged with 
an arrival time by the on-board clock which measures Coordinated Universal Time 
(UTC) with an absolute accuracy 100 fis, and has a relative accuracy of 8 /is. This 
clock also serves the other three instruments aboard CGRO. The spark-chamber 
tracks, energy measurement, and 7-ray arrival time thus comprise the information 
recorded for a single event. There is readout dead time of ~ 100 ms after each 
triggered event. 

The telescope assembly is surrounded by a thermal blanket, which is necessary to 
maintain sufficient temperature control for the anti-coincidence dome. The thermal 
blanket also serves to shield the anti-coincidence dome from micro-meteors which 
might hit the dome and create a light leak, disabling the charged particle anti- 
coincidence system. Micro-meteors are abundant in a low Earth orbit, consisting 
largely of paint flakes from spacecraft. A thermal blanket thickness of 0.17 g cm~^ re- 
duces the chance of a micro-meteor penetration to less than 1% over the course of the 
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Table 2.1: Comparison of the On- Axis Performance of the High-energy 7-Ray Detec- 
tors 

SAS 2 COS B EGRET^ 
Field of View 0.25 sr 0.25 sr 1.0 sr 

Effective Area 100 cm^ 70 cm^ 1200 cm^ 

> 100 MeV 

Angular Resolution 1?5 1?5 0?6 

RMS @ 500 MeV 

Energy Resolution ~ 100% 42% 18% 

FWHM @ 100 MeV 

Point Source Sensitivity 10"^ 10"^ 10"^ 
(photons cm~^ s~^) 

> 100 MeV 
10^ s exposure 



mission and at the same time does not introduce a significant amount of background 
7-rays generated by cosmic-ray interactions in the thermal blanket. 

With its much larger effective area and better angular resolution, EGRET has 
more than an order of magnitude greater sensitivity than SAS 2 and COS B. Some 



of the basic properties of the detectors are summarized in Table 2A. The improved 
angular resolution of EGRET comes from the use of thinner conversion plates as 
well as better track recognition. In addition, SAS 2 and COS B used less sensitive 
Cerenkov detectors to determine the direction of the secondary charged particles. 



2.3 Instrument Calibration 

Much of the EGRET calibration was performed in 1986 at the Stanford Linear Ac- 
celerator Center (SLAC). A beam of electrons with an adjustable energy between 
650 MeV and 30 GeV was back-scattered off pulsed laser photons to produce a beam 
of inverse-Compton scattered 7-rays from 15 MeV to 10 GeV ( [Mattox et al. 1987| ). 
This back-scattered beam was coUimated to a cross-section of 1 cm^ at a distance of 
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170 cm from the interaction region in order to minimize the angular dispersion. The 
beam was then aimed at the EGRET instrument to determine the effective sensitive 
area, as well as the angular and energy dispersion. 

Calibration runs were performed for ten discrete energies (15, 20, 35, 60, 100, 200, 
500, 1000, 3000, and 10,000 MeV) with the beam incident at each of five inclination 
angles (0°, 10°, 20°, 30°, and 40°) and three azimuth angles (0°, 22?5, and 45°). 
Over the course of the calibration, more than 500,000 events were recorded. It was 
found that there is very little dependence on azimuthal angle. A smooth instrument 
response matrix was formed by applying a two-dimensional fourth-order Chebyshev 
fit with respect to energy and second-order with respect to the inclination angle to 
the results of the calibration runs. This fit was stored in three sets of calibration files, 
described below, which form the basis for all analysis of EGRET data. Thompson 
et al. (1993) have performed in-fiight verification tests to ensure that the quantities 
determined during the calibration runs are still valid. 

The sensitive area files contain information on how the instrument effective area 
varies with incident energy and inclination angle. For an on-axis source, the effec- 
tive area is greater than 1000 cm^ from 100 MeV out to 3 GeV. Not surprisingly, 
the effective sensitive area decreases as the aspect angle increases, such that at an 
inclination of 30°, the sensitive area is less than 15% of the on-axis effective area. 

The angular dispersion files describe the instrument point spread function (PSF) 
at a particular aspect angle and energy. Although the PSF does not vary drastically 
with inclination angle, it is a strong function of photon energy. The energy- dependent 
half-angle Oqi of a cone containing ~ 67% of the 7-rays from a source is well fit by 
the function 

^67 = 5?85 {E/m MeV)-°-^=^^ , (2.1) 

where ii^ is the energy in MeV ( [Thompson et al. 1993| ). 

The energy dispersion files describe the energy resolution of EGRET. For incli- 
nation angles less than 20° and energies greater than 100 MeV, the energy resolution 
is ~ 20%. For larger inclination angles or lower energies, 7-rays are much more likely 
to produce pairs for which one of the charged particles misses the TASC, leading to 
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a greater degree of uncertainty associated with the energy estimate. 



2.4 Instrument Operation and Observations 

Upon deployment on 1991 April 7, CGRO was put into an nearly circular orbit of 
455 km, with an orbital period of ~ 93 minutes. The optimal functioning altitude is 
between 350 and 450 km, which avoids excessive drag while remaining below most of 
particle radiation in the South Atlantic Anomaly. Naturally, a gradual deterioration 
of the orbital radius will occur, so CGRO is equipped with an on-board propulsion 
system and carries sufficient fuel to reboost the spacecraft altitude. By 1993 October, 
CGRO had declined to a orbital radius of 345 km and the reboost procedure was 
initiated. The reboost was successfully completed in 1993 December, returning the 
spacecraft to the original orbital radius of 450 km. This procedure can be repeated 
as necessary in the future. 

Due to the low flux of high-energy 7-rays, a typical EGRET observation will 
last from one to three weeks. During this time, the pointing axis of the instrument 
is stable to within 0?5. The CGRO mission has been divided into phases which 
focus on different mission objectives. The observation dates, pointing directions, and 
numerical designations for the flrst three phases of CGRO operation are listed in 



Tables |2.2a| -c. The flrst month of the mission was dedicated to instrument testing, 
including pointings to well-studied objects such as Crab and Vela. Six observations 
from this test phase are of sufficient quality that they are now grouped with the 
Phase I observations. The initial fifteen months of observations that comprise Phase 
I were devoted to carrying out the first complete survey of the 7-ray sky ( [t'^ichtei 



et al. 19^ ). Phase II lasted from 1992 November until 1993 September, revisiting 



some of the more interesting regions found in Phase I, as well as observing some 



targets selected through the Compton Guest Investigator Program (|Bunner 1989| ). 
During Phase III, which began on 1993 August 17 and ended on 1994 October 4, 
more than 50% of the observing time was allocated to Guest Investigators. Due to 
scheduling constraints, the first observation in Phase III occurred between the final 
two observations of Phase II, explaining the gap between the initial observations of 
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Phase III. A record of all of the EGRET observations, as well as periods when the 
instrument was inactive, is kept in the instrument timeline file. 
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Table 2.2a: CGRO Phase I Observations 
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Table 2.2b: CGRO Phase II Observations 
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Table 2.2c: CGRO Phase III Observations 
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In order to minimize the amount of Earth albedo photons which might trigger the 
detector, EGRET has specific modes of operation to hmit its acceptance angles as 
the Earth enters its field of view. Each mode has a corresponding set of calibration 
files which are used to calculate instrument response for that particular mode. For 
each observation, an Exposure History (EXPHST) file is created to keep track of the 
times when the EGRET instrument changes mode and how much livetime is spent 
in that mode. 

In the initial stages of operation, the data from all four instruments was recorded 
on-board the spacecraft and then telemetered to the Earth via the Tracking and Data 
Relay SateUite (TDRS) once every other orbit. Unfortunately, the tape recorders 
failed in 1992 March and all subsequent data had to be telemetered in real time at 
32 kilobits per second during times of TDRS contact, reducing the EGRET data 
coverage to ~ 60%. With completion of an Australian TDRS ground station — named 
the GRO Remote Terminal System (GRTS)— starting in 1994 January, the EGRET 
real time data coverage increased to more than 80%. 

The gas in the EGRET spark chamber modules gradually deteriorates with time, 
causing some degradation in the instrument sensitivity. Additionally, lower energy 
photons have a reduced probability of being detected as the instrument ages, produc- 
ing a drift in the instrument spectral response. This causes sources to be measured 
with harder spectra over time. These effects are corrected for when calculating the in- 
strument response functions by applying energy- dependent scaling factors determined 
from the long-term monitoring of the photon detection rate in five energy ranges: > 
100, 30-100, 100-300, 300-1000, and > 1000 MeV. When the instrument perfor- 
mance degrades below a designated level, the gas in the spark chamber is replaced. 
An on-board gas replenishment system is capable of refilling the spark chamber vol- 
ume five times over the lifetime of the instrument. Gas refills have occurred on 1991 
December 2-3, 1992 December 3-4, 1994 February 8-9, and 1994 November 1-2. 

Flight data have indicated that below 70 MeV, the instrument response differs 
slightly from the behavior determined by the pre-flight calibrations, which were known 
to have limitations at lower energies. In particular, two effects were noted by Thomp- 
son et al. (1993): (1) the back-scattered calibration beam, assumed to consist of 7-rays 
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produced by inverse- Compton scattering, had a bremsstrahlung component extending 
out to calibration energies; and (2) the cahbration measurements were smoothed over 
a broad energy range. In turn, the effective EGRET sensitive area at lower energies 
was over-estimated. Based primarily on the assumption that the Crab pulsed photon 
spectrum has an essentially continuous power-law form over many decades of energy 
( [Strong et al. 1993D , it was determined that the effective area between 30-50 MeV was 
too high by a factor of 2.7 ±0.7, and the 50-70 MeV effective area was overestimated 
by a factor of 1.4 ± 0.2. These factors have been incorporated into the calculation of 
the instrument response functions. 

Secondary pairs formed in the spark chamber are more likely to be detected and 
accepted for analysis if their paths are closer to the instrument axis as opposed to 
those with larger inclination angles. This introduces a small distortion in the derived 
source locations, systematically shifting the determined position toward the EGRET 
pointing axis. Flight data has shown this effect is only consequential for aspect angles 
greater than 20° from the instrument pointing axis ( [i'hompson et al. 1993| ). Based 
on an empirical fit to the deviations as a function of aspect angle, the data have been 
corrected to account for this so-called fisheye effect. 

2.5 General Data Processing 
2.5.1 Standard Data Files 

The events that are determined to be 7-rays are stored in a Primary Database file 
corresponding to that observation. These files contain all of the measured and derived 
quantities for each event. Since these files are too large to be manipulated easily, 
pertinent information such as photon energy, time of arrival, and direction of incidence 
are stored in Summary Database (SMDB) files. The SMDB files are distributed to 
the participating institutions and Guest Investigators for analysis. 

To facilitate the spatial analysis of the data, the photons in the SMDB files are 
binned in 0?5 x 0?5 pixels for the following broad energy ranges: 30-100, 100-300, 
300-1000, > 100, > 300, and > 1000 MeV. The photons are also binned for the 
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following ten standard energy ranges: 30-50 MeV, 50-70 MeV, 70-100 MeV, 100- 
150 MeV, 150-300 MeV, 300-500 MeV, 500-1000 MeV, 1-2 GeV, 2-4 GeV, and 
4-10 GeV. These maps are stored in standard FITS (Flexible Image and Transport 
System) format as described by Wells, Greisen, & Harten (1981). The pixel sizes are 
chosen to obtain a reasonable count rate in each bin while avoiding excessive loss of 
spatial information. 

For each of these photon maps, an analogous 0?5 x 0?5 map of EGRET exposure 
is created. If a detector has exposure £^ to a source with photon flux F, then the 
number of counts which will be measured is 

N = E£ , (2.2) 

with F in units of photons per unit area per unit time. In practice, a detector will 
only detect photons over a certain energy range, so the flux F in equation ( |2.2| ) should 
be defined as 

F{AE) = [ I{E)dE , (2.3) 

JAE 

where AE is the energy range being considered and I{E) is the differential flux as 
a function of energy in units of photons per unit area per unit time per unit energy. 
The differential number of counts that will be detected by EGRET from a source of 
intensity I{E) is 

dN = I{E)A{E)dEdt , (2.4) 

where A{E) is the energy-dependent effective area of the instrument, and dt is the 
differential unit of time. However, due to energy dispersion, these counts will not 
necessarily be measured with their true energy E. Taking into account the energy 
dispersion of the instrument, the correct expression for the number of counts that 
will be measured in the energy range AE from a source of differential flux I{E) is 

N{AE;9,(j),m) =T{9,(j),m) dE' dE I{E) A{E;9,^,m) R{E ^ E';9,^,m) , 

JAE Jo 

(2.5) 

where T{9, 0, m) is the amount of instrument livetime spent observing a source at 
an aspect angle 9 and azimuth while in viewing mode m; A{E] 9, 0, m) is the 
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corresponding effective area at energy E; and R{E E'; 6, 0, m) is the probability 
per unit energy that a photon of true energy E will be measured with an energy 
E'. Substituting equation ( p.3|) and equation (|2.5|) into equation ( pl2|) and solving for 
exposure shows 

S{AE; e,(j)) = Y, T{e, 0, m) A{AE; 6, 0, m) , (2.6) 

m 

with 

A(AE; ^, 0, m) = J^^dEm ' ^ ^ ^ 

Here, 74(Ai?; 9, 0, m) is the average effective area of EGRET over energy range AE 
at an inclination angle 9, azimuth 0, and viewing mode m. Typically, it is assumed 
that the source photon spectrum behaves as a power law, such that 

I{E) = Jo photons cm"^ s"^ MeV"^ , (2.8) 

where a is the spectral index, which is typically ~ —2.0. 

The average exposure in each 0?5 x 0?5 bin is evaluated using equation (|2.6|) . The 
total observing time spent in each mode T{9, 0, m) is found by reading the EXPHST 
file and determining when a source at 9 and is occulted by the Earth. The effective 
area of equation { \^.7\) is calculated via interpolation of the results in the calibration 
files ( Fierro 19941) . 

Once the spatial maps of photon counts and exposure have been made for an 
observation, they can be added directly to results from other observations to produce 
cumulative counts and exposure maps. In this way, the photon intensity maps can 
be obtained for any set of observations by summing their corresponding counts and 
exposure maps, and then dividing the total number of photons in a bin by the total 
average exposure for that bin. The counts, exposure, and intensity maps above 100 
MeV have been summed for all of the observations in Phase I, II, and III, and the 
galactic coordinate Aitoff projections are shown in Figure ^.2[ The three brightest 



sources along the galactic plane are Vela, Geminga, and Crab at galactic longitudes 
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Photon Counts 




Instrument Exposure 




7-ray Intensity 

Figure 2.2: EGRET all-sky maps of photon counts, instrument exposure, and in- 
tensity. These galactic coordinate Aitoff projections were formed from all EGRET 
data above 100 MeV accumulated during Phases I, II, and III. The intensity plot in 
the lower panel is obtained by dividing the 7-ray counts in the upper panel by the 
EGRET exposure in the center panel. 
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of 263?6. 195? 1, and 184?6, respectively. It should be noted from the exposure map 
that instrument exposure is highly non-uniform, with particular concentration on the 
galactic center, the anti-center, and the Virgo region. 



Chapter 3 



Instrumental Background 



Absolute measurements of isotropic backgrounds of any kind are always very difficult 
for the simple reason that there is often no handle with which to distinguish back- 
grounds which constitute real signal from other instrumental noise. The revolutionary 
discovery of the Cosmic Microwave Background (CMB) is an example of a measure- 
ment that was possible only because the experimenters had great confidence that 
they knew the noise characteristics of their detector; in that case a microwave horn. 
The measurement of the IDGRB in the MeV range made using balloon measurements 
( P?rombka et al. 1977|) showed evidence of an 'MeV bump'. However, measurements 
in this energy range are complicated by the prevalence of nuclear activity. A detailed 
discussion of backgrounds in this energy range are discussed in ( [Gehrels 1992| ). More 
recent measurements made aboard CGRO have shown that these early measurements 
contained significant instrumental backgrounds which are likely responsible for the 
'bump' spectral feature ([Kappadath et al. 1995|) . The modulation of this background 
in the earth's magnetic field was of critical importance in rejecting instrumental back- 
grounds. This and other techniques described in the following discussion are used to 
evaluate EGRET's instrumental background. 



3.1 Sources of Instrumental Background 
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3.1.1 Charged Particles 

The earth's atmosphere at sea level provides ~ 20 attenuation lengths of shielding to 
photons with E > lOOMeV. For this simple reason it has always been necessary to get 
above the atmosphere in order to do 7-ray astronomy. The challenge of this approach 
is that the atmosphere also attenuates cosmic rays (~ 27 radiation lengths at sea 
level) so operating a 7-ray telescope in orbit necessitates the rejection of this cosmic 
ray flux which can create showers in a pair conversion telescope such as EGRET. 

The integrated flux of charged particles in the terrestrial vicinity is well approxi- 
mated by, 

j(> E) = lO'^^-i-^^cm-^s-isr^i . (3.1) 

For typical points in the EGRET orbit, the relevant cosmic ray fluxes are on average 
~ 10^ times higher than the corresponding 7-ray fluxes. One must carefully ensure 
that these particles do not contaminate the data. 
a) Downward Moving Gosmic Rays 

As outlined in the previous chapter, the principal subsystem for the rejection 
of downward moving charged particles is the Anti-Coincidence Dome (AC dome). 
Any particle traveling through the spark chamber to the calorimeter must pass first 
through a 2 cm layer of plastic scintillator. The majority of the charged particles will 
pass through as minimum ionizing particles. The energy deposited in this layer is 
very efficiently converted to light which is collected by an array of photo-multiplier 
tubes arranged at the base of the dome. The probability that a charged particle 
passing through the AC dome would escape detection was measured before launch to 
be < 10~^ ([Kanbach et al. 1974] ) 

A more difficult background to reject is caused by the material in front of the 
AC dome. Such material can initiate hadronic showers when bombarded by charged 
particles. A common product in such showers are neutral vr'^'s which can then decay 
into 7-rays which yield triggers indistinguishable from those caused by cosmic 7-rays. 
For this reason, care was taken to minimize the amount of material in front of the 
AC dome in the design of EGRET . It was necessary, however, to cover the AC dome 
with a thermal blanket, which served to stabilize the AC dome temperature, provide 
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a light shield, as well as to shield against micro-meteoroids. The resultant blanket 
has a thickness of ~ 3 cm and a mass per unit area of 0.17g cm~^. 
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Figure 3.1: Example of a proton event interacting in the thermal blanket to yield a 
spurious 7-ray event. 

The background caused by this unavoidable foreground material was measured 

at Brookhaven National Laboratory in 1987. A proton beam was directed at the 
instrument at various incident angles and the resultant count rate was measured. 
These measurements were hampered by the high ambient radiation levels in the test 
beam but an upper limit on this background was estimated to be 



jproton{E > lOOMeV) < 1.3 X IQ-'^cm-^s^'sr"' . 



(3.2) 
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b) Upward Moving Cosmic Rays 

There are additional pathways for cosmic rays to enter the spark chamber that 
must be guarded against. The AC dome shields the spark chamber only against 
downward moving cosmic rays. It is quite likely that a cosmic ray can enter the 
telescope from the calorimeter side and travel upward through the detector. If such 
a particle were to continue out the front side of the detector it would most likely 
trigger the AC dome and thus not be problematic. However, it is also possible that 
such a particle could undergo a hadronic shower and have the products range out in 
the detector. Such an event could conceivably yield a pattern of tracks that could be 
mistaken for a downward traveling 7-ray event. 

6000 1 • •: • I • • • I • • • I • • • I • • • I • • • I 




TOF (ns) 



Figure 3.2: The measured TOF for 17,000 flight events. Triggers are not generated 
for upward moving events with TOF>0. In the limit of infinite time resolution there 

should be no events with -2 ns<T0F<2 ns because these events measure speeds 
greater than the velocity of light. The events in this range arc indications that the 
TOF measurement time resolution is worse than predicted. The spike at T0F=-5 
represents the overflow events which are poorly understood. 

This source of background has been eliminated by the inclusion of a time of 
flight scintillator system between the upper spark chamber and the calorimeter. This 
system is comprised of two arrays of plastic scintillator separated by 60 cm. The TOF 
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of any charged particles traveling in the detector can thus be determined allowing 
the instrument to distinguish between upward and downward moving showers. The 
drawback of this approach is that this system severely limits the field of view because 
it requires a sufficient longitudinal length to allow a TOF measurement of relativistic 
particles. The resultant geometric aspect ratio is correspondingly elongated leading 
to a reduced field of view. 

Figure |3.2| shows the distribution of TOF values for a set of 17,000 fiight triggers. 
The TOF value is measured and stored as a 6 bit word which is telemetered to the 
ground along with the event record. Because events with TOF>0 do not generate 
fiight triggers, the distribution is truncated at this point. The distribution shows 
a peak at TOF=-2.5 ns. The separation of 60 cm between the two scintillating 
layers results in a T0F=2 ns for a relativistic particle traveling vertically through the 
chamber. Slightly larger TOF values correspond to relativistic particles traveling at 
an angle through the TOF system. The events in this peak thus constitute downward 
traveling 7-ray events. There should be a sharp cutoff in the distribution at TOF=- 
2ns. The spread of the observed distribution to smaller TOF values is a measure of 
the time resolution of the TOF system. The flight data indicate a time resolution of 
~ Ins. The pre-flight optimization of the TOF system ([Hunter 1991|) indicated that 
the time resolution of this system was ~ 0.45 ns. This value was obtained after careful 
optimization of cable lengths. It is possible that the light travel times have slightly 
altered over time which has slightly degraded this resolution. This degradation can 
lead to contamination by allowing the tail of the upwardly moving events to trigger 
the telescope with TOF values smaller than 0. There is some direct indication that 
this is happening in that the distribution is seen to be increasing near TOF=0. These 
stray triggers could constitute a signiflcant background and thus must be rejected on 
the basis of their shower patterns in the spark chamber as will be discussed later. 

It should also be pointed out that there are a large number of EGRET triggers for 
which the TOF measurement has overflowed. These triggers are generated by hits in 
the two layers separated by more than ±5 ns. It is not well understood what causes 
these triggers that are rejected on the basis of these anomalous TOF values. A large 
fraction of them could be due to accidental coincidences between the two scintillator 
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planes. These planes have a trigger rate in flight of ~ 2 — 3 kHz. Many of these hits 
will be in coincidence with a veto layer hit however and as such will not generate a 
trigger. 



3.1.2 Neutron Background 

Cosmic ray spallation in the earth's atmosphere produces a background of secondary 
particles. Among the products of such interactions are neutrons. Furthermore, the 
active sun can serve as an additional source of neutrons in low earth orbit. These 
neutrons are problematic in that they are not vetoed by the AC dome. They may, 
however, interact in the spark chamber and produce downward traveling showers. 

The fluxes of these secondary particles are not well measured. Gehrels ( Gehrels 
1992D presents a calculated spectrum in the EGRET energy range. We can use this 



spectrum to get a rough estimate of the number of such events to be expected. His 
spectrum must be viewed as a very conservative upper limit on the neutron flux be- 
cause his calculation includes neutrons produced as secondaries from charged particles 
interacting with a particular spacecraft. These events will be vetoed by EGRET. The 
spectrum given yields .3 neutrons cm~^s~^. If we assume the tantalum nuclei in the 
tracker conversion foils convert these neutrons inelastically causing hadronic showers, 
then the probability of conversion of a neutron into a shower passing through a half 
radiation length of tantalum is, 

P = ^ipA^^rl = ( 6^10'' g-i) X (.2 cm) x (16.6 g cm-3)(180)2/3l0-26cm2 = .02 , 

where we have used the atomic weight of tantalum (A) and its density (p). If we use 
the physical area of the detector as well and assume a live-time fraction of 50% and a 
solid angle of 40° this results in a trigger rate of 4 Hz. The rate seen by EGRET will 
be much smaller than this because the AC dome will veto any cosmic ray shower that 
produces neutrons as well as charged particles when interacting with the spacecraft. 
Nevertheless, there may be some neutron induced hadronic cascades produced in the 
tracker. These events must be rejected on the basis of their shower patterns in the 
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spark chamber. 

3.1.3 Albedo 7- Rays 

The EGRET orbital altitude of 350-450 km ensures that downward traveling 7-rays 
pass unattenuatcd through the residual atmosphere above the telescope. Further- 
more, downward traveling cosmic rays have a similarly small chance of producing 
hadronic showers which could produce a spurious 7-ray signal. However, the amount 
of matter in front of the telescope increases dramatically when the telescope is di- 
rected near to the earth's limb. In these direction, cosmic ray interactions with the 
earth's atmosphere generate copious 7-rays. 

While these albedo 7-rays cause triggers identical to their celestial counterparts, 
they are also imaged in the same way as the signal 7-rays. These spurious events can 
then be rejected simply on the basis of their proximity to the earth's limb. In fact, 
EGRET is capable of operating in several directional modes in which combinations 
of the trigger scintillator array tiles are disabled and as a result events from near the 
earth's limb are not permitted to trigger the telescope. In this way, excessive triggers 
caused by albedo 7-rays are prevented as these events would be eventually rejected 
based on their origin near the earth's limb. 

3.2 Triggering Modes 

The EGRET instrument triggers and generates an event record when the following 
criteria are met: 

• A signal is detected in both the upper and lower arrays of scintillator in the 
TOF system described above. 

• The TOF between these two signals is within an acceptable range. 

• The energy deposited in the TASC exceeds a defined threshold value (nominally 
6MeV). 

• No signal from the AC dome is detected. 
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The trigger rate in this mode of operation in orbit is ~ 0.7 Hz. If one disregards 
the problematic TOF overflow events this rate drops to ~ 0.4 Hz. We can estimate 
the trigger rate caused by 7-rays (both celestial and the more dominant albedo) as 
follows. 

T ■ 1 1 

Rate = VAO(^) ~ 20xl0-^x-x(80cm)^x(27r(l-cos(40°))x- ~ .4 Hz , (3.4) 

for £^ > 30 MeV and we have assumed a photon conversion rate of 33% and used the 
approximate physical area of the tracker. While this estimate is not accurate to within 
a factor 2 it does indicate that the hardware trigger has been successful in reducing 
the trigger rate from the several kHz cosmic ray flux (typical TOF deck rate = 2 
kHz) down to values less than 1 Hz of which 7-rays constitute a large fraction. This 
analysis is not sufficient to determine whether there remains any cosmic ray signal. 
In fact the distribution of TOF values indicates that there is still a substantial cosmic 
ray contamination in the hardware triggered events. Furthermore there may be a 
significant trigger rate caused by secondary neutrons. These must be rejected by the 
characterization of the track patterns in software. 

3.2.1 Search and Analysis of Gammaray Events (SAGE) 

SAS 2 and COS B relied on a tedious cvent-by-event interactive analysis to select 
desired events and determine their structure. With the significant increase in data 
recorded by EGRET , this type of manual analysis becomes impractical, and an au- 
tomated procedure must be implemented. 

For EGRET, analysis is done using a computer program called SAGE (Search 
and Analysis of Gammaray Events). For every event in the EGRET telemetry, the 
spark chamber tracks are analyzed by SAGE in a nine step process. First, the event 
must have generated a given minimum number of sparks in both orthogonal views. 
Second, three sparks must be found within a specified distance of each other, and 
must define a reasonable path. These three spark lines are known as triplets. In 
the third step, triplets which start in the top conversion plate or trace back to the 
instrument walls are rejected. In addition, at least one triplet must form before the 
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upper scintillator array. In the fourth step, all possible tracks from the remaining 
triplets are constructed, with duplicate tracks eliminated. The fifth step selects the 
electron and positron tracks from the candidates calculated in step four. Step six 
checks to see if there is a better choice of vertex between the two paths. In the 
seventh step, the two sets of tracks in the two orthogonal views are correlated. The 
eighth step determines event descriptors such as the number of spurious sparks and 
the curvature of the track. In the ninth step, the events are stored. 



Table 3.1: Summary of SAGE Rejected Events 



Total Triggers (06/03/91-06/04/91) 


55439 


Rejection Condition: 




Too few sparks in given projection 


7022 


TOF=0 


20857 


No triplets in given projection 


1666 


Conversion in first plane 


2150 


Conversion below scintillator plane 


4557 


Wall event 


7827 


Angle exceeds 45° 





Track pattern bad 





Fails to meet minimum point requirement 


5800 


Single rejected events 


760 


Total Rejected E\'eiits 


50639 


Total Accepted Events 


4800 



The surviving candidate events are divided into two groups, definite 7-ray events 
and questionable events. Extensive studies have shown that the events which are 
rejected by SAGE are typically either not 7-ray events or are unsuitable for analysis. 
Roughly one good 7-ray in 700 is rejected by SAGE. Subsets of the events which are 
classified as definite 7-ray events are routinely monitored by trained analysts whose 
determinations agree with the software to better than one part in a thousand. About 
12-15% of all events are initially identified as questionable. The analysts manually 
review these events by examining the event tracks on a graphical display. These 
events are usually either accepted or rejected as structured by SAGE. Approximately 
0.5% of the events are interactively restructured. As a continuing check of the SAGE 



CHAPTER 3. INSTRUMENTAL BACKGROUND 



40 



program, 1% of the EGRET data is randomly selected for graphical examination by 
the analysts. 

Table 1 shows a summary of the cuts made by SAGE on one days worth of data. 
An overall rejection factor of 0.08 results. While many of these rejections represent 
7-rays whose hit patterns were simply unsuitable for analysis (e.g. rejections for too 
few sparks) others indicate likely hadronic shower events (e.g. fails to meet minimum 
track requirement). The selectivity of this process ultimately determines the residual 
cosmic ray induced background. 



3.3 In Flight Determination of Background 

The various sources that contribute to EGRET^s instrumental background have been 
outlined above as have the systems designed to reject them. We now turn to the task 
of using the flight data to determine the efficiency with which backgrounds are being 
suppressed. 



3.3.1 Albedo Rejection 

As mentioned previously, albedo 7-rays from the earth's limb are rejected on the basis 
of the proximity of their reconstructed direction to the earth's limb. This is achieved 
by making a zenith angle cut on the reconstructed 7-ray arrival angle. As shown in 
Figure |3.3| , the zenith angle (a) is the angle between a photon arrival direction and 
the vector from the earth's center to the spacecraft. The earth's limb appears at a 
zenith angle of a ~ 111° in the EGRET orbit. For all zenith angles greater than this 
value, the earth is directly within the field of view. Because the earth's atmosphere is 
illuminated from above by cosmic rays, the limb produces the greatest 7-ray intensity. 
Cosmic rays that interact on the far side of the earth's limb produce showers that 



are directed into the telescopes field of view producing bright hard 7-rays (P?hompson 
et al. 198T|) . The limb spectrum is therefore brighter and harder than the radiation 
from the earth's disk. 

In the limit of infinitely sharp point spread functions, all photons arriving from 
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450 km 



6400 km 



Figure 3.3: Orbital geometry showing the zenith angle to the earth's limb. 



zenith angles less than a = 111° should be celestial in origin. However, the EGRET 
point spread function is several degrees wide. Moreover, the point spread functions 
are known to have tails which result in a finite probability of measuring a photon 
arrival tens of degrees away from the true origin. This can occur if only one of the 
e~^/e~ pair tracks are reconstructed. It is this phenomenon which leads to the so 
called 'fish-eye' effect described in Chapter 2. While the correction for the 'fish-eye' 
effect does lead to correct point source positions, it does not eliminate the possibihty 
of photons from outside the field of view leaking in. The result of this effect is that 
bright features near the edge of the field of view can leak into the data stream. This 
is particularly true of the earth's limb which is the brightest 7-ray feature in the sky. 



Figure below shows the distribution of measured 7-rays in zenith angle, summed 
over many pointings. For small zenith angles, the photon count rate is roughly con- 
stant. This is because these photons are dominated by celestial sources which are 
distributed randomly in zenith angle. As a result these features tend to average out. 
Near the earth's limb the count rate drops. This is because when EGRET is directed 
near the earth's limb, restricted field directional modes are incorporated into the 
hardware trigger so as to reduce the overall trigger rate. This process is not perfectly 
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efficient and there is a resulting exclusion of photons that arrive from directions near 
the limb. The bright feature clearly visible at a = 111° is caused by albedo pho- 
tons from the earth's limb. Despite the directional modes intended to reduce the 
albedo trigger, these albedo photons are a factor 3.5 times as common as celestial 
7-rays. This is a strong background that must be carefully eliminated. The standard 
energy-range dependent zenith angle cuts used in EGRET data analysis are shown 
superimposed on this figure. It is clear that albedo photons contribute to some extent 
inside the adopted cuts. While this does not affect analysis of bright sky features it 
can have a potentially important impact on analysis of the faint high latitude diffuse 
emission. 

104 



O 

" 10^ 



50 100 150 200 

Zenith Angle (degrees) 

Figure 3.4: Distribution of photons in zenith angle. The standard acceptance cuts 
are shown as a vertical line for each energy range. 

In order to determine the appropriate zenith angle cut, the flight data were ana- 
lyzed using a variety of zenith angle cuts. Exposure and counts maps were generated 
in four independent zenith angle ranges and in two energy ranges {E > 100 MeV, 
and E > 1000 MeV). These maps were constructed in celestial coordinates. The 
intensity was then calculated in a set of high latitude bins (Qy = 10° x 10°). 

In order to calculate the albedo contamination, several assumptions were made. 
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First, the intrinsic celestial intensity in each bin is assumed to be constant. This 
implies that measurements of intensity that are not concurrent should still give the 
same result. This is known to be untrue of high latitude point sources which can vary 
on the time scale of days ( |von Montigny et al. 19951) . As a result, regions within 10° 
of a strong point source were excluded from this study. The sources excluded were 
those found in the combined data from Phase I+II. 

Secondly, there was assumed to be negligible albedo contamination in the zenith 
angle range a < 50°. The albedo contamination in each spatial bin, Izij, is thus 
calculated to be the difference between the intensity measured in a particular zenith 
range Z and the intensity calculated using a < 50°. The data will be shown to be self 
consistent with this assumption. The uncertainty in each measurement, cxjj, is esti- 
mated from the counting statistics which allows the weighted mean to be calculated 
over all spatial bins, 

. (3.5) 



v.. ^ 

The uncertainty associated with this value is given by 

1 



aj^ = . (3.6) 



Table 3.2 summarizes the results. There is a clear signature of albedo contamina- 
tion outside zenith angles of 65°. While one could argue that this should determine 
the best value of the albedo cut, it should be noted that fully 38% of the standard 
exposure occurs in the zenith angle range 65° < a < 95.375. As a result, such a 
stringent zenith angle cut eliminates a large fraction of the data. 

A better strategy is to require some minimal signal to noise ratio. At high latitudes 
the diffuse intensity for E > 100 MeV is ~ 1 x 10^^ which sets the standard for 
measurement accuracy. This should be compared with the albedo contamination 
averaged over all included zenith angles. These values are shown in Table 3.2 along 
with estimated signal to noise ratios at high latitudes. The standard EGRET zenith 
angle cut in the energy range E > 100 MeV of amax = 95.375° results in an albedo 



contamination of /aibedo = 1-5 x 10 which is ~ 15% of the high latitude diffuse. A 
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Table 3.2: Summary of Albedo Contamination Results 





E > lOOMeV 


Zenith Angle range 


a < 50.0 


50.0 <a< 65.0 


65.0 < a < 80.0 


80.0 <a< 95.375 


Fractional Exposure 


0.38 


0.24 


0.24 


0.14 


Differential Albedo 
Intensity (xlO^^) 


0.0 


0.4 ±0.4 


1.9 ±0.4 


6.8 ±0.6 


Mean Albedo Intensity 
(a < a^„a;)(xl0-6) 


0.0 


0.15 ±0.15 


0.64 ± 0.2 


1.5 ±0.2 


Signal:Noise 


N/A 


66 


16 


7 




E > lOOOMeF 


Zenith Angle range 


a < 65.0 


65.0 < a < 80.0 


80.0 < a < 95.0 


95.0 < a < 105 


Fractional Exposure 


.60 


.25 


.13 


.02 


Differential Albedo 





-1.5 ±1.0 


3.4 ±1.5 


48.8 ± 6.5 


Intensity (xlO~'^) 










Mean Albedo Intensity 
{a < a^„a;)(xl0-6) 


0.0 


0.0 ±0.3 


0.45 ± 0.32 


1.41 ±0.3 


SignahNoise 


N/A 


> 33 


22 


7 



cut of ajnax — 80.0° results in a mean albedo contamination of /albedo = 0.6 x 10~^ 
while only excluding 14% of the data. This is now less than 10% of the high latitude 
diffuse which is acceptable. It is important to note that on the galactic plane, the 
galactic diffuse emission is at least an order of magnitude greater and as a result the 
standard zenith angle cuts are acceptable. 

Despite the fact that the earth is not locahzed in the sky, the albedo contamination 
it produces is not isotropic. The spacecraft orbital pole in a low eccentricity orbit 
remains at constant zenith angle near to the earth's limb throughout a pointing 
period. As a result, pointing near the orbital pole is avoided whenever possible. 
However, in order to obtain coverage of the full sky it is necessary to point near 
the orbital pole on occasion. The inclination of the EGRET orbit is 28.5° and it 
precesses about the celestial pole. As a result regions near the celestial poles and 
in particular certain regions with declinations near 5 — ±61.5° tend to be viewed at 
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Figure 3.5: Residual intensity in various zenith angle ranges. The vertical axis shows 
the mean additional intensity in a set of high latitude pixels in a given zenith angle 
range compared to the intensity in the reference zenith range: zenith< 50°. The 
horizontal error bar indicates the extent of the particular zenith ranges and the vertical 
error bars are statistical uncertainties. 



higher average zenith angle than regions near the orbital/celestial equator. Figure 3.6 
shows the ratio of the exposure in the rigidity range 80 < a < 95.375 to exposure 
in the zenith range a < 80.0. Note the coverage of the regions near the celestial 
pole occurs more frequently at higher zenith angles than regions near the celestial 
equator. The result is to distribute the albedo contamination preferentially in these 
regions of higher average zenith angle. These regions must be treated carefully to 
avoid confusion between albedo induced features and real celestial signals. 



3.3.2 Rigidity Modulation 

The earth's magnetic field provides the handle that allows us to gain some under- 
standing of the level of charged particle contamination. If we follow a charged cosmic 
ray particle in from infinity, the behavior of this particle when it interacts with the 
earth's magnetic field depends on the particles momentum. Low energy particles 
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Right A&cension 

Figure 3.6: Ratio of exposure in the zenith range 80 < a < 95.375 to the zenith range 
a < 80.0. Regions that show up hght colored in this plot are dominated by high zenith 
angle exposure and as a result contain higher than average albedo contamination 

travel in tight orbits around field lines and drift toward the magnetic poles where 
they are funneled into the trapped radiation belts. Higher energy particles are able 
to penetrate all the way to the earth's atmosphere with little defiection. At a given 
point in the earth's magnetic field away from areas of trapped radiation, the effect is 
to introduce a low energy cut-off in the the observed spectrum of charged particles. 
This cut-off is traditionally expressed in terms of magnetic rigidity, 

^ ^ PC M cos^L 

^'fd + cos^cos^D^/^lf 



where p, z are the particles momentum and charge and M is the earth's dipole mo- 
ment. Particles entering the earth's magnetic field at an angle of 6 to the tangent to 
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the latitude circle, L will not penetrate to R, the geocentric radius, and are assumed 
to be trapped if their rigidity is less than that given by the above formula. 

The end result of this magnetic interaction is that the cosmic ray flux is screened 



to a greater extent in regions of high magnetic field. Figure shows the one day long 
EGRET time history showing both the rigidity variation as well as the AC dome rates 
which are telemetered to the ground as part of the housekeeping data at regular time 
intervals. The anti-correlation between rigidity and AC dome rate is clearly visible. 
The AC dome rate increases by a factor of four in the orbital regions of lowest magnetic 
field. This modulation can be used to check for a modulation in observed gamma ray 
background rate with orbital rigidity to measure any cosmic ray background leaking 
into the gamma ray stream. It is important to note that the AC dome rates do not 
directly measure cosmic ray rates because they include a significant contribution from 
hard X-rays. The well measured hard X-ray background was calculated to result in a 
AC dome trigger rate of ~ 20kHz which is entirely consistent with the observed AC 
dome trigger rates. 

The EGRET timeline was divided into two sets of intervals: periods when the 
rigidity was >14G'V" and periods when the rigidity was <1AGV. The resultant in- 
tervals are ~ 40 minutes long on average. The two rigidity ranges contain roughly 
equal numbers of 7-ray events. As in the investigation of albedo contamination, ex- 
posure and counts maps were constructed in celestial coordinates from the two sets 
of rigidity intervals. This results in a measured intensity for each point on the sky 
in times of high cosmic ray flux and in times of less cosmic ray bombardment. The 
difference is a measure of the instrumental background which is not rejected by the 
various systems. Again it must be assumed that the 7-ray signal is not variable on 
the time scale of a viewing period in regions away from strong point sources. Because 
the albedo intensity is also modulated by rigidity, a zenith angle cut of a < 50° was 
adopted. 

Figure |3.8| shows the distribution of the cosmic ray modulated intensity in a set 
of 10° X 10° pixels at high latitude. The weighted mean is calculated in a manner 
exactly analogous to the technique used in calculating the albedo contamination. The 
resultant flux is calculated to be Icr = —5 ± 4 x 10^^. This implies that there is no 
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Figure 3.7: Rigidity history for one day during VP 1.0. Note the modulation with the 
sateUite orbital period as well with the earth's rotational period. The gaps represent 
times when the instrument was off while passing through the South Atlantic Anomaly. 



measurable cosmic ray contamination in the data. The 2a upper limit on the residual 
intensity is Icr < 3 x 10~^. 

In order to extract the upper limit to the mean cosmic ray background from this 
measurement it is necessary to average the variable instrumental background over the 
history of the mission. 

(3.8) 



IIcR{t)dt 



T 



where Icr is the cosmic ray induced instrumental background and T is the total live- 
time of the mission. Because we cannot assume that the cosmic ray background is 
negligible in the high rigidity range, the differential intensity between rigidity ranges 
does not correspond to a measurement of the absolute cosmic ray background. Rather 
we must model the absolute comic ray flux and use the differential measurement to 
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Figure 3.8: Distribution of instrumental background measurements showing the mean 
and x^/i^. 

calculate an absolute intensity. If we assume the X-ray counting rate in the AC dome 
to be 14 kHz, then we can deduce IcRif) from the housekeeping data. Computing 
IcR in each rigidity range we find that, 

7c^<i4 ~ 2 X 7c^>i4 . (3.9) 

As a result. 



IcRo — IcR<li — IcR>li — '^IcR>U — IcR>lA — IcR>lA 



(3.10) 
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Furthermore, integrating over all viewing periods we find, 

T<i4~T>i4. (3.11) 
This in turn allows us to estimate the mean cosmic ray induced background to be, 

7 ICR<1A + ICR>14. '^^CRo + IcRo i cr r /'Q 1 o^ 

^CR = ^ = ^ = 1-5 X Jci^o ' l^.l/Ij 

where IcRq is the measured residual intensity between the two rigidity ranges. 

The resultant upper limit to the average cosmic ray induced instrumental back- 
ground is IcR < 4.5 X 10~^. This value indicates that the cosmic ray contamination 
is less than an order of magnitude less than the isotropic diffuse background and 
is negligible compared with the already documented albedo background. This limit 
is significantly more stringent than the pre-flight calibration of proton backgrounds. 
The ability of SAGE to reject the cosmic ray events which are not rejected by the 
hardware trigger has been demonstrated. 

While rigidity directly modulates the charged particle flux it modulates secondary 
neutron flux more weakly. A neutron produced at the earth's limb travels through the 
earth's magnetic field unimpeded until it interacts with the the EGRET instrument. 
As shown in Figure ^.3| , the interaction point can be 21° distant in latitude from 
the spacecraft latitude. For the EGRET orbit this implies that the interaction can 
occur at a point at which the rigidity is a factor 2 different from the value near the 
instrument. The net result is to slightly wash out the modulation of neutrons. Instead 
of a factor of 2 difference between the incident flux in the two rigidity ranges used, the 
averaged flux modulation is closer to 1.2 for typical points in the orbit. Nevertheless, 
the lack of a rigidity signal also constrains the neutron induced background. In light of 
the fact that charged particles that get through the TOF system are efficiently being 
rejected on the basis of their hit patterns in the tracker, it is reasonable to assume 
that the neutron showers are being rejected in a similar manner. That is, unless the 
neutron induced trigger rate is significantly larger than the rate of charged particles 
that escape the hardware trigger, the upper limit calculated above also applies to 
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neutron backgrounds. This is likely to be the case. 

3.4 Exposure Calibration 

Certain types of instrumental effects are most easily identified by looking at the data 
in instrumental coordinates. If a dead spot in the AC dome were to arise, this would 
show up very significantly as a bright spot in a particular area of the instrument. 
Furthermore, some inaccuracies in the calibration of the sensitivity of the instrument 
can be identified in flight as a non- uniformity across the instrument. 

Of course the difficulty that must be overcome in looking for this sort of signal is 
that the 7-ray intensity varies across the sky by more than two orders of magnitude 
from galactic center to galactic pole. The nonuniform intensity does not necessarily 
average out to yield a uniform calibration source that can be used to look for any 
instrumental asymmetries. Furthermore, there are many bright 7-ray sources that 
move around in instrumental coordinates from viewing period to viewing period. 

In order to be able to attempt an in flight calibration in instrument coordinates, 
it is necessary construct a rough model of the 7-ray sky. This model consists of 
three elements: a model of the galactic diffuse emission, a set of bright point sources, 
and some level of isotropic diffuse background. The galactic diffuse model will be 
described in detail in the next chapter. The point sources are extracted from the 
data using the Maximum Likelihood technique (see Appendix A), and the isotropic 
diffuse level is chosen to make the data match the model along the pointing axis. It is 
important to note that because a given point source is not localized in instrumental 
coordinates it is necessary to evaluate the strength of the source during each pointing 
in order to properly model its impact on the data. The general strategy is not to 
create a perfect model of the sky but rather to approximate a model of the sky. Any 
inaccuracies in this model should not be localized in instrumental coordinates and as a 
result should tend to become evenly distributed as one averages over many pointings. 
Instrumental effects such as inaccuracies in the exposure calibration are localized in 
this coordinate system and as a result tend to accumulate. 

It is necessary first to define a set of instrumental coordinates. This choice is 
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Figure 3.9: Schematic of the CGi? ^spacecraft showing the spacecraft axes and the 
instrument positions. 

dictated by the fact that we do not want a coordinate singularity in the center of the 



field of view. As shown in Figure p.lO| the coordinates are defined by, 

sin(^) =f-x, (3.13) 

and, 

where x, y, z are the spacecraft axes and f is the photon arrival direction. These 
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coordinates are exactly analogous to galactic coordinates with the spacecraft x-axis 
corresponding to the galactic pole and the i-axis corresponding to the galactic center 
direction. As a result (f), 6 will often be referred to as instrumental longitude and 
instrumental latitude respectively. It should be noted that the i-axis is given by 
{(t),9) — (0,0), the X-axis is given by (0,^) = (0,90), and the y-axis is given by 
{4),d) = (90,0). Also note that the COMPTEL instrument lies along the +x-axis 
while the — x-axis is directed toward the edge of the spacecraft. 

X (towards COMFTEL) 




Figure 3.10: Definition of the instrumental coordinates used in this section. They are 
chosen so as not to have a singularity along the pointing direction. 6, (j) are analogous 
to /, b with the x-axis representing the polar axis and the z-axis representing the 
galactic center axis. 

For each pointing, the coordinate transformation matrix is constructed to convert 
a vector in galactic or celestial coordinates into instrumental coordinates. This matrix 
is essentially the product of two rotation matrices. A map of exposure in instrumental 
coordinates is constructed from the standard map in galactic or celestial coordinates 
using the following technique. The center position of each bin in the instrumental 
map is transformed into galactic coordinates. The exposure per steradian at that 
point is calculated from the interpolation between the nearest bins in the celestial 
exposure map. The exposure is assumed to vary slowly across a .5° x .5° bin so 
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that the exposure is simply the bin size times this value. This same technique is 
employed in order to construct a map of sky intensity in instrumental coordinates. 
The maps are multiplied to produce a map of predicted counts in each instrumental 
bin which can be compared to the data which can be directly binned in instrumental 
coordinates. 



Figure p.llj and Figure |3.12| show the exposure and counts maps in instrumental 
coordinates for a set of 62 coadded viewing periods. The 'latitude' and 'longitude' 
angles are those defined in Figure p.lO| . The departure of the exposure map from 



azimuthal symmetry is caused by the embedded square structure of the spark chamber 
and calorimeter. The following cuts have been made on the photons in the SMDB: 

• Measured Energy E > 100 MeV. 

• Energy in TASC greater than threshold of 6 MeV. 

• Photon does not arrive during an excluded interval. 

• Zenith angle a < 80°. 

The total number of photons included in this study is 174,000. 

The measured intensity is shown in Figure |3.13| . The modeled intensity is shown 
in Figure p.l4| . In order to minimize the effect of the counting statistics, both maps 



are convolved with a Gaussian smoothing function with a 2° half-width. The viewing 
periods for this study were selected based on two criteria. All viewing periods within 
80° of the Galactic Center were excluded because the galactic model is known to be 
inaccurate in this region. Furthermore, no viewing periods were included in which 
the galactic ridge lay in the x — z or y — z plane. This was done so as to decouple any 
symmetries in the celestial features with the instrumental symmetries. Unfortunately, 
this eliminated many galactic viewing periods because due to the constraints of the 
OSSE] detector which points in the xz plane, the galactic ridge is often polarized in 
one of these two planes. 

The instrumental effect is measured as the fractional residual at each location in 
instrumental coordinates, 

i?.. = (3.15) 
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Figure 3.11: Exposure map in instrumental coordinates. The energy range is £^ > 
lOOMeV. 
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Figure 3.12: Measured counts in instrumental coordinates. The energy range is £^ > 
IQQMeV. 
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Figure 3.13: Measured intensity in instrumental coordinates. The energy range is 
E > lOOMeV. 
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Figure 3.14: Modelled intensity for the same set of viewing periods shown above. The 
model used is the Bertsch et al. model (see Chapter 4). A time variable set of point 
sources have been included as has an isotropic diffuse background. 
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where Cij, Mij are the measured and modelled counts in the bin ij. Figure 3.15 shows 
the normalized residual map for the data set shown above. There is clear evidence of 
features correlated with the instrumental structure. Figure |3.16| shows the normalized 
residuals for a Poisson sample of the model in order to indicate the level of statistical 
fluctuation present in this map. The features seen in the data clearly cannot be 
statistical fluctuations. 

There are two possible explanations for these features: background contamina- 
tion or exposure calibration inaccuracies. If the former explanation holds true, the 
excesses seen in the residuals are due to some sort of contamination which is particu- 
larly concentrated in certain areas of the field of view. Because this would imply that 
the contaminating intensity is uncorrelated to the 7-ray signal, the residual intensity 
should be independent of the 7-ray intensity of the sky being observed. This assump- 
tion can be tested by seeing whether the residual intensity at a particular position 
in the field of view correlates with the residual intensity at the same instrumental 
position measured using a very different set of pointings which give different intrinsic 
sky intensities at that point. On the other hand if the residual intensity is due to a 
miscalibration of the exposure, the excess intensity will depend on the intrinsic ce- 
lestial 7-ray intensity. In this case the normalized residual should correlate between 
different sets of pointings. 

This test was done using a set of high latitude pointings and a set of galactic 
plane pointings. The two data sets were completely independent and the modelled 
intensities were much higher for the galactic plane pointings. In each case the resid- 
ual intensities were calculated in instrumental coordinates as were the normalized 



residuals in instrumental coordinates. Figure 3.17 shows the instrumental residual 
intensity for the galactic plane pointings plotted against the same quantity from the 
high latitude map. The correlation is very poor. In contrast. Figure p.l8| , shows the 
correlation using the normalized residual maps. The correlation is much better. This 
test provides strong evidence in favor of the source of the residuals being an exposure 
miscalibration. 

The first feature of note is the pronounced hole along the ±y-axis and the +x-axis 
at inclination angles of > 30°. In order to understand this, the projection of the top 
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Figure 3.15: Normalized residual intensity in instrumental coordinates. The values 
plotted are defined in equation 3.15. A clear instrumental pattern is evident. The 
squarish profile shows the projection of the top of the spark chamber into these 
coordinates to emphasize the instrumental symmetry. 
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Figure 3.17: A correlation between the instrumental residual intensity maps for two 
independent data sets. The vertical axis represents the measured residual intensity 
in a set of galactic plane viewing periods. The horizontal axis represents the same 
values calculated using high latitude viewing periods. 



— 0.40 




Instrumenlal Normalized Residual (High Latitude) 

Figure 3.18: A correlation between normalized instrumental residual maps for two 
independent data sets. The vertical axis represents the normalized residual intensity 
in a set of galactic plane viewing periods. The horizontal axis represents the same 
values calculated using high latitude viewing periods. The correlation is clearly su- 
perior to that shown above indicating that the exposure calibration is more likely to 
be a source of the inaccuracies than contaminating radiation. 
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plane of the spark chamber into these coordinates is shown in Figure 3.15| as the outer 



squarish contour. It is evident that there is a correlation between the symmetry of the 
detector and the instrumental residuals. Furthermore, there are residual excesses that 
form a sqarish feature at inclinations greater than ~ 15% which is most pronounced 
along the lines which lie at 45° to the x and y axes. This effect is at worst ~ 20%. 

The prefiight calibration of EGRET was intended to be completed to yield accu- 
racies within ~ 1%. The source of the features that are apparent in Figure |3.15| are 
poorly understood. Some suggestions follow. 

The instrument was only calibrated in one octant. The calibration for the re- 
maining 7 octants was obtained through symmetry operations. There is an apparent 
variation from corner to corner in the residuals of about 10%. This could be due 
instrument asymmetries that were not detected. 

The most likely source of this calibration inaccuracy is the convolution of the 
model with the PSF. The correct calculation of the measured count rate as a function 
of instrumental position, C{9, 0), is as follows. 



where I{6, 0) represents the intrinsic sky intensity, S{6, 0) represents the instrumental 
exposure in a given direction, and PSF{6', 0'; 6, 0) is the appropriate point spread 
function. It is, tedious, however to carry out this convolution for every instrumental 
pointing and thus the analysis is expedited by calculating one convolved background 
model, 




(3.16) 




(3.17) 



leading to a modelled counts map. 




(3.18) 
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In the limit of a featureless background, the resultant discrepancy between mea- 
sured counts and modelled counts is given by, 

R{9, <P) = (3.19) 

In other words, the residuals map is given by the difference between the exposure 
map convolved with the PSF and the unconvolved exposure map. In the limit of a 
symmetric PSF, this does not introduce substantial error, however, the EGRET PSF 
is not symmetric because of the 'fish-eye' effect. The 'fish-eye' correction ensures that 
point sources are correctly positioned, but it does distort the background. 

Figure Figure |3.19| shows the residuals calculated using equation 3.16. For the 
purposes of this demonstration, a skewed Gaussian PSF is used. The features ap- 
parent in the flight data are reproduced using this method. It should in principal be 
possible to exactly reproduce these features given knowledge of how the PSF skew- 
ness varies across the instrument. For the purposes of this study, it is sufficient to 
correct the exposure maps in order to remove this instrumental effect. This technique 
is sometimes called 'flat-fielding'. It is important to note that this is only applicable 
for high latitude data for which the sky intensity is roughly uniform. 



3.5 Summary 

The in-fiight instrumental performance of EGRET has been studied in detail. Instru- 
mental background due to cosmic ray and neutron backgrounds in the energy range 
E > 100 MeV have been shown to be smaller than 4.5 x 10~^. This amounts to better 
than an order of magnitude signal to noise ratio for the detection of the IDGRB. This 
demonstrates the excellent selectivity with which the various EGRET hardware and 
software systems reject the high flux of cosmic rays in low earth orbit. 

The albedo 7-rays from the earth's limb are more problematic. Using an albedo 
cut of a < 80° minimizes this problem but still results in a background intensity 
above 100 MeV of 6 x 10"''. This background is concentrated in certain regions near 
the celestial pole requiring particular care when analyzing the diffuse background in 
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Figure 3.19: Residual introduced by convolution with a skewed Gaussian PSF. The 
plot represents a slice through the residuals map along the y-axis in instrumental 
coordinates. The features apparent in the measured data residuals are reproduced. 



these areas. 

Finally, the exposure calibration has been reviewed using the flight data. Signif- 
icant discrepancies between the flight data and the prefiight calibration have been 
revealed. These discrepancies are most hkely due to PSF skewness but can be cor- 
rected for through the 'flat-fielding' of the exposure maps. 



Chapter 4 

Diffuse 7 Radiation 



4.1 Galactic Diffuse Radiation 

The measurement of a diffuse source of galactic 7-rays dates back to the earhest days 
of 7-ray astronomy. A bright feature in a narrow band along the galactic plane was 
first seen by 0503 (P:^raushaar et al. 1972| ) as well as by a high altitude balloon 



instrument ( [b'ichtel et al. 1972| ). The subsequent instruments (S'AS' 2 and 005* B) 
were able to verify the correlation between 7-ray features and galactic structural 
features QFichtel et al. 1978| , [Strong et al. 1988|) . 



This radiation was quickly interpreted based on its spatial distribution and spec- 
trum as being associated with cosmic-ray/matter interactions and to a lesser extent 
cosmic-ray/photon interactions. The most compelling evidence for this origin of the 
galactic diffuse radiation was provided by the strong correlation between galactic 
matter tracers such as 21 cm radio emission and diffuse 7-ray emission ( [Fichtel et al 



1978| ) . It was also realized 7-rays could provide valuable information about the distri- 
bution of cosmic-rays in the galaxy because of the relative transparency of the galaxy 
to 7-rays. 

In this chapter, a model of the diffuse galactic emission will be presented. This 
model was developed at Goddard Space Flight Center by a team of scientists led by 
D.L.Bertsch. This model will then be used to extract a measurement of the IDGRB. 
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4.2 Cosmic-ray Interactions and 7- Ray production 

As discussed previously, cosmic-rays are relativistic charged particles. Most evidence 
points to a galactic source of the cosmic rays at least up to energies of a few TeV. 
Most theories speculate that supernovae are the predominant sources of these particles 
( [Biermann 1995| ). 

The interactions between such particles and galactic matter and photons are well 
understood in the EGRET energy range from terrestrial accelerator measurements. 
The dominant processes which yield 7-rays are: 

• hadronic showers yielding pions which decays to 7-rays. 

• cosmic ray electron bremsstrahlung. 

• inverse Compton interactions between relativistic electrons and soft photons. 

These three processes are characterized by the three 7-ray production functions 
qpion{E),qbrem{E),qjc{E) which measure the differential 7-ray energy spectrum pro- 
duced per target nucleon (or photon in the case of the inverse Compton function). 
These functions are determined using an assumed cosmic ray spectrum which is mea- 
sured locally. 



4.2.1 Nuclear Interactions with matter 

The dominant source of 7-rays above ~ 70 MeV is the decay of neutral and charged 
pions. Pions are produced in hadronic showers with center of mass energies greater 
than the rest mass of a pion. Neutral pions decays directly into 2 7-rays whereas 
charged pions decay into electrons and positrons which in turn annihilate to produce 
7-rays. Several calculations of the 7-ray production function from the decay of pions, 
Qpion, have been made ( [(Jaravallo 1971| , [Stecker 1970| , permer 1986| ). This function 



has a broad maximum at the rest energy of the pion, 68 MeV. 
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4.2.2 Electron Bremsstrahlung 

The electrons also contribute to the production of 7-rays through bremsstrahlung. 
The bremsstrahlung production function, qbrem, is difficult to ascertain because this 
interaction becomes significant at lower energies where the local electron spectrum 
is difficult to measure due to the solar modulation. Fichtel et al. ( [I'lchtel et al 



1991D have estimated the local electron spectrum using 7-ray data matched with 
synchrotron data. Their resultant bremsstrahlung production function becomes the 
dominant source of 7-rays at energies below ~ 20 MeV. 

4.2.3 Inverse Compton 

Cosmic ray electrons can also boost soft photons up to 7-ray energies through the 
inverse Compton process. 7-Rays produced through the inverse Compton process 
have a mean energy which is a function of the product of the electron energy and the 
photon energy, 

E-y = -7 -S^photon = -{ n) -^photon , (4.1) 

3 3 nieC^ 

where E-y, E^, -Ephoton are the energies of the resultant 7-ray, incident electron, and 
incident photon respectively. The two factors of 7 result from the Doppler shift of 
the target boost of the photon into the electron's rest frame and then the subsequent 
Doppler boost back into the observer's rest frame. In order to produce a 100 MeV 
7-ray using cosmic ray electrons from with energies from 1 GeV to 100 GeV requires 
target photons with wavelengths from ( 50 nm — 500 fim). As a result the important 
photon distributions include the microwave background. Far Infra- Red (FIR), Near 
Infra- Red (NIR), optical, and UV. 

4.3 Galactic Diffuse 7-Ray Model 



The EGRET team has developed a galactic diffuse 7-ray emission model ( [Bertsch] 



et al. 1993 ). This model was revised and analyzed spectrally by ( Hunter et al. 199B| ). 



The model will be outlined here. The galactic diffuse 7-ray intensity in a given 
direction (/, b), is given by the line of sight integral. 
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j{E, l,b) = — X j [Ce(p, /, b)qi,rem{E) + Cn{p, I, b)qpion{E)] 

x [riHiip, I, b) + nniiip, I, b) + nn^ip, i, b)]dp 

+ ^T.J('e{p,l,b)qic{E,p)uic{p,l,b)dp. (4.2) 

i 

In the above equation, Qpion, Qbrem, Qic represent the 7-ray production functions for 
nuclear, bremsstrahlung, and inverse Compton processes discussed above. These func- 
tions measure the 7-ray production per target atom when exposed to the local cosmic 
ray flux. The number densities of atomic hydrogen, ionized hydrogen, and molecular 
hydrogen are denoted by nHi,nHii,sind uh.^- These densities are the matter densities 
at the position (p, l,b). The cosmic rays densities arc also functions of position in the 
galaxy. The functions Ce{p,l,b),Cn{p,l,b) are the ratios of the electron and nucleon 
cosmic ray intensities relative to the local intensities. These functions could also 
depend on energy if the spectral indices vary with location within the galaxy. The 
second line of sight integral accounts for the inverse Compton contributions. The pro- 
duction functions arc calculated for four independent photon energy bands resulting 
in a summation over photon energy band. The functions Ui{p, I, b) reflect the photon 
densities in the i*^ band. 



4.3.1 Matter Distribution 

Equation 4.2 provides a recipe for calculating the 7-ray intensity given knowledge of 
the matter, cosmic ray, and photon distributions everywhere in the galaxy. In this 
section the state of knowledge of the matter distributions will be discussed. 

The bulk of the galactic mass is in the form of atomic hydrogen, HI. The distri- 
bution of this component is traced using the 21 cm hyperfine transition of hydrogen. 
Galactic rotation induces Doppler shifts which allow the three dimensional structure 
of this gas to be determined. Knowledge of the galactic rotation curve allows one to 
place the atomic hydrogen along a given line of sight at a given distance based on 
the projected velocity of that position along the line of sight. In the inner galaxy 



CHAPTER 4. DIFFUSE 7 RADIATION 



67 



(i.e. within the solar galactocentric radius) there is a near-far ambiguity that must 
be resolved. Matter at a given velocity can lie at either of two intersection points 
between the line of sight and the relevant line of constant galactocentric radius. In 
this case the ad hoc assumption is made to split the gas evenly between each position. 

In the inner galaxy there is a considerable of matter in the form of molecular 
hydrogen, This gas is not directly observable. The best available tracer for 

galactic molecular matter is the 2.6mm J = 1 —>■ rotational transition in CO 
molecules ( [Kutner fc Leung 1985|jMaloney fc Black 1988|) . Observations of this line 
can be deconvolved in the same way as the 21cm data in order to construct a model 
of the CO distribution in the galaxy. There remains the task of converting this 
measurement into a map of the H2 density. The conversion factor between CO line 
intensity and H2 density is given by the X-f actor, X. This value is a free parameter 
in the model presented here. 

The radio survey data used to construct the galactic matter distribution are de- 
scribed in ( [Bertsch et al. 1993| ). The deconvolved matter distributions that are con- 
structed using this data are shown in Figure which is reproduced from ( Hunter 



et al. 1996 ). The spiral arm structure of the galaxy is clearly visible. The dominance 
of the molecular hydrogen in the inner galaxy is also evident. 



4.3.2 7-Ray Production in the Galactic Plane 

In order to use this matter distribution to predict the galactic diffuse 7-ray emission, 
the distribution of cosmic rays throughout the galaxy must be known in order to 
perform the integral in equation 4.2. Because there is no direct measurement of cosmic 
ray density available except for the local measurement, some assumption about cosmic 
ray distribution must be adopted. The assumption used in generating the diffuse 7-ray 
model is that the cosmic rays are in dynamic balance with the magnetic fields and the 
gravitational attraction of the galaxy ( [Bertsch et al. 19*9^ ). The result of this balance 
is a coupling between cosmic rays and matter in the galaxy. The cosmic ray densities 
used to predict the diffuse 7-ray emission are thus derived through convolving the 
matter distribution with a characteristic smoothing function which is the second free 
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Figure 4.1: A reproduction of a figure from Hunter et al. (1996) showing tlie two 
dimensional plots of the matter distributions as seen looking face on our galaxy. The 
scale indicates the distance from the solar position which is marked with a cross. 
Figure (a) shows the distribution of HI, figure (b) shows the distribution of H2 and 
figure (c) shows the distribution of ionized hydrogen, HII. 
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parameter in the model. 

Figure [4.2| shows the resultant cosmic ray distribution which is derived fitting the 
two free parameters in the model to the EGRET data. The details of this analysis are 
described in ([Hunter et al. 19961) . The similarity between the map of cosmic rays and 
the input matter distribution is the direct result of the assumption about a coupling 
between these two quantities. 

An alternative approach to the determination of the cosmic ray distribution is 
to instead assume an azimuthally symmetric cosmic ray distribution in the galactic 
plane. The gradient of this distribution can then be determined from the 7-ray data. 
This approach was used by Strong et al. ( [Strong et al. 1988| ) in analyzing the COS B 
data. Figure shows a comparison between the radial cosmic ray distribution as 
determined under the assumption of dynamic balance from the ECRET data, and the 
cosmic ray gradient determined by Strong et al. using COS B data. In general the 
azimuthally symmetric method yields a smaller cosmic ray gradient than is predicted 
by dynamic balance. 



4.4 High Latitude Diffuse Emission 

The problem of modeling the diffuse intensity at high latitudes presents a slightly 
different set of challenges. [Sreekumar et al. 1996] have extended the galactic diffuse 



emission model to latitudes | & |> 10°. The density of neutral galactic gas (molecular 
and atomic) falls off rapidly away from the plane of the galaxy. As a result, the neutral 
gas along lines of sight with h > 10° is dominated by gas within a few hundred pc 
of the solar position. Local cosmic ray conditions are assumed to be valid at this 
range. This assumption obviates the need to deconvolve the matter distributions at 
high latitudes. The contribution from neutral gas is thus readily calculated using 
equation 4.2 where we use local electron and nucleon cosmic ray spectra and assume 
the value for the X parameter to be 1.5 consistent with the derived value from the 
galactic plane. 

The cosmic ray distributions undoubtedly extend to much larger distances from 
the galactic plane than does the neutral gas. As a result, the diffuse emission along 
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Figure 4.2: A reproduction of a figure from Hunter et al. 1996 showing two dimen- 
sional plots of the cosmic ray distributions as seen looking face on our galaxy. The 
scale indicates the distance from the solar position which is marked with a cross. 
Figure (a) shows the distribution of cosmic rays using the best fit parameters from 
the model. Figure (b) shows a comparison between the galactic radial distribution 
of cosmic rays using matter/cosmic ray coupling (solid line) and assuming azimuthal 
symmetry to the cosmic rays (dotted line). 
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high latitude hnes of sight contains significant contributions from 7-rays produced in 
interactions well above the galactic plane. There are three relevant interactions that 
need to be modeled: interactions between cosmic ray nucleons and ionized hydro- 
gen (HII), interactions between cosmic ray electrons and HII, and inverse Compton 
interactions between cosmic ray electrons and the soft galactic radiation field. 

In order to model the contributions from these large scale height components it is 
necessary to construct a model of the cosmic ray distribution above the galactic plane 
as well as the distribution of the various targets. The cosmic rays are assumed to be 
galactic in origin with most models predicting the source of all but the highest energy 
cosmic rays to be supernovae. The diffusion perpendicular to the plane of the galaxy 
is characterized by a cosmic ray scale height, Zp,Ze. A lower limit on the electron 
scale height can be extracted from galactic synchrotron radio data ( psborne 1995| ). 
This upper limit has been shown to be ~ 1 kpc. In an attempt to make a minimal 
set of assumptions in constructing our model, the cosmic ray electron scale height Ze 
is assumed to be at the lower limit of the possible range: Ze = 1 kpc. Furthermore, 
the cosmic ray nucleon scale height is assumed to be identical: Zp = Ze- Lastly 
the electron and nucleon cosmic ray spectra are assumed to be identical. A three 
dimensional model of the cosmic ray distribution is then constructed by projecting 
the cosmic ray densities, inferred from the galactic plane data, to higher latitudes 
using the scale height of 1 kpc. Uncertainties associated with these assumptions will 
be discussed in more detail below. 

The only available constraints on the distribution of ionized hydrogen come from 
dispersion measures of radio pulsars. Taylor & Cordes ( [Taylor fc Cordes 1993| ) have 
used a population of pulsars to constrain a model for the distribution of ionized matter 
in the galaxy which is adopted for this work. 

As discussed previously, the important the important photon distributions include 
the microwave background. Far Infra- Red (FIR), Near Infra- Red (NIR), optical, and 
UV. With the exception of the microwave background radiation, the photon distri- 
butions are spatially dependent and not very well constrained. The following input 
assumptions constitute the assumed photon radiation field, 

• Blackbody radiation: isotropic at a temperature of 2.7 K. 
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• FIR: the cold dust emission model of Cox, Krugel & Metzger (1996) with the 
total FIR luminosity of the galaxy normalized to 1.5 x IO^^Lq. The resultant 
distribution is radially symmetric with a peak at ~ 4 kpc and a scale height of 
100 pc 

• (NIR+Optical+UV): These components are stellar in origin and the model of 
Chi & Wolfendale (1991) is used. 

Given these models for the cosmic ray distribution, matter distribution, and pho- 
ton distribution, equation 4.2 is again used in order to generate diffuse 7-ray intensities 
along each line of sight. 

4.5 Model Predictions 

Figure [4.3| shows the three principal components to the galactic diffuse 7-ray emission. 
The first panel shows the contribution from neutral atomic and molecular hydrogen. 
Because we are only interested in the high latitude diffuse emission the color scale 
saturates on the galactic plane but reaches values as high as 8 x 10~^ in the galactic 
center. The second panel shows the contribution from HII and the third shows the 
inverse Compton component. All three panels show the integrated flux with E > 
100 MeV. 

4.6 Measurement of the IDGRB 
4.6.1 Likelihood Analysis 

The combined EGRET data set from the first four years of observations was used 
for the diffuse analysis described below. The photons used in this study were se- 
lected using the same criteria used in the instrumental calibration. This selection 
encompasses the standard selection criteria for E > 100 MeV except that the more 
restrictive albedo cut of a < 80° is used to avoid albedo contamination. Furthermore, 
only photons within 30° of the instrument pointing were included. A total of 630,416 
photons resulted from these cuts. 
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Figure 4.3: Tfie galactic contributions to tfie diffuse 7-ray emission. Tfie scales are 
in units of lO^^ph cm~^s~^sr~^ 



CHAPTER 4. DIFFUSE 7 RADIATION 



74 



These photons were binned in celestial coordinates using .5° x .5° bins. Celestial 
coordinates were chosen so as to avoid coordinate singularities at the galactic poles. 
The exposure maps were calculated for the appropriate data cuts. These maps were 
combined to generate a full sky exposure map. Two such maps were generated: one 
using the preflight calibration and a second using the in-flight correction outlined 
in the previous chapter. The diffuse model outlined above was convolved with the 
EGRET point spread function to generate an observational background model. 

These three maps were analyzed using the EGRET mmaximum likelihood soft- 
ware (LIKE V 4.5) ( pVlattox et al. 1996| ) (see also Appendix A). This package searches 
for point sources using the maximum likelihood method. To this end, a 3 free pa- 
rameter fit was made to the data at every point on the sky. The 3 parameters that 
were estimated using this method were the source flux (5) of a putative source at a 
given location (Zo,&o) as well as 2 diffuse parameters (/iDGRB,5'm) which measure the 
IDGRB intensity and the galactic diffuse emission model scaling, respectively. 

/(/, b) = Sx PSF{1, b; lo, bo) + gmG{l, b) + Jidgrb • (4.3) 

These 3 parameters were optimized at a given point using data within 15° of that 
point. 

Following the procedures used to generate the 2nd EGRET source catalog ( P?homp- 
son et al. 1995|) , this process was carried out until all sources with significances a > 4 
were found. The number of sources isolated using this method was 147. The point 
source detections will be discussed in detail next chapter. Maps of the 2 diffuse 
parameters were also generated. 

If the 2 diffuse sources were totally independent, the maps of the diffuse parameters 
would provide measurements of the IDGRB intensity. However at high latitudes, the 
relative featurelessness of the galactic diffuse emission causes the diffuse parameters 
to be highly correlated. As a result, it is only possible to extract the sum of the two 
components. 
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4.6.2 Residual Emission 

In order to find the level of isotropic diffuse emission that is consistent with the 
EGRET data, the galactic diffuse model is directly correlated with the measured 
intensity. To properly account for the effect of the point sources, the point sources 
isolated using the likelihood technique are added to the galactic diffuse model using 
the appropriate point spread function. Figure shows the measured intensity as well 



as the modeled intensity including the point sources reconstructed using maximum 
likelihood. The intensities in each case have been smoothed using a Gaussian with a 
5° width in order to decrease the Poisson noise and the galactic plane has again been 
saturated in order to reveal the structure at high latitude. 



The first panel in Figure |4.6| shows the difference between these two quantities. 
This is a map of all the intensity in the sky that is not associated with resolved point 
sources and is not accounted for in the galactic diffuse model. There is clear evidence 
of excess diffuse emission at high latitudes. 

In order to ascertain the level of the IDGRB a set of 6° x 6° in celestial coordinates 
are used. The measured intensity is correlated to the galactic diffuse model intensity 
in these pixels: 

lij = QmGij + /iDGRB • 

(4.4) 

The parameters Qm, /idgrb are estimated by linear regression. A one parameter fit is 
also performed using a fixed value of = 1. 

While in theory, one could use data from the entire sky to contribute to this 
analysis, this approach is complicated by several factors. Firstly, the galactic plane 
emission is very intense compared with the corresponding high latitude intensities. 
Small inaccuracies in the galactic model at low latitudes could severely bias the anal- 
ysis and as a result regions within 30° of the galactic plane have been excluded. 
Secondly, regions in the direct vicinity of strong point sources are dominated by the 
emission from such sources and the residual emission is likely to be dominated by 
small inaccuracies in the fitting of the point source. As a result, regions within 3° of 
any of the 147 resolved point sources are excluded. 

It is evident from inspection of the residual intensity in Figure ffTBI that there is a 
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Figure 4.4: Measured intensity and modeled galactic diffuse intensity. The likelihood 
reconstructed point sources with significances a > 4 are also included. The scales are 
in units of 10~^ph cm~^s~^sr~^ 
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Figure 4.5: Measured diffuse intensity versus modeled galactic diffuse intensity for a 
set of high latitude pixels. Regions within 90° of the galactic center are excluded. 
The extrapolation to zero galactic diffuse intensity yields a measurement of /idgrb = 
1.33 ± .03 X 10"^. The departure of x^/'^ from 1.0 indicates that local diffuse model 
inaccuracies remain. 



prominent feature surrounding the galactic center which is not the result of isotropic 
emission. In order to gauge the effect of this feature on the measurement of the 
IDGRB, regions within ^gc of the galactic center are excluded from the study, ^^gc 
is allowed to vary from 40° to 110° in order to determine whether the effect of this 
feature extends to large viewing angles. 

Table 4.1 summarizes the results of this study. For values of ^^gc < 90° the 
values of Qm depart significantly from 1, indicating that the galactic diffuse model is 
incorrectly predicting the 7-ray intensity. However, the fact that the values of x^/^ 
depart significantly from 1 in these regions indicates that the data are inconsistent 
with the additional emission being isotropic. For values of ^qc > 90° the influence of 
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Table 4.1: Residual intensity as a function of sky region. 





1 Free Parameter 


2 Free Parameters 




-^IDGRB 
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xlO" 
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xlO" 
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1.23 ± 
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1.30 ± 


.04 


.98± .05 


1.6 
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1.26 ± 


.01 


1.4 


1.34 ± 


.03 


.95± .05 


1.4 


90 


1.28 ± 


.01 


1.6 


1.33 ± 


.03 


.99 ± .04 


1.6 


80 


1.29 ± 


.01 


1.9 


1.32 ± 


.03 


1.02 ± .04 


1.8 


70 


1.32 ± 


.01 


2.3 


1.30 ± 


.03 


1.07 ± .04 


2.3 


60 


1.32 ± 


.01 


2.5 


1.30 ± 


.02 


1.10 ± .03 


2.4 


50 


1.35 ± 


.01 


2.6 


1.27± 


.02 


1.18 ± .03 


2.5 


40 


1.37 ± 


.01 


2.7 


1.24 ± 


.02 


1.24 ± .03 


2.5 



this additional emission is not important. In these regions, the gas- map appears to 
correlate well with the measured intensities as evidenced by values of Qm very close 
to 1. Furthermore, values of x^/^ close to 1 indicate that the additional emission is 



consistent with the hypothesis that this emission is isotropic in nature. Figure ^"5 
shows the explicit correlation for Oqc = 90°. The measured value of /idgrb is (1.33 ± 
0.03) X 10^^ ph cm-2s-isr~\ 



4.6.3 Likelihood Fit to the Diffuse Background 

In order to extract the galactic bulge feature, the likelihood fit to the diffuse param- 
eters is used. While these parameters cannot be directly separated due to their high 
degree of correlation, they do give a measure of the best fit to the diffuse intensity 
within 15° of each point on the sky. These values can then be used to construct a 
maximum likelihood diffuse map in the following way. 



-^ij 9inijCrij ~\~ I] 



IDGRBij 



(4.5) 



This method gives a model of the total diffuse emission which has been fit to the data 
on large scales. 

In order to analyze the large scale features, the galactic diffuse model is subtracted 
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Residual Intensity using the high latitude galactic diffuse model 




Galatiic LongFtude 

Small Scale Residual Intensity 

Figure 4.6: The residual intensity after subtracting a galactic diffuse model and re- 
solved point sources. The scales are in units of 10~^ph cm~^s~^sr"~^ 
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from the maximum likelihood diffuse model as is the IDGRB as measured above. 
What remains is the residual emission on large scales. A map of this emission is 



shown in the second panel of Figure ^ 



4.6.4 Possible Sources of Large Scale Residual Diffuse Fea- 
tures 

The most prominent diffuse feature that remains after subtracting our model of the 
7-ray intensity is a prominent feature around the galactic center. The position of this 
feature in the galactic bulge is strongly suggestive that it is the result of imperfectly 
modeled inverse Compton emission. A few of the uncertainties in this calculation will 
now be discussed and their impact on the data assessed. 

There is much uncertainty associated with the photon distributions used in this 
model. The FIR photon distribution was taken to be 100 pc ( Sodroski et al. 1994| ). 



It is possible that these values are biased by the fact that photons that escape the 
galactic disk are not observed directly despite the fact that they can serve as target 
photons for cosmic ray electrons above the galactic plane. A scale height of 1 kpc 
for these photons would increase the inverse Compton contribution by ~ 20% at 
mid-latitudes. 

Furthermore, the absolute normalization of the FIR photon distribution is un- 
certain. Sodroski et al. ([Sodroski et al. 1994|) quote a galactic FIR luminosity of 



1.6 - 2.2 X 1O^°L0 using COBE/DIRBE data whereas Wright et al. ([Wright et al. 



T99TD report a value of (1.8 ± 0.6) x IO^^Lq from an analysis of COBE/FIRAS. It 



is clear that the total FIR luminosity could be increased by 20% without being in 
conflict with the available FIR data. 

Similar uncertainties exist in the work of Phi fc Wolfendale 1991] . Both the photon 



scale heights and normalizations are uncertain to within 30%. 

Finally, the cosmic ray scale height itself is very uncertain. The value adopted for 
this analysis was the lower limit from the deconvolution of the galactic synchrotron 
data. Youssefi & Strong ( [Youssefi fc Strong 1995[) used the opposite approach in 
analyzing the COS B diffuse 7-ray data. In this work, the authors use the 7-ray 
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data to constrain the inverse Compton emission leaving the cosmic ray scale height 
as a free parameter. The reported electron cosmic ray scale height is found to be 
3 — 5 kpc. The same authors have proposed a cosmic ray diffusion model ([Strong 



fc Youssefi 1995|) . Under assumptions about the galactic magnetic field distribution 
as well as the cosmic ray injection spectrum and distribution, an energy dependent 
cosmic ray electron scale height is derived. This scale height peaks at ~ 1 GeV after 
which inverse Compton losses tend to contain the cosmic rays to a tighter distribution 
about the galactic plane. While more data is needed in order to determine the electron 
scale height accurately, it is not unlikely that much of the residual surrounding the 
galactic center is the result of additional a inverse Compton emission. 

More puzzling, perhaps, than the bright residual in the galactic bulge is the North- 
South anisotropy at high latitudes. In the North Galactic polar region there are 
several large scale residual excesses. The region extending from the Virgo region to the 
North Galactic pole appears to be significantly brighter than modelled. Furthermore, 
the South Galactic pole intensity is significantly underpredicted by the model. A 
third region near (/, b) = (70, 60) has been reported as a region of excessive diffuse 
emission (|Chen et al. 1995|) . However, the proximity of this region to the orbital pole 
in the principal observation of this region means that this observation was most likely 
contaminated by albedo photons. Reanalysis of this region using the more restrictive 
zenith angle cut also reveals a point source in this region (see next chapter). The 
remaining excess in this region is most likely a statistical fluctuation. 

The asymmetries at such high latitudes argues in favor of a local source of these 
anisotropies. A correlation between the northern excesses and the large radio super- 
nova remnant features (Loop 1 and Loop 3) has been suggested ( Ps borne 1995|) but 
the evidence is not very compelling. The most likely source of these features is a local 
variation in the cosmic ray density within a few hundred parsecs. 

The underprediction of the model near the South Galactic pole is potentially due 
in part to the sensitivity drift of the EGRET instrument which is caused by gas 
degradation. This effect is accounted for by looking at overlapping diffuse regions 
and fitting the sensitivity trend versus time. Several of the dominant observations 
of the South Galactic Pole exhibit apparent sensitivities which are below the trend 
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indicating that there may have been slower than average degradation during these 
observations. 



4.6.5 Effect on the Measurement of the IDGRB 

In hght of the problems with the galactic diffuse model the impact on the measurement 
of the IDGRB must be assessed. In the limit of infinite cosmic ray scale heights, 
the total isotropically distributed 7-ray emission that would result from interactions 
between cosmic rays and microwave background photons well above the galactic plane 
amounts to ~ 20% of the IDGRB (|Sreekumar et al. 1996|) . Thus it is clearly impossible 
to attempt to explain the IDGRB as the result of an extended cosmic ray halo. 
However, it is clear that the systematic uncertainties in the foreground model limit 
our ability to measure the IDGRB intensity. These uncertainties are significantly 
larger than the statistical uncertainties. 



4.7 Spectral Measurement of the IDGRB 

While this study is primarily concerned with extracting the spatial signals from the 
diffuse background, the spectrum of the IDGRB is also of critical importance. A 
measurement of the IDGRB spectrum will be outlined here. This work follows the 
basic techniques outlined in ([Kniffen et al. 1995|). A more detailed analysis of the 



IDGRB spectrum is to appear in ( [Sreekumar et al. 1996| ). 



Maps of the summed EGRET counts through Viewing period 404.0 were gen- 
erated in ten energy ranges: {E (MeV): 30-50,50-70,70-100,100-150,150-300,300- 
500,500-1000,1000-2000,2000-4000,4000-10000). Corresponding exposure maps were 
also generated. To minimize the albedo contamination in making this sensitive mea- 
surement, an albedo cut of a < 80° was again adopted for all energy ranges. The 
sky at latitudes | & |> 30° was divided into equal solid angle (20° x 20°) bins in 
celestial coordinates. Regions within 5° of a strong point source were excluded from 
the analysis. In order to avoid the regions in which the unmodelled inverse Compton 
emission is important, regions within 90° of the galactic center are excluded from the 



CHAPTER 4. DIFFUSE 7 RADIATION 



83 



analysis as well. 

The galactic diffuse contribution to the intensity of each pixel was calculated 
using the extrapolation technique described above. The scaling of the galactic diffuse 
model was left as a free parameter in each energy range which was determined by 
the correlation between measured intensity in each pixel and the modelled intensity 
in that pixel. In each energy range this scaling was consistent with unity except in 
the highest three ranges in which the galactic model significantly underpredicted the 
diffuse emission. For a discussion of this effect see Hunter et al. ([Hunter et al. 1996| ). 

The number of residual non-galactic counts in each bin in each energy range 
were thus determined. The spectrum was then deconvolved from this data using 
the SPECTRAL program ( [Nolan et al. 199^ ). This program performs a maximum 
likelihood fit to the measured counts amongst a space of power law input spectra 
which are folded through the EGRET energy response functions. 



Figure shows several of the measured spectra across the sky. The map indicates 
the regions for which the spectra are shown. Even with the stringent data cuts that 
are made in order to obtain clean a clean data set, well constrained spectra are 
available. All the spectra are good fits to a power law. There is a slight indication 
that there is an excess in the 2000-4000 bin and a deficit in the 4000-10000 bin in all 
the spectra. This is due to an inaccuracy in the energy determination of the Class B 
events in which an insufficient amount of energy was deposited in the TASC to pass 
the photo-multiplier threshold. These events had to have their energies calculated 
from their track patterns. The power law indices are all within one or two standard 
deviations of each other. 



The distribution of the spectral index across the sky is shown in Figure O. The 
variation in the error bar size indicate the variable statistical errors in each mea- 
surement due to exposure variation as well as the truncation of some of the analysis 
regions by their proximity to point sources. The weighted mean is also given in this 
figure. The best fit value is a = 2. 22 ±.01. Also shown is the minimum x^/^ value for 
the data set. The elevated value of this number is indicative that systematic errors 
in the galactic diffuse model dominate over statistical errors. There is no evidence 
for a smooth variation of the spectral index across the sky. 
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Figure 4.7: Four independent measurements of the IDGRB spectrum. The regions 
used to calculate each spectrum are shown in the top panel. Each region is 20° x 20°. 
Regions within 5° of strong point sources are excluded. There is good consistency 
between the best fit spectra. 
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Figure 4.8: Several independent measurements of the IDGRB spectral index (a). The 
error bars represent the statistical errors from the deconvolution. The high value of 
X^/z/ indicates that the systematic errors in the galactic diffuse model are dominant. 



The cosmic background spectrum can now be placed in its spectral context. Fig- 



ure |4.9| shows the best available measurements of the cosmic background radiation 
from hard X-rays through the EGRET range. The X-ray data is from HEAO-1 as 
analyzed by Kinzer et al. ( Kinzer et al. 1996| ). The data in the MeV region is from 
COMPTEL ( [Kappadath et al. 1995|) . The dramatic difference between this and other 
measurements of the background in this spectral region results from the fact that 
instrumental background was able to be subtracted through its rigidity modulation. 
Previous balloon measurements were performed at a single rigidity, prohibiting the 
use of this powerful technique. 

Also shown on in Figure [4.9| is the contribution of X-ray Seyferts to the back- 
ground ( [Zdziarski et al. 1995|) . Essentially all of the hard X-ray background is now 
believed to be due to unresolved X-ray Seyferts. However, these spectra are well fit 
by thermal comptonization spectra which are exponentially suppressed in the 7-ray 
region. Indeed no X-ray Seyferts have been detected by EGRET . These sources do 
not contribute substantially to the background above 1 MeV. 
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The remaining high energy spectrum is a fairly consistent power law throughout 
the measured spectral region. There is no evidence of spectral features that would 
indicate the presence of a high energy particle decay. Furthermore, the absence of 
the 'MeV bump' eliminates the strongest evidence in favor of an annihilation feature 
in that spectral region. 

4.8 Discussion 

EGRET has proved to be a powerful tool in the study of high energy diffuse 7-ray 
emission. The large effective area, good spatial resolution, and low instrumental 
background have allowed measurement of both galactic and extra-galactic diffuse 
emission with unprecedented accuracy. 

A model of the diffuse 7-ray emission has been presented. This model with only 
two free parameters is an excellent fit to the observed 7-ray data. At intermediate 
galactic latitudes toward the galactic center there remains a clear deficiency in the 
model. The high residuals in this region are almost certainly the result of unaccounted 
for inverse Compton interaction between galactic cosmic ray electrons and soft back- 
ground photons. The poor constraints on the distribution of cosmic ray electrons 
remains a difficulty in the construction of more accurate models. The 7-ray data can 
be used to constrain this distribution but such work is difficult because of the many 
uncertainties involved. 

Despite these difficulties, there is clear evidence of an IDGRB that cannot be 
explained by the known galactic sources of diffuse emission. It has been shown that 
the intensity of this emission above 100 MeV is, 

/iDGRB = (1.33 ± .03 X 10"^ photons cm^^ s"^ sr"^ . (4.6) 

The statistical error quoted is small compared to the systematic uncertainties involved 
in the calculation. Most importantly the poorly constrained cosmic ray scale height 
allows for up to 20% of this emission being galactic in nature. 
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Spectral measurements of the IDGRB have shown that this emission is well de- 
scribed by a power law with spectral index of a = 2.22 ± .01. This spectrum does not 
appear to vary across the sky lending credibility to the hypothesized extra-galactic 
origin. The spectrum is harder than measured by SAS 2 but the measured value is 
within the uncertainty of that measurement. This spectrum matches well to the new 
Comptel spectrum. The combination of the harder high energy spectrum and the 
disappearance of the 'MeV bump have most likely ended speculation about Baryon 
Symmetric Big Bang theories of the origin of this radiation. The fact that this spec- 
trum is now consistent with measured AGN spectra has fueled speculation that this 
emission is due to unresolved members of that population. This notion will be ex- 
plored in detail in the subsequent chapters. 
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Figure 4.9: The cosmic background radiation spectrum from hard X-rays through 
the EGRET range. The X-ray data is from recently re-analyzed work by Kinzer et 
al. 1996 in which new measurements of the instrumental background were obtained. 
The McV data are from the COMPTEL instrument. The pronounced spectral bump 
at ~ 1 MeV is not evident in this data after a careful measurement of the rigidity 
dependence of the background (Kappadath et al. 1995). Also shown is the spectrum 
of calculated Seyfert contribution to the background (Zdziarski et al. 1995). 



Chapter 5 

High Latitude Point Sources 



Prior to the launch of EGRET, the only extra-galactic source detected in the high en- 
ergy 7-ray band was the quasar 3C 273 ( |Bignami et al. 1981| ). One of EGRET^s most 
dramatic results has been the detection of dozens of sources at latitudes greater than 
10°. The majority of these sources have been identified as being extra-galactic based 
on their positional co-incidence with bright Flat Spectrum Radio Quasars (FSRQ's) 
( [von Montigny et al. 1995 ). This discovery has focussed the discussion about the 
origin of the IDGRB on the contribution from active galactic nuclei (AGN). In this 
chapter, the nature of the high latitude point sources detected by EGRET will be 
examined and several models to explain their distribution will be discussed. 



5.1 High Latitude Point Source Analysis 

The second EGRET catalog ( P?hompson et al. 1995|) summarizes the point source 
analysis of the EGRET data through August 17 1993. In this study the maximum 
likelihood technique is used to locate point sources and to determine their fluxes. This 
technique is applied to each viewing period individually as well as to the summed all- 
sky data. The point source analysis technique is discussed is detail in ( [Mattox et al. 
|1996|) and is summarized in Appendix A. 

The point source and diffuse model parameters are fit to the data within a 15° 
radius of a given test point. In this way local inaccuracies in the diffuse model such as 
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the large scale underprediction of the model in the galactic bulge region, do not bias 
the analysis. The point source significance is then estimated using the distribution of 
the likelihood statistic in the null hypothesis (see Appendix A). This analysis allows 
the simultaneous determination of the value of the diffuse backgrounds at a given 
point as well as the significance and flux of a putative point source at this location. 

In order to simultaneously analyze the many point sources that EGRET has 
discovered, it is necessary to account for the coupling between sources that is the 
result of a flnite point spread function. This is done using an iterative technique. A 
likelihood map is constructed and sources above a given significance threshold are 
extracted. The three point source parameters {I, b, S), for each source are then itera- 
tively adjusted and the effect on the likelihood is determined. Because simultaneous 
optimization of all three parameters for each source would be computationally infea- 
sible, the source fluxes are adjusted flrst and then held constant while the positions of 
these sources are adjusted. When a simultaneous maximum has been achieved, these 
sources are then added to the diffuse map as background features and the data are 
reanalyzed to see if any additional sources that were previously obscured are thus re- 
vealed. All sources arc then simultaneously optimized. This process is then repeated 
until no sources above a given significance threshold remain. 

The above process is reliable so long as the detected point source density is suf- 
ficiently low so as not to create false convergences under optimization. Simulations 
have shown that sources with significances of greater than or equal to Aa can be reli- 
ably reconstructed using this technique (sec Appendix A). It is, however, impossible 
to push this technique to lower significance as the results tend to change with the 
order in which the sources are added. 

The point source analysis of the summed data carried out for the Phase I-I-II 
catalog has been redone with the following modifications: 

• Data through 11 April 1995 (Viewing Period 419.1) are included. 

• Only regions with |6| > 20° are analyzed. 

• A zenith angle cut of 80° is adopted. 
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5.2 Point Source Results 

Table 5.2 summarizes the results of the high latitude point source analysis. A total 
of 55 sources were detected with \b\ > 20°. Of these 33 are blazars which were also 
detected in Phase I+II (P12) data. Another 4 can be identified with blazars that 
were undetected in the P12 data. The only other identified high latitude source is the 
Large Magellanic cloud (LMC) which is also seen in the P12 data (jSreekumar et al" 



19921) . The remaining 17 sources are unidentified. Of these 6 were detected in the 
PI 2 data. The remaining 11 sources are new unidentified sources. The fact that most 
of these sources are near the significance threshold is indicative of the fact that they 
are sources rendered visible through increased exposure. At least one of these sources 
( / = 73.09,6 = 68.18) is visible in the P12 data if one adopts the more restrictive 
albedo cut made in this analysis. This is a result of the fact that this source appeared 
near the orbital pole during its observation which increases the albedo background as 
discussed in Chapter 3. This source contributes to the diffuse emission enhancement 
reported by (|Chen et al. 1995| ). 



5.2.1 Time Variability 

Before it is possible to analyze the source flux distribution, it is necessary to under- 
stand the time variability of high latitude sources. The time variability of blazars is 
most dramatically exemplified by the 7-ray flare that occurred in June 1991 in the 
source 3C279 ( [Knitfen et al. 1995| ). In this observation the 7-ray flux was seen to 
change on the time scale of several days. This short time scale variation in such a 
luminous source is compelling evidence that the source is compact and that a jet is 
involved in the acceleration mechanism. 

The data from the Phase I+II catalog are used here to address the question of time 
variability for the blazars as a class. Figure ^]l] shows the distribution of fluxes in each 
independent observation of a source normalized by the mean flux calculated using the 
summed data. Only observations in which the source was detected with cr > 4 are 
shown. This results in an asymmetry about unity. The error bars are calculated by 
compounding the uncertainty of a given measurement with the uncertainty of the 
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Table 5.1: Results of the likelihood analysis performed using all data up to VP 419.1 
using an 80° zenith angle cut. 



Category 


I 


b 


S(E > lOOMeV) 


a 


P12? 


ID 


Blazars 


23.52 


41.02 


27.0 ± 5.1 


6.7 


X 


1606-1-106 


Detected in P12 


35.98 


-24.68 


22.8 ± 4.2 


6.7 


X 


2022-077 


(=33 sources) 


55.81 


46.02 


53.1 ± 7.7 


9.9 


X 


1611-1-343 




61.69 


42.16 


40.4 ± 6.7 


8.0 


X 


1633+382 




77.71 


-39.09 


20.9 ± 3.0 


8.8 


X 


CTA 102 




79.43 


31.60 


28.4 ± 5.1 


7.1 


X 


1739-H522 




85.98 


-37.97 


68.5 ±4.4 


22.2 


X 


30 454.3 




143.54 


34.55 


12.1 ± 2.1 


7.0 


X 


0836-1-710 




143.55 


28.12 


19.0 ± 2.3 


10.5 


X 


0716-1-714 




147.65 


-44.57 


10.0 ± 3.3 


3.6 


X 


0202+149 




150.17 


-28.73 


12.1 ± 2.9 


5.0 


X 


0234+285? 




156.29 


-38.83 


32.5 ± 4.4 


9.4 


X 


0235+164 




159.67 


47.54 


9.0 ± 1.8 


6.1 


X 


0954+556 




169.96 


32.09 


18.0 ± 3.0 


7.5 


X 


0804+499 




177.00 


44.28 


13.8 ± 2.2 


7.9 


X 


GRO J0916+43 




179.89 


64.89 


14.0 ± 2.0 


9.0 


X 


Mrk 421 




-172.88 


-20.34 


20.2 ±3.1 


7.5 


X 


0446+112 




-165.21 


-32.17 


19.8 ± 3.7 


6.3 


X 


0420-014 




-160.11 


31.36 


23.7 ± 4.1 


7.6 


X 


0827+243 




-168.94 


-26.36 


10.4 ± 2.8 


4.2 


X 


0468-020 




-118.99 


-31.54 


19.3 ± 3.8 


6.3 


X 


0521-365 




-109.97 


-30.84 


24.4 ± 3.5 


9.0 


X 


0537-441 




-83.78 


-61.78 


71.2 ± 6.7 


16.1 


X 


0208-512 




-83.05 


43.53 


10.2 ± 2.5 


4.7 


X 


1127-145 




-69.30 


64.07 


19.5 ± 1.9 


12.3 


X 


3C 273 




-55.06 


56.98 


81.2 ± 3.0 


39.8 


X 


3C 279 




-51.52 


28.09 


17.2 ± 2.9 


7.0 


X 


1313-333 




-25.73 


50.48 


26.6 ± 2.9 


11.9 


X 


1406-076 




-9.00 


40.39 


15.8 ± 4.0 


4.6 


X 


1510-089 




-105.92 


81.70 


15.2 ± 2.0 


9.4 


X 


1222+216 




168.19 


-77.32 


13.6 ± 3.7 


4.6 


X 


0130-171 




145.92 


44.10 


9.2 ± 1.8 


5.8 




0954+658 




107.12 


-41.72 


11.7 ±3.3 


4.2 


X 


2356+196 


New Blazars 


-162.64 


-28.84 


20.1 ± 3.8 


6.4 




2251+158 


(=4 sources) 


-161.24 


83.43 


11.2 ± 2.0 


6.7 




1219+285 




27.03 


20.65 


17.2 ± 4.3 


4.6 




1725+044 




17.76 


-52.26 


12.4 ± 3.4 


4.5 




2155-304 


Other Identiiied 


-80.50 


-32.28 


14.4 ± 2.8 


6.1 


X 


LMC 


Unidentified 


88.63 


25.05 


69.8 ± 5.4 


18.2 


X 


GRO J1837+59 


(=6 sources) 


163.14 


29.11 


13.9 ± 2.8 


6.1 


X 


GRO ,10744+54 




179.64 


-22.86 


22.7 ± 3.3 


8.0 


X 


GRO .10421 + 15 




-66.22 


67.09 


10.7 ± 1.7 


7.. 3 


X 


c;rc) ji2:s7+04 




-44.!)!) 


-24..S1 


15. li ± 5.1 






( :R< ) .Ti7:il-78 




147.40 


5:i.:il 


l).l ± l.(i 


4.(i 


X 


GRCJ ,11047-58 


New Unidentified 


15.93 


24.59 


15.1 ± 3.8 


4.5 






(=11 sources) 


73.09 


68.18 


15.6 ± 4.8 


4.2 








-176.46 


-31.74 


12.8 ± 3.3 


4.4 








-119.60 


-21.53 


13.2 ± 3.4 


4.7 








-63.99 


58.23 


7.5 ± 1.7 


5.0 








10.71 


23.62 


14.5 ± 3.5 


4.6 








-83.78 


-56.59 


12.6 ± 3.6 


4.2 








-152.87 


35.39 


11.3 ± 3.2 


4.3 








-48.29 


57.30 


7.2 ± 1.8 


4.4 








4.25 


-24.03 


11.0 ± 2.9 


4.2 








-117.68 


-47.12 


9.5 ± 2.7 


4.3 







mean. The distributions are shown for blazars as well as for unidentified sources. 

The blazars show clear evidence of time variation on the time scale of weeks (the 
average duration of a viewing period). This is evidenced by high value of jv- The 
unidentified sources show weaker evidence of variation. This may be evidence that 
some fraction of these sources are not blazars but some foreground population of time 
independent sources. However, this conclusion is not straightforward because of the 
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upper limits are not included in this method. The distributions shown are truncated 
at the detection threshold. The effect of this is to make the distributions of weak 
sources cluster about the detection threshold. Because the unidentified sources are 
weaker on average than the blazars, the effect is to introduce a bias which tends to 
remove the evidence of time variability for these sources. 

The observations of blazars in which upper limits were obtained can be collec- 
tively analyzed to ascertain whether there is evidence of continued quiescent flux 
during these periods. This is done by examining the distribution of the likelihood 
test statistic at the positions of known sources during periods in which only an upper 
limit was obtained for their flux. This distribution is then compared to the a similar 
distribution of test points which were randomly selected in the sky avoiding regions 
near known point sources. The result is shown in Figure p]^. There is clear evidence 
for the continued emission during these 'off' periods of the blazars. 

There are several ways to calculate the source flux distribution in light of this time 
variability of the sources. The simplest method which is described above is to sum 
all the data together and calculate the source flux. This can be called an exposure 
weighted mean in the sense that observations near the the instrumental axis and 
during times of high instrumental efficiency (when the effective area, Ai, is high) are 
more highly weighted in the mean flux as they provide a greater fraction of the total 
counts C, 

-Q X/j C-i X/j ^i^i X/j SiAjt^ -\\ 

•J Summed c c T"I ' l'-'--'^/ 

The result is that some sources are dominated by one or two observations and the 
flux is not averaged over a very long period. The more meaningful average flux is the 
time averaged flux, 

-S, = . (5.2) 

The difficulty in calculating this sort of average is the problem of how to treat upper 
limits. For the purpose of this study only upper limits that are inconsistent with the 
exposure averaged flux are considered. These observations are assigned a flux equal 
to half the upper limit flux. 

Figure |5.3| shows the integral source distributions {logN — logS distributions) for 
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Blazar Observation 



Blazar flux variations. 



10 20 
Unidentified Source Observation 



Unidentified source fiux variations. 

Figure 5.1: Distributions of flux in a given viewing period normalized by the mean 
flux of that source. Only detections of greater than Aa are shown. The asymmetry 
about 1 is caused by the absence of the upper limits for which a flux measurement is 
impossible. The high value of x^/^ for the blazars demonstrates strong time variabil- 
ity in these sources. There is weaker evidence for this variation in the unidentified 
sources although this may be result of truncation at the detection threshold for these 
comparatively weaker sources. 
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Figure 5.2: The distribution of tiie likelihood test statistic (TS) for a set of obser- 
vations of blazars for which only upper limits were obtained. This is compared to a 
set of control points away from any strong sources. The result clearly indicates that 
these sources do not completely turn off but rather exhibit some quiescent flux. 

the high latitude blazars calculated using different flux measures. The summed flux 
distribution is shown as a gray line. The time averaged flux distribution is shown as 
a black line. Also shown for reference is the distribution of peak fluxes. The two flux 
averages are similar but there is some evidence of a steepening of the distribution as 
one averages over longer times. 

5.2.2 Spatial Distribution 

If the sources are extra-galactic in origin, their spatial distribution should be uniform. 
In order to verify this, it is necessary to be able to account for the variations in 
detection threshold across the sky. The detection threshold is determined jointly 
by the exposure and background at a given point. To a lesser extent, the spatial 
distribution of the background also affects a source's ability to be detected (e.g. it is 
hard to detect a point source on top of a spike in the background than a comparably 
intense flat background). A simple model of the relationship between flux, S, and 
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Figure 5.3: Integral source distributions for tlie high latitude blazars. Three different 
measures of source flux are shown. The distribution to the right uses the peak flux 
measured during the mission. The other two distributions represent two measures of 
mean flux. The gray line shows the distribution using the exposure weighted (summed 
data) flux, while the black line shows the distribution using the time averaged flux. 
There is a evidence that the distribution steepens as one averages over longer times. 

significance, a, of a source is given by the statistical fluctuation of the background. 



where Cs,Cb are the counts from the source and background respectively. For low 
signiflcance sources, the signiflcance is dominated by the background fluctuations 
Cb > Cs- In this assumption, the signiflcance of a source of strength S to which 
there is exposure £ and which sits on a diffuse background of strength B is. 



Cs 




a — 






where /c is a constant which is determined by the point spread function. In this way 
source signiflcance is shown to be proportional to its flux near the detection threshold. 
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This relationship was fit to the data from Phase I+II. The parameter k was evalu- 
ated for all sources with 4.0 < cr < 10.0. The exposure was taken from the calculated 
exposure map and the background was assumed to be represented by the best fit to 
the diffuse background (galactic and extra-galactic). The resulting distribution of 
k is shown in Figure |5.4| . The distribution is sharply peaked about k = 0.23. The 
several orders of magnitude variation in the different constituent parameters has been 
successfully scaled out. 



S 40 ^ 
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Figure 5.4: The distribution of the quantity k defined in equation 5.4 for sources 
detected in Phase I+II. a,S,B,S are respectively the source significance,flux,diffuse 
background at the source location, and exposure. The tight distribution of this param- 
eter indicates that the several orders of magnitude variation in the input parameters 
have been successfully scaled out. 



This model allows us to determine a 4a detection threshold map for the data set 



in question. Figure q]5| the detection threshold map for the summed data through 
VP 419.1. Note that the non-uniform exposure creates almost an order of magnitude 
variance in the detection threshold at high latitude. 

The source distribution can now be analyzed to check for uniformity. For a given 
range of |6| the predicted number of sources can be calculated using some expected 
source flux distribution. A Euclidean source distribution with A^(> 5*) = AS^^'^ 
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Figure 5.5: The Aa detection threshold map for the combined data through VP 
419.1 calculated using the above method. The values at high latitudes range between 
5 X 10"^ and 2 x 10"^. 



is used. The source distributions shown above appear to be consistent with this 
assumption. The total number of sources, Nf,, expected in a band of latitude 6\b\ is 
thus given by, 

m = [ As^tl^,,^dn . (5.5) 



S\b\ 



The ratio of this number to the number of sources detected in the same latitude range 
is a measure of the true source density on the sky. 



Figure ^]6| shows the corrected source density as a function of latitude for the 
blazars as well as for the unidentified sources. The distribution of blazars appears to 
be uniform as expected. The unidentified sources show a possible excess at medium 
latitudes between 20° and 30°. Kanbach ( [Kanbach et al. 1996| ) has used similar 
analyses to suggest that on the basis of the spatial distribution as well as the time 
variability differences between the two populations of high latitude sources, there is 
evidence for a local galactic component in the unidentified sources. The evidence, 
however, inconclusive. 
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Figure 5.6: The detected source density as a function of galactic latitude. The source 
density has been corrected for the variable exposure and background. The blazar dis- 
tribution appears uniform, however, there is slight evidence for an excess of uniden- 
tified sources at medium (20° — 30°) latitudes. 

5.3 Source Distribution Models 



The first step towards understanding the extent to which unresolved members of this 
population of sources contributes to the IDGRB involves gleaning information from 
the distribution of detected source fluxes. As will be seen shortly, the number of 
sources (A^) with a fiux greater than some value (S) is most often given by a power 
law and because power laws are most often displayed on log — log plots, the flux 
distribution A^(> S) is often called a logN — logS relation. 

If an observer sits at the center of an infinite distribution of standard candles of 
strength Lq, the observed flux of these candles will depend on their distance as. 



(5.6) 



while the number of sources observed with a given flux or greater is given by the 
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volume of the sphere enclosed of radius R, 

N{> S) = —pR^ = kS-^'^ . (5.7) 
3 

This simple form of the logN/logS relationship is called the Euclidean distribution. 
The Euclidean distribution applies just as well to any population of sources with a 
given luminosity distribution so long as this distribution does not evolve in time (and 
therefore space). One consequence of this type of distribution is that the integrated 
diffuse flux from such a distribution, 

I = — ^SdS = k'\S~'/^]^, = oo , (5.8) 

diverges in the limit of infinitely large distributions. Such distribution must be trun- 
cated at some minimum flux Smin in order to avoid this problem. This is commonly 
referred to as Olbers' paradox. 

Figure |5.7| shows the measured logN — logS for the detected blazars. Also shown 
on this plot is the source flux distribution for the unidentified sources. While the 
statistics are poor, there is no evidence for there being any detectable difference 
between these two groups of sources. The Kolmogorov-Smirnov probability that 
these source distributions are drawn from the same underlying distribution is shown 
in the figure to be 85%. Hereafter these distributions will be considered jointly. 

The simplest parameterization of the high latitude point source flux distribution 
is given by a truncated power law. 

diV ^ f A(5/5o)-T S > S^in 



and A, 'J, Smin are considered free parameters. Figure ^]8| shows the measured high 
latitude logN — logS function as well as the best fit power law parameterization. 
The best fit parameters of this fit are A = 33.1,7 = 1.45. This nearly Euclidean 
logN — logS function leads to divergent values of the total diffuse flux if it is integrated 
down to arbitrarily small source fluxes. It must break or roll over somewhere below 
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Figure 5.7: The measured logN/logS function for the identified AGN high latitude 
sources and the unidentified high latitude sources. The Kolmogorov- Smirnoff (KS) 
statistic is calculated and shows that the distributions are consistent with each other. 

the detection threshold. The roll over visible at the smallest fluxes is most likely due 
to the non uniform detection threshold due to the variable exposure. 

A more realistic parameterization allows sources to be distributed not only in 
space but also over a range of intrinsic luminosities. The number of sources per unit 
volume per unit luminosity, 

dN 

<'(r,L) = — , (5.10) 

is called the luminosity function. Knowledge of the luminosity function allows one to 
calculate the observed logN — logS function and can be integrated over all space and 
luminosities in order to calculate the total intensity contributed by a source class. 
In the case of blazars, the sources in question are located at cosmological distances 
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Figure 5.8: The measured log N/ log S function for the combined high latitude sources 
with I 6 |> 20°. The best fit power law logN — logS relationship is also shown. This 
best fit function is shghtly fiatter than Euclidean. The roll over below 1 x 10~^ is 
most likely caused by the nonuniform detection threshold. 

(i.e. redshifts z ~ 1). In order to obtain a logN — logS function from a luminosity 
function a cosmology must be assumed. For the purposes of this analysis a fiat 
universe will be assumed (Q = 1,A = or alternatively qQ — 0.5). Furthermore the 
Hubble constant, Hq, which relates recession velocity (most commonly measured by 
redshift,^;) to distance will be taken to be 75 km s~^ Mpc~^. 

Using this cosmological model allows a source fiux S to be calculated given a 
luminosity L, 



^=(71)12' (5-11) 



.1 L 

where di is the luminosity distance which for our choice of cosmology is given by 
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( IWeedman 1986|) , 

2c 

di = , , (5.12) 

Ho{l + z-Jil + z)) 



and we have introduced the concept that the radiation emitted by a source may be 
beamed into a sohd angle 6. One further comphcation is introduced if one is not 
deahng with bolometric flux but rather some discreet spectral window. Because the 
frequency of a given photon depends on the redshift at which it is observed, one must 
adjust for the shifting of a spectral window with redshift. This correction is called 
the K-correction. Given power law energy spectrum with spectral index a the flux 
becomes, 



,1,L(1 + 
f 

We can now compute an observed logN—logS function given a luminosity function 
by substituting the relation dL = dS{z) from equation 5.13 into the integral, 

dN _ r^rna. dV dN (1 + zf-"" 

Is ~ io l^dVdL 9df ■ ^^'^ ^ 

where dV/dz reflects the fact that the volume element per unit redshift changes as 
a function of redshift as the proper distance changes with redshift. This relation 
can in turn be integrated in order to determine the contribution of the point sources 
described by a given luminosity function to the diffuse background. 



'Sources 



where Vl is the total surveyed solid angle and the factor ^ accounts for the fact that 
as sources are more tightly beamed there is a higher probability that the emission 
will be beamed away from the observer. Note that the beaming angle cancels out 
of this equation making the total contribution to the background insensitive to the 
beaming of the sources. Smin is determined by the minimum luminosity as well as by 
the maximum redshift. 
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5.3.1 Direct Luminosity Function Constraints 

The identified EGRET blazars constitute a set of objects for which we now have 
measured 7-ray fluxes as well as measured redshifts from optical measurements. This 
population can be used to constrain the luminosity function. This has been done 



the EGRET blazars was deduced from the distribution of detected blazars using the 
V/Vmax test. For a given detection threshold and a model cosmology, an accessible 
volume of space can be defined as the region of the universe within which a source of 
given luminosity is able to be detected by an instrument at earth. If a given source 
remains at constant luminosity throughout its lifetime, it is equally likely to be found 
in the nearby half of this volume as the distant half. As a result, the distribution of 
the statistic V/Vmax should have a mean of 0.5 in the limit of zero source evolution. If 
a source tends to get brighter with time it will more likely be found in the near half of 
the total volume and as a result the mean V/Vmax would be smaller than 0.5. Chiang 
et al. measure V/Vmax to be 0.7 indicating that 7-ray blazars were more luminous 
in the past. This hypothesis is known as pure luminosity evolution. The alternative 
hypothesis is that the blazars were more numerous in the past which is known as pure 
density evolution. 

The evolution was extracted by assuming functional form for the evolution. 



and adjusting the free parameter /3 until V/Vmax was restored to 0.5. This was 
achieved for /3 = 2.6 ± 0.3. 

Having measured the evolution, the sources were all de-evolved into the same 
epoch {z = 0) and the luminosity function was determined by fitting a broken power 
law model to the luminosity function. 



by Chiang et al. ( [Chiang et al. 1995[ ). In this work, the luminosity evolution of 




(5.16) 




B 



(5.17) 
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where Lb is the break luminosity determined to be 10^^ erg s ^ and the power law 
indices were constrained by the data to be, 

71 = 2.9 ±2.0 

72 = 1.6 ±0.4 

(5.18) 

The large uncertainty in 71 reflects the small number of sources used to constraint 
the low luminosity end of this luminosity function. 

Integrating over this luminosity function allows one to estimate the contribution 
of the blazars to the diffuse background. In order to do this it is necessary to adopt a 
cut-off luminosity below which no sources exist. This is required in order to prevent 
divergence of the total flux from these sources. Because no such cut-off is directly 
observable in the 7-ray source distributions, the best that can be done is to adopt 
the lowest de-evolved luminosity as this cut-off. This choice results in a lower limit 
for the total diffuse contribution because the true cut-off must be smaller than this 
value. The resulting constraints are quite weak, 

/^e^ = {0.5 - 14.0} X 10^^ . (5.19) 

All that can be said is that these sources contribute signiflcantly to the diffuse back- 
ground. 



5.3.2 Scaled Radio Luminosity Function Constraints 

The large numbers of high latitude point sources that have been associated with 
flat spectrum radio quasars (FSRQ's) suggests that an appropriate test luminosity 
function to use would be the luminosity function for these sources which has been 
inferred from the large population of detected FSRQ's. Dunlop & Peacock ( Dunlop 
fc Peacock 1990| ) have analyzed several hundred flat spectrum quasars for which the 
redshifts are either directly measured or inferred from the relation between fC-band 
flux and redshift for faint galaxies. The data constrain the luminosity under a variety 
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of parameterizations quite well to redsliifts of about 2. Beyond this there is some 
indication of a redshift cutoff at 2 ~ 5. If one again adopts the assumption of pure 
luminosity evolution, the luminosity function can be parameterized as, 

MU, z) = lo-'-'nij^/-'' + (^/'T' , (5.20) 

where, 

logio Lciz) = 25.26 + l.lSz - 0.28z^ , (5.21) 

and the units of radio luminosity Lr and co-moving density pr are, respectively, 
W Hz~^sr~^ and Mpc~^(unit interval of log^gLr)^^. 

In order to understand whether this luminosity applies to the 7-ray luminosity of 
AGN's, it must be ascertained whether the radio and 7-ray fluxes correlate. Sala- 
mon & Stecker ( |Salamon fc Stecker 199^ ) report a strong correlation between radio 
luminosity and 7-ray luminosity. However, the correlation they observe is between 
luminosities is strongly biased by their choice to correlate luminosities rather than 
fluxes. Mattox ( pVIattox 1995| ) has correlated fluxes using a larger sample of marginal 
EGRET detected blazars using a rank order statistic and finds a significant correla- 
tion. However, recent work by ( |Muecke et al. 1996| ) finds no evidence for a correlation 
between concurrently measured 7-ray and radio flux. A possible explanation for the 
weakness in the correlation is due to the fact that radio emission is believed to origi- 
nate both in the relativistic jet as well in radio lobes in which decelerating electrons 
can create large radio fluxes. While the core radio flux may well be correlated with 
the 7-ray emission, the radio lobe emission should not be and depending on the rela- 
tive strengths of these two radio components the correlation could be made weaker or 
stronger. Furthermore the large time variability of these sources makes the choice of 
what fluxes to correlate somewhat ambiguous. While the evidence of this correlation 
is weak at best, it will be taken as an ansatz which will allow us to proceed. More 
correlated VLBI radio observations (e.g. [Grandi et al. 19^ , [Maraschi et al. 1994D 
will help to resolve this issue. 

The second challenge to the direct extension of the RLF to the 7-ray blazars 



CHAPTER 5. HIGH LATITUDE POINT SOURCES 



107 



is the relatively few radio loud AGN's that have been observed to emit in the 7- 
ray spectral region. In order to produce such intense sources of high energy 7-rays 
the AGN emission must be beamed in order to avoid 77 attenuation which would 
prevent the 7-rays from escaping the source. Salamon & Stecker thus propose that 
the 7-ray emission is beamed into a smaller solid angle than is the radio emission. 
This leads to the consequence that the earth can lie within the radio cone and outside 
the 7-ray cone to explain the discrepancy in the counting statistics. Dermer ( permer 
1995|) has provided a theoretical justification for this assumption by demonstrating 
that synchrotron self Compton radiation will indeed be more tightly beamed in a 
relativistic jet than would synchrotron radiation. 

Taken together, these two assumptions provide a two parameter prescription for 
defining a 7-ray luminosity function (GLF), 

p^{L^, z) = (0^/0,) V((^r/10^), z) , (5.22) 

This luminosity can be used as outlined above to generate an observed luminosity 
function. Figure |5.9| shows a surface plot of this luminosity function to help visualize 
the various distributions. The surface has been truncated at the approximate EGRET 
detection threshold for reasonable values of ^ in order to illustrate the visible region 
of the universe. 

Figure ^.10| shows a family of logN/logS functions generated using this scaled 
FSRQ luminosity function for several choices of the parameters ^, j3. The effect of 
decreasing ^ is to increase the 7-ray intensity for a given radio flux. As a result, 
a greater fraction of the universe becomes detectable and the logN — logS function 
shifts to higher values. There is a characteristic knee in the source distributions 
in this parameterization. The other result of changing the parameter is to shift 
the position of this knee. The second parameter (3 merely shifts the logN — logS 
distribution vertically. As a result it is possible to adjust both parameters so as to 
keep the number of sources above some detection threshold constant as the slope at 
this threshold is adjusted. In this way, the detected sources restrict the acceptable 
parameter space to a one dimensional band in ^, (3 space. The shape of the detected 
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Figure 5.9: The 7-ray luminosity function for blazars modelled as discussed in this 
chapter. The surface plot shows the density of blazars (p), as a function of 7-ray 
luminosity (L^) and redshift (z). The superimposed contour plot shows the contours 
of source flux as seen by the EGRET telescope. The luminosity function is truncated 
at the EGRET detection threshold to demonstrate the region of (L^, z) space which 
is inaccessible to EGRET. 

logN — logS function above the detection threshold in principal provides an additional 
constraint but there is not much leverage in this range and the statistics are poor. 

A more useful constraint is provided by the observed redshift distribution of 
the detected blazars, shown in Figure p.ll| . The redshifts were compiled from the 
NASA/IPAC Extra-galactic Database (NED). Shown superimposed are the redshift 
distributions calculated using the scaled FSRQ luminosity function and assuming a 
detection threshold of 1 x 10"''. Such distributions were calculated for a range of 
values of ^. The parameter (3 does not effect the redshift distribution. As mentioned 
above, low values of ^ mean that more of the universe becomes visible to EGRET 
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Figure 5.10: A family of logN — logS functions generated using the scaled FSRQ 
luminosity function and several choices of the parameters ^,(3. Also shown are the 
detected point sources with | 6 |> 30°. 

leading to a larger mean redshift. 

The combination of these two constraints limits the parameter space available to 
9.0 < ^ < 11.0. The second parameter (3 is then adjusted correspondingly in order 
to conserve the number of detected sources. The resultant range is .0001 < (3 < 9.0. 
These limits will be revisited in detail next chapter. 

Again it is possible to translate these constraints into contributions to the IDGRB. 
Integrating over the luminosity functions using the scaled radio cut-off luminosity as 
the 7-ray luminosity cut-ofT yields, 



/^Gjv = {0.3 - 2.6} X 10"^ 



(5.23) 
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Figure 5.11: The measured redshift distribution for the detected blazars. Also shown 
are the redshift distributions calculated using the scaled FSRQ luminosity function. 
The mean redshift of the distribution constrains the acceptable range of ^. 

This is slightly more restrictive than the constraints due solely to the 7-ray data but 
not sufficient to determine whether the bulk of the IDGRB is due to unresolved point 
sources. 



5.4 Discussion 

The summed EGRET data from the first 3| years of the mission have been reana- 
lyzed using the more stringent albedo cuts motivated by Chapter 3. 55 sources were 
detected with cr > 4. Of these 37 of these sources have been identified with FSRQ's. 
Another 17 are unidentified and the last is the Large Magellanic Cloud. At least one 
of these sources is revealed by the ehmination of the albedo background which had 
obscured it in the P12 analysis. 

A series of tests have attempted to ascertain whether the unidentified should 
be grouped with the identified blazars. Results of the analysis of their spatial dis- 
tribution, time variability, and source distribution are inconclusive. There may be 
an indication that there is a galactic, time invariant component to the unidentified 
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sources. This has already been suggested by ([Kanbach et al. 1996|) . If these sources 
are indeed blazars the question of why no radio counterparts are detected must be 
answered. The time variabihty of these sources may provide the answer. It may be 
that these sources were observed in a relatively radio quiet state. There is, however, 
evidence that these sources are less variable in the radio than in the 7-ray band. This 
makes such an explanation seem less likely. Without strong evidence to the contrary, 
the unidentified sources will be treated as being due to the same blazar population 
for the rest of this discussion. 

The observed sources place constraints on the luminosity function which models 
the distribution of these sources in space and luminosity. The work of ( [Chiang et al. 



1995D has shown that while this is possible, the resulting luminosity function con- 



straints are not very sensitive to low luminosity sources due to the small number of 
detections. 

This situation can be compared to the situation in the X-ray band in which the 
Einstein Observatory discovered many hundred high latitude sources. A uniform 
sample of 427 of these sources were associated with spectroscopically identified AGN 
and used to constrain a luminosity function ( [Maccacaro et al. 1991 ). These sources 
were sufficient to not only extract the evolution parameters but also to constrain 
the luminosity function over four orders of magnitude of luminosity. This data was 
sufficient to resolve not only a break in the luminosity function power law but also 
to show that pure density evolution was not consistent with the data. As a result 
the X-ray Seyferts have now been determined to be the dominant source of the hard 
X-ray background ( [Zdziarski et al. 1995| ). 

Adopting a scaled radio luminosity is problematic as well. This resulting correla- 
tion between radio and 7-ray flux is not seen very strongly. While it is possible to fit 
this luminosity function to the data assuming a smaller beaming angle for the 7-ray 
sources, the total contribution to the IDGRB is still not well determined. 

Because of the poor statistics of the detected blazars it is necessary to attempt to 
constrain the luminosity function by other means. In the following chapter, indirect 
constraints on the blazar luminosity function will be evaluated using the technique of 
fluctuation analysis. 



Chapter 6 



Fluctuation Analysis 



The observations described in the preceding chapters have estabhshed two results: 
an isotropic diffuse 7-ray flux of 1.33 x 10~^ is present, and a distribution of high 
latitude, isotropically distributed sources has been measured. The logical question 
that follows is , "Do unresolved members of this source population constitute the the 
IDGRB or is there some other truly diffuse source of this radiation?". As shown in 
the preceding chapter, the detections alone are not sufficient to answer this question 
directly. 

Similar problems confronted scientists working on the origin of the X-ray back- 
ground. While it is now generally accepted that this background is due to unresolved 
Seyfert galaxies, this was not clear from such missions as HEAO-1. One technique 
that was derived in order to deal with this problem in the X-ray band was the use of 
fluctuation analysis. 

Fluctuation analysis is a technique to search for the signature of excess variance 
in an ensemble of observations of a quantity: in this case the intensity of the IDGRB. 
This excess variance is used as a tracer for the inevitable random clumping of discreet 
sources. The method was first developed for radio data ( [Scheuer 19571 ) and has been 
used extensively in the analysis of the X-ray background ( [Shafer 1983| , |Hamilton| 
1987|) . The approach presented here is directly adapted from the analysis used by 
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6.0.1 Distribution of Surface Brightness 

The goal of this section is to derive the distribution of surface brightness, P{I) in 
the presence of a population of point sources as seen by an instrument with a known 
point spread function. Given a test log N/ log S relationship (^), we first determine 
the probability of picking a source at random from the sky and having it contribute 
an intensity / to a pixel subtending a solid angle 0. ( Note that / is a measure of 
surface brightness whereas 5 is a measure of source flux.) 

m) = ^ , (6.1) 

where H = J n{l)dl represents the total number of sources which contribute to the 
intensity in the pixel and n{I) represents the number which contribute and intensity 
between / and / + dl. 

Knowledge of the point spread function allows us to obtain n{I) in the general 
case in which sources are not perfectly localized. 

r dN 

<^)= / -r^m,oc,S))d^, (6.2) 

where S'(/, a, 5) is the flux of a source located at a, 5 that contributes an intensity / 
to a pixel at ao, ^o, 

^) ^ 7 PQPfi^j; A ^^o ' ^^-^^ 

and PSF{a, 6, aQ, Sq) is the appropriate EGRET point spread function. 

The general effect of an extended point spread function is to allow more sources 
to each contribute a smaller intensity to a given pixel. As a result wider functions 
tend to wash out fluctuations. 

We now have a method for analytically determining the probability of obtaining 
an intensity measurement of / if one source is near enough to a pixel to contribute to 
its intensity. In order to find the probability of making such a measurement with an 
arbitrary number of sources allowed to contribute, we proceed by adding one source 
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at a time. The probability of two sources jointly contributing an intensity / to a given 
pixel is, 

P2(/) = Pi{I - I')Pi{I')dI' . (6.4) 

Jo 

This is simply a convolution integral which is more easily evaluated in Fourier space. 

nP2{i)) = Hpmnpm = Hpm' , (e.s) 

where J^{P) is the Fourier transform of the distribution P. It follows that the prob- 
ability of measuring an intensity / given n contributing sources is given by, 

Pn{i) = T-\nPi{i)r) ■ (6.6) 

The probability of n sources being near enough to contribute to the intensity of 
a given pixel is determined by the Poisson probability with mean /x. Putting all this 
together we can now calculate the distribution of intensities in a given pixel with an 
arbitrary number of sources allowed to contribute. 

oo -At n 

ni) = E —^^-^((^(n(/) //.))-) . (6.7) 


Exploiting the linearity of the Fourier transform allows us to simplify this sum, 

l2-T{n{I)/i^))- 







= e-'^^-i(e^("(^») 

= ^-i(e(^("(^))-^)) . (6.8) 
P{I) describes the variation in surface brightness across the sky due to the set of 



point sources. Figure |6.1| shows a series of P{I) distributions for various logN/logS 
relations. In each case the total diffuse intensity is constant but the fraction of that 
total which is due to unresolved point sources increases from 10% to 100%. The 
increasing variance in these distributions is due to the increasing numbers of point 
sources that contribute to these distributions. It is this signal that indicates the 
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presence of unresolved point sources. 




UpMuft'' el""' s"' si"'! 



Figure 6.1: P{I) for three different logN/logS models with increasing point source 
contributions and decreasing true diffuse intensity. Point source contributions of 10% 
(solid line), 50% (dashed line), and 100% (dotted line) are shown. 

Because what is actually measured in a given pixel on the sky is a number of 
observed 7-rays and not a continuous variable /, the next step is to calculate the 
distribution of observed counts, n.^, in a given pixel in the sky for which the galactic 
background intensity is modelled to be Gij and to which the there is exposure Sij. 
This is simply the distribution of surface brightness P{I) convolved with the Poisson 
probability of measuring riij photons in that pixel. 



It should be noted that the distribution Pc{nij)ij bears an extra set of subscripts 
to denote the fact that this distribution only applies to a particular pixel. This 
distribution must be calculated for each point in the sky to account for the varying 
exposure and galactic foreground emission. 




(6.9) 
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6.0.2 Fluctuation Constraints 

The above discussion presents a methodology for obtaining a distribution of counts 
in a set of pixels given a test logN/logS relationship. We now turn our attention 
to solving the inversion problem that is presented by the actual data. We wish to 
constrain model logN/logS relationships given a set of measured counts in each pixel, 
{n°j}. We use the technique of maximum likelihood for this purpose. If the logN/logS 
relationship is parameterized by a set of parameters {ofc} then the likelihood is defined 
as, 

LiiKh = UMK; ■ (6.10) 

ij 

Again it is more useful to deal with the test statistic Ts, 

Ts = 2/nL({n° }; {a J) = 5]2Zn(PcK; {a J).,) . (6.11) 

This quantity can then be maximized with respect to the parameters of the logN/ logS 
relation in order to find the best fit. 

The advantage of this method is that the behavior of Ts in the null hypothesis 
is known. As discussed in Appendix A, Ts is distributed as about the best fit 
value, where v is the number of free parameters ( [Eadie et al. 19711) . This allows one 



to calculate confidence regions for the parameters at- 

Maximizing Ts reveals the most likely logN/logS model given the data but it 
does not ensure that the absolute fit to the data is an acceptable one. This is only 
the case if the parameterization includes the true distribution. Because each data 
point is selected from a different distribution, it is not possible to simply compare the 
collection of data points to any one distribution to test the absolute fit. Instead it is 
necessary to construct the distribution of the following quantity. 



P{x')ijdx' . (6.12) 



This new statistic has the property that it is distributed uniformly on the interval 
y = {0, 1} if the values Xij are distributed as P{x)ij. This can be seen easily by 
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changing variables, 



P'{y^J) = '^n^^3)^3 (6-13) 



r>( \ \-^i-J 



= 1 



In our case the independent integral is over counts and so the integral becomes a sum, 

Vij = Pc{n; {ak})ij + ^Pc{n1j\ , (6-14) 

^ 

where the last term is divided by two in order to deal with the discontinuity at the end 
point of the integral. Comparing the distribution of {vij} with a uniform distribution 
gives a direct measurement of whether the absolute fit to the model is an acceptable 
one and as a result whether the parameterization is valid. 



6.0.3 Additional Constraints 

The analysis that has been described thus far does not include the constraint that any 
choice of logN/logS relationship that is used must match the observed distribution 
above the detection threshold. Because this constraint should be independent of the 
fluctuation constraint, the joint likelihood is multiplicative. If we detect No^s sources 
above the detection threshold, 

L{{nij}: Nobs] oik) ^ L{{n'^j};ak) X -Pdetections(^o6s; Qifc) ■ (6.15) 

The probability of observing Nobs sources above the detection threshold given a 
logN /logsS relationship which predicts N^^, sources above the detection threshold, is 

Pdetectionsi^ak) — 77 — ) (6.16) 
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So the likelihood function that describes both constraints is 

Ts = 2lnL{n, N; Ok) = 2{Y, ln{Pc{ni,; ak)ij)) - 2N^^ + 2NobMNaJ , (6.17) 

where we have dropped model independent terms. This function is optimized with 
respect to the free parameters a^. 

6.1 Implementation 

The implementation of the above technique will now be outlined. The EGRET pho- 
tons are binned in celestial coordinate bins such that each bin has equal solid area. 
This is done by increasing the latitude interval with increasing latitude by a factor 
l/cos{6). For a given choice of input logN — logS parameters, the distribution n{I) 
is calculated by numerical integration using equation 6.3. In general the tail of the 
point spread function causes n(/) to get very large at small values of /. In other 
words, a great many distant sources contribute small intensities to a pixel. Compu- 
tationally it is impossible to account for all of these sources explicitly, however, their 
effect is indistinguishable from a true diffuse flux because of the large numbers of 
these sources. Consequently, all sources that contribute an intensity smaller than Jq 
are converted into a set of sources with intensity equal Jq such that the total intensity 
these sources contributes remains unchanged. If Jq is chosen to be sufficiently small, 
the fluctuation analysis should remain unaffected by this computational approxima- 
tion. In the following analysis Jq is chosen to be 1 x 10~^ which is four orders of 
magnitude smaller than the diffuse intensity. The analysis is completely insensitive 
to the form of n(J) at this point. 

The distribution n{I) is then used to calculate P{I) using equation 6.9. The 
Numerical Recipes subroutine df ourl is used to perform the necessary Fourier trans- 
forms (press et al. 199^ ). A discreet imaginary array 2^^ elements long was used to 



define n(/) and P{I)- The step size used was 1 x 10 . The accuracy of the transform 
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and inverse transform is gauged by checking that the quantity, 

7=j2Pim, (6-18) 

i=0 

is indeed consistent with the input value of the diffuse intensity due to point sources. 

For each pixel, the distribution Pc{n)ij was constructed. The quantities Eij, Gij 
were read from the exposure and background maps respectively. The mean counts in 
a given pixel is given by, 

n = (/iDGRB + Gij)Sij . (6.19) 
The integral over / in equation 6.10 is carried out up to a maximum value of 




(6.20) 



The severity of this approximation is gauged again by check the integral, 

k='£ Pc{n) (6.21) 

n=0 

If this value deviates from 1.0 by more than a few percent then Umax must be increased 
because there is a significant probability that n exceeds nmax- 

Finally the Ts is calculated using equation 6.19. The whole process is repeated for 
a different set of logN/logS parameters and the best fit set of parameters is selected. 
Confidence regions are ascertained using the contours of Ts — {Ts)raasi- 95% confidence 
regions imply Ts drops of 5.8 for 2 free parameter models and 7.7 for 3 free parameter 
models. 

A consistency check is then performed to ensure that the best fit model is a 
absolute good fit. This is done by constructing the distribution of and looking at 
the /v probability. If it is not close to unity, the parameterization is inappropriate 
to the data. 
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6.1.1 Spatial Scales 



The spatial scale at which one searches for fluctuations is of critical importance to 
this analysis. It is important to optimize the sensitivity of the analysis to changes in 
the point source distribution while minimizing any contaminating signal from other 
sources of emission. 

It is self evident that the fluctuation signal is washed out on large spatial scales 
because too many sources are allowed to contribute to the intensity in a given pixel. 
However, the spatial scale of the test bins cannot be arbitrarily reduced as the analysis 
starts to become dominated by the Poisson noise associated with counting a small 
number of photons. The most sensitive scale is therefore a compromise between these 
considerations. 

This issue was investigated for the EGRET data using simulated data. Simulated 
EGRET counts maps were generated using the galactic diffuse model and the com- 
bined exposure maps (see Appendix A). A model of the sky is constructed including 
diffuse galactic emission, diffuse isotropic emission, and a list of point sources. This 
set of point sources was sampled randomly from an input Euclidean point source 
distribution, distributed randomly across the sky, and convolved with the EGRET 
point spread function. The resulting maps were Poisson sampled in order to generate 
a typical EGRET observation of this model sky consistent with the EGRET exposure 
map. This data can then be treated in the same manner as the flight data. 

For the purposes of this simulation a Euclidean logN/logS relationship was used. 
This logN/logS relation was then extended below the detection threshold and trun- 
cated at a specific value. 



Another way of parameterizing the same set of functions is using the integrated 
contribution to the diffuse flux. 




-1.5, 



srcs sr 




(6.22) 



(6.23) 



-^IDGRB 
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In the simulation Smin was chosen to be 2 x 10~^. This choice gives e = 10%. The 
data were then analyzed using a set of different values for e and using a variety of 
pixel sizes. The goal of this simulation is twofold: to provide a consistency check of 
the fluctuation calculation, and to find the spatial scale which shows the maximum 
sensitivity. 



5 




Figure 6.2: Likelihood curves as a function of e on several spatial scales. The input 
value of e was 0.1 which is consistent with the maximum likelihood value for all spatial 
scales except for the smallest. The fluctuation analysis is shown to be most sensitive 
at 3.5°. 

Figure |6.2| shows the results of this simulation. The first important result is that 
the likelihood (T5) is maximized at the input e. This verifies that the calculation of 
P(/) and Pc{n)ij is accurate. While this is not true in the case of the smallest spatial 
scale (1.5° x 1.5° pixels), the likelihood function at this scale is so fiat for e > 0.1 that 
this does not constitute a significant discrepancy. 

The sensitivity of the analysis is measured by the steepness of the likelihood 
function around the best fit value of e. This steepness reaches a maximum in the case 
of 3.5° X 3.5° pixels. At this scale, the average bin contains ~ 50 photons. Smaller 
pixels start to become dominated by Poisson fluctuations and larger pixels contain 
too many sources. 
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Figure 6.3: The distribution of y as defined in equation 6.16 for the best fit model 
to the simulated data. The value of x^/^ shows that the distribution is entirely 
consistent with a uniform distribution. 

The distribution of y provides the final consistency check of this method. Fig- 
ure |6.3| shows the distribution of y for Smin = 2 x 10~^ as well as the x^/ ^ assuming a 



this data is drawn from a uniform distribution. The result that x^/ ^ is nearly unity is 
evidence that the calculated distributions Pc{n)ij are accurate. This also implies that 
higher order effects that are not being taken into account are not seriously impairing 
the analysis. The principal complication that is not treated in this analysis is the 
effect of the non-uniform exposure on the sky on P{I)- In calculating P{I) exposure 
is assumed to be uniform and the detection threshold is not allowed to vary across the 
sky. Because the point spread function reduces the effects of distant sources this ef- 
fect is not critical as evidenced by the relatively good agreement between simulations 
using the actual EGRET exposure maps and the calculated distributions. 



6.1.2 Foreground Emission Uncertainty 

The above simulations are carried out using a perfect model for the the foreground 
galactic diffuse 7-ray emission. When analyzing the flight data, the potential effects 
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Figure 6.4: The distribution of y as defined in equation 6.16 for the best fit model to 
the flight data using the galactic diffuse model derived from radio data. The value of 
X^/z/ shows that the distribution is inconsistent with a uniform distribution and thus 
that the diffuse model is inconsistent with the data. 

The diffuse gas model was discussed in the previous chapter. The large scale 
excesses at high latitude cannot be explained by point source fluctuations. This 
point is emphasized in Figure |6.4| . This figure shows the distribution of y for an 
input logN — logS in which the point sources make up the entire IDGRB. This model 
represents the best fit to the data and yet it is a demonstrably poor fit to the data as 
evidenced by the bowing of the y distribution. The data shows more variance than can 
be accounted for using point sources alone. In order to perform meaningful fluctuation 
analysis, it was necessary to use a diffuse model that gives isotropic residuals on large 
spatial scales. Such a model is generated by the maximum likelihood fit to the diffuse 
background described in the previous chapter. 

Using the data to constrain the diffuse model introduces potential biases to the 
analysis that must be guarded against. The process of fitting the diffuse model to 
the data can erase real fluctuations on the spatial scale of the likelihood fit. Thus if 
the diffuse parameters were fit on the scale of a few square degrees, the fluctuations 
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would be absorbed into the foreground diffuse model and the fluctuation analysis 
would be biased toward smaller point source contributions. Figure |6.5| shows the 
relative fluctuations using 15° radius bins compared to the fluctuations on a 3.5° 
scale. The relative absence of a signature due to point sources on the larger scale 
ensures that no significant point source signal is being absorbed into the galactic 
foreground model. 




Figure 6.5: P{I) distributions calculated using (15°)^ pixels (dashed curve) and for 
(2.5°)^ pixels (solid curve). A Euclidean extrapolation of the detected source dis- 
tribution was used to calculate these curves. The sharpness of the dashed curve 
demonstrates the absence of significant fluctuations on large spatial scales. 

A more important bias is caused by small scale fluctuations in the galactic fore- 
ground diffuse model that are unaccounted for when the likelihood scaling is done. 
If the unmodelled source of emission that is removed through likelihood analysis ex- 
hibits spatial structure on the scale of the fluctuations, this structure could bias the 
analysis toward models containing higher point source contributions. 

If it is assumed that the additional non-isotropic diffuse radiation beyond what is 
directly modelled is produced by some interaction between cosmic rays and photons 
or gas, then its spatial distribution should be similar to that of the modelled gas. 
Under this assumption the bias caused by small scale fluctuations in this radiation 
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can be directly estimated. The integral in equation 6.10 becomes a double integral, 

poo poo 

Pc{n^,)r, = / ^ T P{I)P\G,,)dG,,dI , (6.24) 



=0 riij. 



where P'iGij) expresses the uncertainty in the foreground emission, Gij, on the spatial 
scales under investigation. The galactic diffuse model can be used to directly measure 
P'iGij). This is done by calculating the residual intensity between a a pixel of solid 
angle 3.5° x 3.5° and the mean intensity within 15° of this pixel for a set of independent 
pixels across the sky. This distribution is well approximated by a Gaussian with a 
relative error of 14%. This distribution is used in calculating Pc{nij)ij in order to 
account for this bias. 



6.2 Results 



The data analyzed using these techniques is shown in Figure 3^. This plot shows 
the residual intensity from the likelihood fit to the diffuse model. Regions that are 
excluded due to their proximity to detected points sources or the galactic plane are 
shown in black. 

In order to understand how the various contributions to the likelihood function 
affect the constraints Figure shows the likelihood contours of the various compo- 
nents of the likelihood function separately as well as jointly. The plot shows likelihood 
contours in the space of power law logN — logS relations parameterized by A, 7 as 
defined by equation 6.24. For the purpose of this plot, Smin was fixed at 1 x 10^^°. 

The contributions of the fiuctuations to the overall constraints is clearly visible 
in this plot. The constraint imposed by the detected point sources is manifestly 
insensitive to the slope of the logN — logS function below the detection threshold. 
As a result this constraint is linear in the function parameter space. The fiuctuation 
analysis is sensitive to both parameters but is more sensitive to the overall slope of 
the logN — logS function. The joint constraint, shown as a solid line in the contour 
plot, is thus bounded in both dimensions. 

The minimum point source contribution which is allowed by the data is determined 
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Figure 6.6: The residuals from tiie likelihood fit to the diffuse model. Regions that 
are excluded from the analysis are shown in black. These are regions near strong 
point sources and the galactic plane. 



by exploring the three dimensional parameter space in which the added dimension 
represents the position of the break in the logN — logS function. The 95% confidence 
region in this parameter space achieves its minimum point source contribution at the 
point given by (A, 7, = (20, .6, < 1 x 10~^). The analysis is not sensitive to 

values of Smin < 1 X 10"^. This is because there are a sufficiently large number of 
sources below this flux to wash out their fluctuations. 

The constraints on steeper logN — logS functions are more interesting. It should 
be noted that Euclidean logN — logS functions are strongly excluded as the maximum 
allowed value of 7 is 1.27. It is clear that that the fluctuation analysis indicates a 
flattening of the source distribution below the detection threshold. 

Any source distribution that approximates a power law is constrained according to 
the analysis presented above. If the scaled radio luminosity function for flat spectrum 
radio quasars is used to generate a logN — logS function as described last chapter, the 
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Figure 6.7: The likelihood contours in the space of power law source distributions. 
The dashed lines indicate the constraints imposed by the directly detected sources. 
Not surprisingly these constraints are only restrictive in one dimension as they are 
insensitive to the slope of the power law below the detection threshold. The dotted 
lines indicate the fluctuation constraints. These constraints are sensitive to both 
parameters but are more sensitive to the slope parameter. The joint 95% confidence 
interval is bounded by the solid line. 

result is not well approximated by a power law. This is because there is a characteristic 
knee in the distribution. For this reason this parameterization was treated separately. 
As described last chapter, this method leads to a two parameter model of the logN — 
logS function, where. 



Lr = lO^L^ 

f3 = {e./Orf 



(6.25) 
(6.26) 
(6.27) 
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The constraints on these two parameters are highly correlated so the likelihood con- 
tours in this parameterization are shown in Figure in the space of ^, B{^, j3) where 
B is chosen to minimize the correlation between the uncertainty in ^ and the uncer- 
tainty in B. Contours of (3 are shown overlayed on this plot. It is important to 
note that if the knee in the logN — logS function occurs at fluxes higher than the 
detected sources, the slope of the logN — logS function does not vary as one changes 
the parameters. At this point the parameterization becomes essentially a one dimen- 
sional power law model in which the normalization of the logN — logS relationship 
is the only free parameter. For this reason the constraints are unbounded at the low 
end. Physically what this represents is the fact that the scaling proportional constant 
between the 7-ray and radio luminosities can be made arbitrarily large as long as 
the 7-ray beaming factor is made sufficiently small to retain the correct number of 
observed sources. 

It is more instructive to display both these constraints by projecting the confidence 
interval back into the observed distribution space rather than the parameter space. 
Figure shows the fiuctuation constraints as a gray shaded area on a logN — 
logS plot. All the source distribution functions allowed by the fiuctuation analysis 
presented above lie within this region. It should be noted that these boundaries are 
curved because they represent the locus of tangent lines to the bounding logN — logS 
functions as a function of source fiux. The constituent logN — logS functions are 
nonetheless straight power law lines. The fiuctuation analysis allows one to constrain 
the logN — logS distribution an order of magnitude below the detection threshold. 
In this representation the fiattening of the logN — logS distribution suggested by the 
fiuctuation analysis is manifest. 

The constraints on the source distributions generated using a scaled radio luminos- 
ity are shown as a vertically hatched region. The fiattest allowed source distribution 
is comparatively steeper than the corresponding power law constraint. This is a result 
of the fact that the parameterization does not allow fiatter logN — logS functions. 
The luminosity function itself determines the lower bound on fiat distributions. The 
steeper models are, however, constrained by the fiuctuations. Here the upper bound 
is higher than for the power law constraints. The steepening of the source distribution 
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Figure 6.8: The likelihood contours in the space of scaled radio luminosity function 
source distributions. The 68% and 95% contours are shown as solid lines. Because 
the two parameters ^, P are highly correlated the parameter space is displayed as a 
function of a third variable B which minimizes this correlation. The contours of f3 
are overlayed in order to illustrate the range of this parameter allowed by the data. 
The data are unbounded for small values of ^ because the minimum slope of the 
logN — logS relationship is determined by the luminosity function. 

when parameterized in this way allows models with comparatively large numbers of 
sources below the detection threshold to remain consistent with the observed sources 
at the detection threshold. 

Also shown in Figure |6.9| are the family of source distribution models generated us- 
ing the pure luminosity evolution model of ( [Chiang et al. 19951) . The region spanned 



by these models is shown as the horizontally hatched region in Figure |6.9| . These 
models are steeper than the corresponding functions generated using a scaled radio 



CHAPTER 6. FLUCTUATION ANALYSIS 



130 




Flux (E> 100 MeV) (ph cm^'^s"') 



Figure 6.9: logN/logS constraints. The jagged line represents the flux distribution of 
the detected sources with a > 4. The light gray region shows the region within which 
the extrapolated power law logN/logS function must lie. The horizontally hatched 
region displays the constraints on this function if it obeys a scaled version of the 
radio luminosity function. The directly deconvolved 7-ray luminosity function under 
the assumption of pure luminosity evolution is constrained by the vertically hatched 
region. The black region indicates the areas in this space in which the integrated 
intensity from the point sources exceeds the measured IDGB. 

luminosity function. Since the source distributions generated by this model are es- 
sentially power laws, the gray areas on the plot represent the region of this space 
that is allowed for this model by the fluctuation analysis. There is however a zone of 
consistency between the two regions. 

The final constraint illustrated in Figure |6.9| is the limit imposed by the direct mea- 
surement of the IDGRB. The summed contribution of the unresolved point sources 
must not exceed the total flux of the observed sky. Because this constraint is on the 
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integral flux, 

-^Sources = / ^~Jq^^ ' (6.28) 

this region is not uniquely determined in this plot. The region shown in dark gray 
represents the region within which a power law extending directly from the detected 
sources at the detection threshold yields more than the observed IDGRB. Because 
most functions discussed here are near power laws it is a good approximation to the 
region excluded by direct observation of the diffuse. It should be noted that most of 
the family of curves predicted by the pure luminosity evolution model are excluded 
by direct observation of the diffuse background. 

6.3 Discussion 

The fluctuation analysis described here allows the next decade of point source flux to 
be probed without having to resolve the sources individually. The limitations of this 
technique stem principally from the difficulty encountered in trying to separate vari- 
ations in sky intensity caused by foreground emission from point source fluctuations. 
Furthermore, there is no way to treat the identifled and unidentifled sources inde- 
pendently because both contribute to the fluctuation analysis in an indistinguishable 
way. 

Despite these difficulties, this analysis is able to place interesting constraints on 
the sources below the detection threshold. Table 6.1 summarizes the constraints that 
result from this analysis. 



Table 6.1: Summary of Fluctuation Constraints 



Parameterization 


Parameter Constraints 


Diffuse Contribution 


Power Law 


7 = {0.61. L27} 
-4 = {69.2. 177.8} 


Srnin = : {6%. 100%} 

= 10-^ : {1%. 15%} 


Radio Luminosity Function 


e = {9.1, 10.4} 
/3 = {.005, .4} 


{21%, 120%} 


Pure Luminosity Evolution 


71 = {1.0,1.9} 

72 = (2.0, 2.6} 


{5%, 14%} 
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The weakest constraints are those provided by using only the direct hmits on 
the logN — logS relationship without using the measured redshifts to constrain a 
luminosity function The fluctuation analysis has revealed the existence of at least 
an additional 100 point sources in the decade of flux below the EGRET detection 
threshold. These sources make up 4% of the measured IDGRB. Furthermore, the 
analysis described above has limited the number of sources in this decade to be less 
than ~ 300 making up 15% of the IDGRB. In the absence of a means of determining 
the cut-off in the source distribution, the fluctuation analysis permits a large range of 
models in which the entire IDGRB is composed of unresolved point sources. By itself 
it thus provides only comparatively weak constraints of the point source contribution 
to the IDGRB. 

If these source distributions are put in the context of a cosmological model through 
the introduction of a luminosity function, more restrictive constraints are possible. 
Two luminosity functions have been explored in the last two chapters. In each case, 
the directly observed sources provide the most restrictive lower limits on the point 
source contribution while the fluctuation analysis provides the more restrictive upper 
limit. 

The first model described uses the inferred luminosity function from the FSRQ's 
measured at radio wavelengths. Under the ansatz that there is some scaling law 
between the radio fiux and the 7-ray flux this luminosity function can be used to 
describe the 7-ray blazars provided the 7-ray emission is assumed to be beamed into 
a smaller solid angle than the radio emission. 

There arc three independent constraints on such a luminosity function:thc source 
flux distribution of the detected sources, the redshift distribution of the detected 
sources and the fluctuation constraint on the sources below the detection threshold. 
As shown in the previous chapter, the redshift distribution constraint restricts the 
scahng parameter ^ to the range 9.0 < ^ < 11.0. Models with smaller values of ^ tend 
to predict sources at a higher mean redshift the is observed in the detected blazars 
and conversely models with larger values of ^ tend to have distributions with smaller 
than detected mean redshift. 

The detected source distribution imposes a a relationship between the two model 
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parameters ^, p. In order to predict the correct number of sources above the detection 
threshold of ~ 1 x 10~^ as ^ decreases, the beaming angle for blazar 7-ray emission 
relative to the radio emission must be similarly reduced. The result is that the 
luminosity functions tend to predict larger numbers of sources with high 7-ray flux. 
In order to make the number of sources above the detection threshold remain constant, 
these additional sources must be assumed to be beamed away from the earth which 
places a tighter constraint on their beaming angle. 

The fluctuation analysis is insensitive to the lower bound on ^. This is because 
the shape of the logN — logS relationship is invariant to changes in ^ once the knee 
in the distribution has been shifted to fluxes larger than the detection threshold. On 
the other hand, large values of ^ are efficiently constrained using this analysis because 
as the knee shifts through the detection threshold, the resulting number of sources 
predicted by the luminosity function below the detection threshold varies rapidly with 
increasing ^. The analysis presented above shows that the upper limit on ^ can set 
at ^ < 10.4 using this analysis. 

The real power of this approach lies in the fact that the constraints on the high flux 
tail of the luminosity function can be used to determine the two scaling parameters 
allowing the comparatively well constrained radio luminosity function to be used. This 
luminosity function can then be integrated over the entire distribution which extends 
many orders of magnitude below the detection threshold in order to estimate the total 
blazar contribution to the IDGRB. The results of this integration are shown in Table 
6.1 for consolidated constraints provided by three independent measurements. These 
joint constraints cannot exclude models in which unresolved blazars contribute 100% 
of the IDGRB, however, this contribution can be as low as 20%. Significant amounts 
of true diffuse emission cannot be ruled out on the basis of this analysis. 

The second luminosity function investigated was determined through the direct 
measurement of the 7-ray blazar luminosity evolution. This model assumed pure 
luminosity evolution was responsible for the departure of the value of V/Vmax from 
0.5. Under this assumption a broken power law luminosity function was constrained 
using the de-evolved source luminosities of the detected blazars with measured red- 
shifts. This luminosity is also parameterized by two free parameters: the power 
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law index above and below a break luminosity. As shown in the preceding chapter, 
the power law index below the break luminosity is poorly constrained by the detected 
sources because very few low luminosity sources are detected by EGRET. As reported 
previously this index is restricted to the range 1.0 < 71 < 4.9. The point source con- 
tribution to the diffuse background is, however, dominated by low luminosity sources 
and result reflects this uncertainty. 

The fluctuation analysis has again added an upper limit to the steepness of the 
logN — logS distributions that result from this model. In general the family of 
logN — logS distributions which result from this model tend to be steeper than those 
resulting from the luminosity function previously discussed. Most of the parameter 
space allowed by the detected sources is excluded when one considers the fluctuation 
constraints. The fluctuation constraints limit the parameter 71 to the range {1.0, 1.9}. 
This reflects the fact that the fluctuations indicate a flattening of the logN — logS 
distribution below the detection threshold. 

As mentioned in the preceding chapter, a minimum luminosity must be assumed 
in order to evaluate the detection threshold. Because this model does not assume a 
correlation between the radio and 7-ray flux, the cut-off luminosity is determined by 
the minimum observed luminosity. As a result the integrated intensity due to unre- 
solved sources in this model is really an upper limit to the point source contribution 
as the cut-off luminosity may be far smaller than the minimum detected luminosity. 
As shown in Table 6.1, the integrated point source contribution for this luminosity 
function after applying the fluctuation constraints is e = {5%, 14%}. These are much 
smaller than the estimates made using the detections alone which predicted that the 
entire IDGRB would be made up of sources in the next 3 decades of flux below the 
detection threshold. These model were biased by the steep logN — logS distribution 
above the detection threshold. This may be the result of the biasing of sources near 
the detection threshold toward higher detected sources. This is caused by sources just 
below the detection threshold contributing to the flux of detected sources (|Schmitt 



1986| ). This effect is demonstrated in simulations described in Chapter 7. 



There remain many uncertainties associated with the point source contribution 
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to the IDGRB. The preceding analysis suggests that 7-ray blazars contribute signifi- 
cantly to the diffuse extra-galactic 7-ray background, however, it is not possible to say 
that they constitute the entire IDGRB. Several outstanding questions remain. These 
include: the nature of the unidentified sources, the correlation or non-correlation be- 
tween radio and 7-ray fluxes, the validity of pure luminosity evolution, and the cutoff 
luminosity of the luminosity functions. Many of these outstanding issues will be ad- 
dressed in the following chapter which describes what might be possible with a next 
generation 7-ray telescope. 



Chapter 7 

The Gamma-ray Large Area Space 
Telescope 



EGRET has provided an unprecedented view of the 7-ray universe. It has brought 
7-ray astronomy into the mainstream of astrophysical research. This discussion of 
the diffuse background as measured by EGRET highhghts some of the hmitations 
of this instrument. While some of these hmitations reflect the inherent limitations 
of satellite based 7-ray astronomy, others can be overcome in the future using the 
knowledge gained from the EGRET mission. 

This chapter will first outline some of the scientific issues pertaining the study of 
the study of the diffuse 7-ray background that are beyond the reach of the EGRET 
instrument and the requirements these issues make on a future 7-ray telescope. The 
rest of the chapter will be dedicated to the description of one such possible future 
instrument: the Gamma-ray Large Area Space Telescope (GLAST). The design of 



this telescope has been developed by a large collaboration of physicists ([Michelson 



et al. 1995| ) led by groups at Stanford University and the Stanford Linear Accelerator 



Center (SLAG). The design of this instrument will be summarized here as a prelude to 
a discussion of the scientific capabilities of such an instrument. A virtual calibration of 
this instrument has been carried out through Monte Garlo simulation by W.B.Atwood 
at SLAG. This work provides the basis for the evaluation of the instrument's scientific 
performance. 
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7.1 Unresolved Scientific Issues 

7.1.1 The Nature of the Isotropic Diffuse Gamma-ray Back- 
ground 

EGRET has confirmed the existence of a high energy isotropic diffuse 7-ray back- 
ground. The intensity of this emission has been measured with unprecedented ac- 
curacy. Considerable uncertainty about its intensity remains due to the remaining 
uncertainty associated with our hmited understanding of the foreground diffuse emis- 
sion. 

As explained in the preceding chapters, the physical explanation of this back- 
ground now centers on the radiation produced in relativistic jets emanating from 
massive black holes in active galaxies. Such sources have been positively detected 
by EGRET and unresolved members of this population undoubtedly contribute to 
the IDGRB. Because only the brightest examples of these sources are observed, there 
remains considerable uncertainty as to their total contribution to the total IDGRB. 
Fluctuation analysis has placed an upper limit to this contribution as being > 5% of 
the total but there remains the possibility that much of the observed intensity is the 
result of a truly diffuse source. In order to further constrain this contribution it will 
be necessary to achieve an increase in point source sensitivity in order to place more 
meaningful constraints on the 7-ray blazar luminosity function. 

In light of the difficulty of exactly accounting for the integrated blazar contribution 
to the IDGRB, a definitive discovery of a truly diffuse 7-ray background will most 
fikely be discovered through spectral signatures. As discussed in Chapter 1, the 
most intriguing possibility for such a discovery would result from the annihilation 
of a WIMP with some characteristic energy. This sort of particle would result in a 
spectral bump (in the case of an unstable WIMP) or a spectral fine (in the case of 
a stable WIMP) in the high energy diffuse spectrum. Interesting constraints on such 
theories could be obtained with improved counting statistics at high energies. 
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7.1.2 Understanding the Active Galactic Nucleus Accelera- 
tion Mechanism 

In addition to measuring the contribution of active galaxies to the IDGRB, a detailed 
understanding of the physics of the acceleration mechanism in these sources remains 
beyond the scope of EGRET. These sources are distributed over a wide range of 
redshifts (.003 < z < 2.28) and have can broadly be classified as blazars (a term 
encompassing BL Lacs, OVV and HPQ quasars). Many are associated with superlu- 
minal radio sources. The huge flux emitted in the 7-ray energy band as well as the 
time variability observed on time scales as short as a few days, all seem to point to 
massive black holes as the central acceleration engine. 

While there is general agreement that these extraordinary sources that the energy 
source is accretion onto a massive black hole, there is little agreement on the specific 
models used to explain the multiwavelength emission as well as the time variability. 
The latter points strongly to a beamed highly relativistic jet being involved but 
disagreement persists as to the composition of such a jet (electon-positron jet versus 
electron-proton jet) ( [Mannheim 1992 , Blandford fc Levinson 19"75D . 



The most powerful constraints on the physics of the acceleration mechanism arise 
from correlated multiwavelength studies that span the huge spectral space in which 
these sources emit ( [von Montigny et al. 1995]) . There is currently some controversy 
about the correlation of 7-ray flares with optical and radio flares ( [Wagner et al. 1995a , 



Wagner et al. 1995E| ). While the EGRET team has led multiwavelength campaigns in 



order to address particular sources with simultaneous observations, it is impossible to 
predict the occurrence of such events due to the fairly small and uncertain duty cycle 
of blazars. A large solid angle 7-ray monitor would allow much more powerful insight 
into the emission mechanisms as ground based observers in other spectral bands could 
be alerted to the occurrence of a blazar flare allowing correlated multiwavelength 
studies to be performed on a large population of flares. 
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7.1.3 What is the Nature of Unidentified Gamma-ray sources? 

Perhaps the most intriguing question to emerge from the EGRET mission is the 
question of the nature of the many unidentified point sources that have been identified 
based on their 7-ray emission. There are 80 sources reported in the second EGRET 
catalog with no counterparts in other wavelength bands. Of these 32 lie within 10° 
of the galactic plane and the remaining 48 are at higher latitudes. 

The identified sources fall into two categories: nearby pulsars on the galactic 
plane, and blazars and BL Lac's off the galactic plane. The only identified sources 
that do not fall into this category is the Large Magellanic Cloud (LMC). It remains 
an open question whether the unidentified sources represent further examples of these 
two classes for which the corresponding radio and optical counterparts have yet to be 
discovered or whether they constitute an entirely new class of object. 

The Geminga pulsar provides one clue as to the nature of the unidentified galactic 
plane sources. While this source is the second brightest point source in the 7-ray sky, 
it shows no evidence of a radio counterpart. This puzzle was resolved when EGRET 
detected pulsations consistent with the period of the X-ray pulsations discovered using 
Einstein Observatory data ( Palpern et al. 1992| ). It is now believed that the radio 
and 7-ray emissions are beamed into different solid angles allowing the 7-rays from a 
source to be detected without a radio counterpart ( [Romani et al. 1995| ). If this is the 
case then it is possible that some fraction of the unidentified sources on the galactic 
plane are simply pulsars whose radio emission is simply beamed away from the earth 
( [Yadigaroglu fc Romani 1995 ). A search for pulsations in these sources is underway 
(e.g. [Brazier et al. 1996| ). Furthermore, it is possible that a faint radio counterpart 
has not yet been detected in the radio band. Several searches of EGRET error boxes 
have been performed ([Nice et al. 1994|) . 

Another possible class of objects that may contribute to this population is that of 
Supernova Remnants (SNR). ([Sturner et al. 1995)) have found a statistical correlation 
between these objects and EGRET sources. An increase in resolution will greatly 
constrain this correlation. 

The high latitude unidentified sources have been discussed in some detail in Chap- 
ter 5. There is little observational evidence with which to distinguish these sources 
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from the blazars. It is possible that they exhibit a smaller degree of time variability 
but this is most likely due to a bias caused by truncation at the detection threshold. 
As presented in Chapter 5 there is some suggestion of a slight excess of sources near 
the galactic plane. This could be an indication that a fairly local population of galac- 
tic sources contributes to these detection but the evidence for this is statistically weak 
due to the small number of sources detected. An increase in resolution is critical in 
narrowing the list of candidate counterparts to these sources. 



7.1.4 Closing the Gap with Ground Based Gamma- ray As- 
tronomy 

7-rays with energies greater than ~ 300 GeV cause sufficiently large electromagnetic 
showers when they interact with the earth's atmosphere, that the Cerenkov light from 
these showers can be imaged by ground based 7-ray telescopes ( |Weekes et al. 19"89| ). 
While the rejection of cosmic ray backgrounds is challenging, these telescopes have 
large effective areas which allow them to overcome the statistical challenge presented 
by the falling spectrum of sources in this region of the spectrum. 

To date 4 7-ray sources have been detected by ground based sources. The Crab 
nebula remains the brightest source in the sky at TeV energies. The EGRET pulsar 
1706-44 has also been detected. The spectra of these sources are well joined to the 
EGRET data by synchrotron self-Compton models of these nebulae. No evidence of 
pulsation has been detected at these high energies ([Weekes et al. 19891) . 



Perhaps more intriguing is the detection by the Whipple observatory of Mrk 421 
( [Punch et al. 199^ ) and subsequently Mrk 501. These are both relatively nearby 
EGRET detected BL Lacs with hard spectra. Here again the spectra seem to match 
the extrapolation from the EGRET data. These detections suggest that the 7-ray 
production mechanism in blazars extends out to TeV energies. The non-detections of 
the stronger but more distant blazars such as 3C273 and 3C279 ( [Kerrick et al. 19"95D 
raises the possibility that these sources are extinguished at higher energies by 7 — 7 
interactions between these high energy 7-rays an the intergalactic IR radiation field 
( [Stecker 19921) . 



CHAPTER 7. THE GAMMA-RAY LARGE AREA SPACE TELESCOPE 



141 



Unfortunately, the cut-off in these sources occurs below the region of the spectrum 
accessible to ground based telescopes and above the the maximum energy to which 
EGRET is sensitive. This hole in the explored spectrum extends from ~ 10 GeV to 
~ 300 GeV. Not only is this an important region of the electromagnetic spectrum 
to explore for the reasons given above but the total lack of data in this band leaves 
open the door for unexpected new discoveries. 

7.2 Design Criteria 

The above considerations can be used to inform the design of a future 7-ray telescope. 
The most obvious requirement of such a telescope is that it possess a greater active 
area than EGRET. This is of crucial importance in improving the counting statistics 
which remain a limiting constraint on 7-ray astronomy. This in turn would increase 
the point source sensitivity as well as extending the accessible high energy spectral 
range, both of which are crucial performance criteria. Because of the inherent limita- 
tions on absolute detector size imposed by the difficulty of putting a large instrument 
in orbit, the effective area goal should be to design an instrument of geometrical area 
comparable to EGRET whose detection efficiency approaches unity. The maximum 
EGRET efficiency is of order 10%. 

A second generic requirement of a future 7-ray telescope is that the point source 
localization be improved. To some extent this will be a consequence of larger effective 
area and improved statistics, but it is more directly a function of the instrumental 
point spread function. A reduction in the PSF width would allow more accurate 
source positions which will in turn aid the search for counterparts in other spectral 
bands. Furthermore, improved resolution will allow a more detailed analysis of the 
structure of the galactic diffuse emission which will in turn allow more accurate models 
of this radiation to be constructed. In order to achieve an improvement in resolution 
for high energy 7-rays in a pair conversion telescope, the tracker element spacing must 
be reduced or the tracker depth increased or both. For low energy photons whose 
resolution is limited by multiple scattering it is necessary to attempt to minimize the 
thickness of the converter foils in order to preserve the directional information as the 
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electron/positron pair travels through the tracker planes. 

In addition to these rather straightforward telescope improvements, the EGRET 
data has added new constraints on future instruments. The unexpectedly high degree 
of time variability of many of the detected sources places a premium on the ability 
to monitor many sources continuously. This requirement argues strongly in favor of 
an instrument with an enlarged field of view in order to maximize the chances of 
observing flaring sources at all phases of their emission cycle. In a pair conversion 
telescope, this capability requires a flattened aspect ratio in order to permit oblique 
photons to travel through a tracker to a calorimeter. 

Furthermore, the intriguing detections and non-detections of EGRET sources at 
TeV energies emphasize the importance of extending the high energy reach of a future 
telescope. The self-veto problems along with the falling spectrum of sources at high 
energy combine to limit EGRET^s range to < 10 GeV. If this problem could be 
overcome in a larger telescope it should be possible to close the gap with the ground 
based telescopes at around 100 GeV. 

7.3 Gamma-ray Large Area Space Telescope (GLAST) 
Baseline Design 

7.3.1 Introduction 

In the more than twenty years since the design of EGRET there have been major 
technological advances that allow the basic principals of the EGRET design to be 
upgraded to obtain significantly increased performance. Most important among these 
are the improvements in silicon strip technology driven by the application of this 
technology to high energy physics detectors, and the vast improvements in available 
microprocessor power. A collaboration between scientists at Stanford University, the 
Stanford Linear Accelerator Center (SLAG), the Naval Reseach Laboratory (NRL), 
and other institutions ( [Michelson et al. 1995| ) has incorporated these technological 
advances into a prototype 7-ray telescope design. This instrument has been given the 
name: Gamma-ray Large Area Space Telescope (GLAST). 
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Figure 7.1: Schematic of the GLAST instrument showing one tower's internal struc- 
ture. 
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Figure 7A_ shows a schematic view of the proposed GLAST instrument as well 
as a similar view of EGRET to scale. The GLAST detector consists of a 7x7 set of 
towers. Each tower consists of three familiar subsystems. The surface of the tower 
is covered with an anti-coincidence scintillator layer. The readout for each tower's 
anticoincidence layer is independent of the others. Beneath this layer in the tracker. 
This tracker consists of 12 planes of silicon strip arrays. Each plane in the tracker has 
two sets of strip detectors that can accurately measure the charged particle tracks in 
two orthogonal dimensions. The first ten planes are covered thin layer of lead to serve 
as converter material to initiate the electromagnetic showers. Finally, the showers are 
collected in a calorimeter at the bottom of the tower. The design of this calorimeter 
is still in progress but for the purposes of this discussion it will be assumed to be an 
array of Csl crystals each 10 radiation lengths long. 

The schematic views of EGRET and GLAST illustrate several major advantages of 
GLAST over EGRET. First, the added size of GLAST results in increased effective 
area. Furthermore the aspect ratio of GLAST is much flatter than EGRET. This 
results in a greatly increased field of view which as discussed above, is of critical 
importance for monitoring time variable sources. This increase in field of view is 
made possible by the absence of a time of flight system that is used by EGRET to 
reject upward moving cosmic-ray showers. This innovation is made possible by the 
technological advances already mentioned and will be discussed in detail below. 



7.3.2 Veto Layer 

The EGRET veto layer functioned very well in rejecting cosmic ray backgrounds (see 
Chapter 3). The drawback in the design of the EGRET AG dome is the lack of 
segmentation. All fluorescent light in the AG dome is collected by a set of phototubes 
whose signals are coadded. As a result no positional information is available from the 
AG dome. This is problematic because high energy 7-ray events create large showers 
in the calorimeter. It is not unlikely that a charged particle from such a shower should 
travel back through the spark chamber and trigger the AG dome. The end result of 
this is that EGRET tends to self-veto on high energy events. 
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This effect, called the 'backsplash effect', has been minimized in the GLAST design 
by segmentation of the veto layer. Each individual tower retains its own layer of plastic 
scintillator. This plastic scintillator layer is read out with waveshifting optical fibers 
and Avalanche Photo Diodes (APD's). If an event originates in a particular tower 
whose veto layer received a hit then the event is rejected. However, if this veto hit 
occurs in any other tower, it is assumed to be backsplash and the event is processed. 

7.3.3 Tracker 

As mentioned previously, developments in the use of silicon strip technology for parti- 
cle detection provide the main technical stimulus for the GLAST tracker. These solid 
state devices replace the obsolete spark chamber technology used aboard EGRET. 
The principal advantages of silicon strip technology over spark chambers are: 

• Silicon strips can achieve ~ 10 times better position resolution. 

• They exhibit no appreciable dead-time. 

• There are no consumables that limit mission lifetime. 

• They do not require high voltage. 

• The absence of a chamber gas and a replenishment system gives improved safety 
and reliability. 

Taken together, these advantages amount to a very significant improvement in per- 
formance. 

The tower tracker planes are constructed from industry standard silicon wafers of 
6 cm X 6 cm with a thickness of 300 — 500 /im. A tower plane consists of two 3x3 chip 
arrays on top of each other; one layer each for the x and y position measurement. 
Double sided devices could theoretically be used but the increase in cost does not 
justify the added performance. Each chip has diode strips implanted on them to 
establish a readout pattern. The strip separation (pitch) is chosen as a compromise 
between resolution and channel number which is limited by power requirements as 
will be discussed below. The baseline design uses a strip pitch of 240 iim. This pitch 
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does not present a manufacturing challenge. Strip pitches of 50 /jm are routinely 
used in particle physics experiments. The chosen pitch will result in an rms position 
resolution of ~ 69 /xm in x and y. 

All but the final two tracker planes are covered with a 0.05 radiation length lead 
foil. The omission of converter material from the final two layers is a result of the 
requirement that at least three planes register hits in order to reconstruct a track. Any 
conversion that would take place in the final two layers would be unable to produce 
a valid trigger. The inclusion of converter material in these two layers simply serves 
to degrade the resolution of prior conversions through multiple scattering. 

One critical design feature of the tracker is that the converter material be very 
close to the first silicon plane. It is apparent from the examination of a shower 
development that this is important to preserve low energy position resolution. A 
conversion in a lead foil will immediately give a hit in the first silicon layer. This 
will give a very accurate measurement of the initial conversion point. After traveling 
through the second converter, the secondary pair will undergo multiple scattering. 
However, this scattering will be of minimal importance because the secondaries will 
have a high probability of interacting in the second silicon layer. Because this second 
layer is very close to the converter, there is a very small lever arm for the scattering 
to degrade the position. As a result the first two layers after conversion retain a large 
amount of directional information at low energies. 

At higher energies the directional resolution is no longer limited by multiple scat- 
tering but rather the geometric resolution resulting from finite pitch size. In the 
current design in which the planes arc separated by 3 cm the rms angular resolution 
approaches 0.03° in the limit of infinite energies. 

7.3.4 Calorimeter 

The photon energy is collected and read out in the calorimeter. The calorimeter must 
provide good energy resolution in the energy range from 20 MeV to 100 GeV. As in 
the case of the veto layer, pixelation is introduced in order to give a more complete 
image of a shower which will aid in background rejection. 
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The current GLAST baseline design calls for the use of thallium doped cesium 
iodide crystals. CsI(Tl) provides excellent energy resolution at modest cost, has a 
short radiation length (1.86 cm), provides a fast signal, and is radiation hard. It 
is also a rugged material that has been used in space previously. The calorimeter 
thickness is a compromise between weight constraints and high energy resolution. 
The current design achieves this compromise by having 10 radiation length crystals. 
Each tower site on top of a 64 element array of 3 cm x 3 cm crystals. 

The calorimeter described above can resolve 7-ray energies out to ~ 60 GeV. At 
energies higher than this the bulk of the shower leaks out the back of the calorimeter. 
Because GLAST sensitive areas are large enough to generate a significant number 
of such events, it would be desirable to be able to glean so information from these 
events. This is particularly important in view of the fact that one of the major 
goals of a future satellite borne instrument will be to bridge the gap with air shower 
experiments. 

It is possible to retain some energy information from these massive showers if 
the shower can be imaged longitudinally within the calorimeter. Several advanced 
calorimeter designs arc under consideration which explore the feasibility of this idea. 
These include longitudinal layering of the Csl crystals, and lead/scintillating fiber 
'spaghetti' calorimeters. 

7.3.5 Computing and Power Requirements 

A recurring theme in the design outlined above is the segmentation of many of 
EGRET^s monolithic systems. This segmentation increases the effective area by 
allowing more selective triggering. The price that must be paid for this additional 
area is a non-trivial on-board computing requirement as well as a restrictive power 
requirement. 

In order to improve the tracking resolution, GLAST has increased the channel 
number dramatically from the numbers used on EGRET . A quick channel number 
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indicates that 3 million channels will be included. If the strips are ganged together 
across an entire tower the number of channels is reduced to ~ 10^. Each of these 
channels requires power for a preamplifier and readout chip. Current state of the art 
techniques in low power electronics indicate that a strip of this capacitance can be 
read out using optimized electronics requiring only 350 /iW. This places the power 
requirement for the entire tracker at ~ 350 W. This represents a significant increase 
over EGRET which consumes only ~ 150 W for the entire instrument, however, it 
does not pose a serious power generation challenge. 

Another issue raised by this power requirement is the ability to dissipate the heat 
generated by this power. No high Z materials such as copper are permitted inside 
the tracker because they would cause scattering. Current schemes use aluminum to 
conduct away the heat. Berrylium is also under investigation. 

Perhaps the most fundamental departure of the GLAST design from that of 
EGRET is the trigger scheme. The EGRET instrument was designed to have a 
very selective trigger. The requirements of a valid TOF measurement, a threshold 
energy in the calorimeter, and the absence of a veto from the AC dome combine to 
limit the trigger rate to 0.7 Hz. This was done because the spark chamber gas de- 
grades with each trigger. Even with this selective trigger, the EGRET instrumental 
sensitivity had dropped to 40% of its pre-flight value before the latest spark chamber 
gas replenishment ( [Sreekumar 1996| ). Furthermore, the spark chamber requires a 100 
msec recharging time which implies that a much faster trigger rate would result in 
severe dead time complications. 

Because GLAST will not share these spark chamber limitations, the GLAST trig- 
ger can be much less selective. This is reflected at a design level by the complete 
absence of a TOF system. The philosophy of the GLAST design is to record as many 
events as possible while relying on pattern recognition in the hit patterns to reject 
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background. 

As mentioned previously, the cosmic-ray background in low earth orbit is typically 
10^ times greater than the corresponding 7-ray signal at a given energy threshold. In 
a particle physics environment, it would be possible to record both background and 
signal events which would later be analyzed in order to distinguish the hadronic events 
from the electromagnetic ones. In a satellite environment, however, the trigger rate 
is limited by the availability of ground communication bandwidth. If a conservative 
upper limit to the communication rate is adopted to be ~ 100 kbits then we can 
estimate the maximum allowable trigger rate based on the average event size. The 
~ 10^ channels are addressable with a 20-bit word. If an average event contains ~ 100 
hits, then 250 bytes are required just to describe the tracker pattern. This should 
almost be doubled in order to account for calorimeter data as well as other signals 
(e.g. current pointing knowledge, etc.) to give an average event length of ~ SOObytes. 
The maximum tolerable trigger rate is then ~ 25 Hz. 

The challenge in implementing this sort of trigger scheme is that it requires non- 
negligible on-board computing power. It is estimated that this trigger scheme could 
be implemented using a 15 MIP processor. The technological advances in micropro- 
cessor have resulted in space qualified chips capable of delivering this sort of comput- 
ing power. Radiation hard versions of the commercial RISC processors ( e.g. IBM 
RAD6000) will most likely be available for a GLAST launch. This is one of the main 
technology improvements that lead to the increased power of GLAST. 

7.4 Monte Carlo Simulations 

In order to gain insight into the performance of the GLAST instrument, the powerful 
technique of Monte Carlo simulation has been employed. Software developed for the 
simulation of particle physics detectors has been used to simulate both 7-ray and 
hadronic showers in a virtual GLAST and the results are then used to characterize 
the detector. The simulations described below have been carried out at SLAC by 
W.B.Atwood and at University of Washington by T.Burnett. 

The simulations described below use the state of the art simulation package named 
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Gismo ( Atwood et al. 1992| ). This program incorporates the QED interactions from 



the EGS4 code (pair production, bremsstrahlung, Compton scattering, MoUer scat- 
tering, Bhabba scattering, pair annihilation, etc.)) and hadronic interactions from the 
Gheisha code ( proton-nucleus interactions, vr-nucleus interactions, 7r° decay, charged 
Ks decay, neutral Ks decay, etc.). Showers are calculated down to an energy of 200 
KeV at which point the energy is assumed to be absorbed. To date more than 10^ 
photons showers have been simulated by Gismo using the GLAST baseline geometry. 



Figure \l.2\ and Figure show a comparison between a simulated electromagnetic 
shower caused by a 1 GeV 7-ray and a hadronic shower caused by a 15 GeV proton 
which enters GLAST from such an angle as to miss the veto layer. The particle tracks 
are shown in the first image and the resultant hit pattern in the detector is shown in 
the second frame. The comparatively tighter, more coUimated shower is produced by 
the 7-ray. This provides one of the principal background rejection techniques as will 
be described below. 

7.4.1 Background Rejection 

The GLAST trigger design is multi-level. The level 1 trigger is simply the requirement 
that the silicon tracker show hits in three consecutive planes. The Monte Carlo 
simulations described above using a realistic cosmic ray spectrum indicate that this 
level one trigger rate will be ~ 5 kHz. The level 2 trigger requires the best fit track 
to be reconstructed. A clean 7-ray track will have no hits prior the conversion point 
whereas a hadronic shower is likely to have scattered hits above the best fit track 
conversion point. The simulation results indicate that a cut on the distance of closest 
approach of a hit prior to the best fit conversion point to the best fit track, can 
give a background reduction ratio of 250:1 while only eliminating a few percent of 
the good 7-rays ( [Michelson et al. 1995 ). The second level trigger reduces the trigger 



rate to ~ 20 Hz which is sufficient to allow the data to be funneled through the 
communications bottleneck. 

Down-linked data can then be subject to further pattern recognition cuts in much 
the same manner as EGRET data. A series of cuts on the Monte Carlo data have 
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Figure 7.2: 1 GcV photon simulation event. The particle tracks are shown in white 
and the energy deposited in a given strip is shown as a vertical black line. The energy 
deposition pattern in the calorimeter is also shown as a histogram over calorimeter 
elements. The lower panel shows a close-up of the hit pattern alone. Notice the 
coUimation of the tracks. 
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Figure 7.3: 15 GeV proton simulation event. Again both particle tracks and resultant 
hit pattern are shown. Such a proton would miss the veto layer and thus would have 
to be rejected on the basis of the track pattern. Notice the track scatter compared 
to the electromagnetic shower produced by the 7-ray event. 
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successfully demonstrated the ability to reject cosmic rays at the level of 5 x 10~^ 
while preserving 80% of the 7-ray induced triggers. While this software can certainly 
be improved in order to give still further selectivity, the background rejection already 
demonstrated leads to an cosmic ray induced background an order of magnitude less 
than the extra-galactic diffuse radiation. 

It is worth highlighting the fact that such method of background rejection is al- 
ready in use on the EGRET data. The discussion in Chapter 3 of the EGRET 
instrumental background indicated that there is almost certainly a significant con- 
tribution of cosmic rays to the EGRET hardware trigger rate. Nevertheless, these 
events appear to be very efficiently rejected by the SAGE software package developed 
by the EGRET team. 



7.4.2 Performance Parameters 

All the tools necessary for characterizing the performance of the GLAST instrument 
are now in place. A point source of 7-rays is simulated using the Monte Carlo package. 
The resultant events are subject to the flight trigger and passed through the track 
recognition software. A virtual instrument calibration has been performed. The 



results are shown in Figure 7.4. 



7.5 Simulated Scientific Performance 

The techniques described above have provided a parametric description of the GLAST 
instrument. Based on these parameters, the capability of GLAST to answer some of 
the relevant scientific issues will be examined in this section. 



7.5.1 Simulation Package 

Several improvements in the GLAST performance lead to increased point source sen- 
sitivity. The increased sensitive area improves the statistics. The extended energy 
range allows one to make full use of the high energy photons which carry the most 
positional information, the sharper point spread function increases the signal to noise 
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Angular Resolution vs Energy 
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Figure 7.4: GLAST calibration as derived from the Gismo simulations. Also shown 
for comparison are the EGRET parameters. Note in particular the dramatically 
improved sensitive area at high energy as well as the greatly expanded field of view. 



ratio, and the wider field of view increases the average coverage of a source which in 
turn leads to improved statistics. The most direct way of quantifying the degree of 
improvement is again through Monte Carlo simulations. 

The first step in the construction of a simulation package is the calculation of the 
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appropriate broad energy PSF for GLAST. 

1 



PSF{e) = — dE PSF{9,E) / E'-'^dE' EDF{E',E) SAR{E') , (7.2) 

where PSF{6), EDF{E' , E), SAR{E') are the point spread function, energy disper- 
sion function, and sensitive area function respectively. In general, the spectral index 
a is chosen to be 2.0. 



Figure [7. 5| shows the GLAST PSF calculated using the parameters taken from the 
Monte Carlo results. It should be noted that the increased importance of the high 
energy photons results in the greatly enhanced 'spike' at the center of the PSF. 

This PSF is then used to generate the GLAST background model. For this the 
galactic diffuse model described in Chapter 4 was convolved with the GLAST PSF. 
The results are shown below compared with the EGRET background model. The 
improvement in resolution is readily visible in this comparison. This added resolution 
will provide powerful new data with which to constrain such models. The construction 
of more accurate galactic models will be a crucial part of the future improvements in 
measurements of the IDGRB. 

The exposure map can also be generated using the parameters that have been 
established. The following assumptions are made about the viewing mode. The 
telescope operates in a zenith pointed mode in which the pointing direction is swept 
across the sky in such a way as to always face away from the earth. In this way the 
total exposure to the sky is maximized because dead time spent looking toward the 
earth is eliminated. Furthermore, assuming a long enough integration, this exposure 
should be equally spread across the sky. The total exposure to a source of spectral 
index a is then, 

^ _ Jo^{Ao cos 9)271 sin9d9 _ AqTu^ 

where 

J^^-- SAR(E)E-^dE , ^ 

= — , 7.4 

§^-^E-^dE ^ ' 

is the spectrally weighted on-axis sensitive area. Using the sensitive area results 
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Figure 7.5: The simulated GLAST PSF compared to the measured EGRET PSF. 
The top panel shows the probability of detecting a photon per unit steradian whereas 
the bottom plot shows the probability per degree. The central spike in the GLAST 
PSF is very pronounced. 
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Figure 7.6: Comparison between tlie diffuse galactic emission with E > 100 MeV as 
seen by GLAST (above) and EGRET (below). 

from above we find that for E > 100 MeV, = 7300cm^ and consequently after 
a 1 year full sky survey the average exposure to any given point on the sky will be 
£ — 5.78 X lO^'^cm^s. This can be compared with the average exposure obtained by 
EGRET Regret = 4.3 x lO^cm^s. This amounts to a factor of 130 improvement over 
EGRET. If we consider where this improvement arises we find that there is about a 
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factor 10 improvement in Aq, about a factor 4 improvement in field of view, and an 
additional factor 3 in live time improvement. This last factor comes partly from the 
fact that GLAST can be operated in a zenith pointed mode because of its wide field 
of view whereas EGRET must spend quite a bit of time looking at the earth. An 
additional consideration is that the failure of the flight tape recorders increased the 
dead time by a factor of two above the nominal value. 

In order to analyse point sources in simulated GLAST data the EGRET maximum 
likelihood analysis package was adapted to accept GLAST parameters. The only 
substantive difference between this package and the EGRET package was that in 
the case of GLAST a large amount of positional information is contained in the 
E > 1 GeV photons and it is important to preserve this information by not binning 
the data together in wide energy ranges as is done for EGRET. Instead analysis is 
done separately in the energy ranges 100 < E < 1000 MeV and E > 1 GeV and the 
two independent results are used in order to ascertain significances and point source 
locations. Because the two data sets are independent the likelihoods are multiplicative 
and thus their logarithms are additive. As a result the likelihood values in the two 
energy ranges can be summed to yield a joint likelihood measurement. 



7.5.2 Point Source Sensitivity 

The most critical measure of telescope performance is its point source sensitivity. This 
is best expressed as a relation between point source strength and the significance of 
its detection. As shown in Chapter 5, in the regime in which background counts 
dominate source counts, 

a = , ^ Ms , (7.5) 

/£iS + B) V B 



the significance, a, grows linearly with source flux S. Figure |7.7| shows the relation 
between average a and point source strength for two background conditions varying 
from high latitude to near the galactic center. Fifty simulations were analyzed in 
the two energy bands for each point in order to ascertain the mean source detection 
confidence. 
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Figure 7.7: Point source significance versus point source flux. The relationship is 
shown for typical points on the galactic plane and at high galactic latitude. 



The results indicate that the high latitude point source 5a detection threshold is 
~ 2.6 X 10^^ whereas near the galactic center it is ~ 8 x 10^^. This amounts to a 
factor 40 improvement over EGRET performance. 



7.5.3 Point Source Location Accuracy 

Correlated multiwavelength observations provide a powerful tool in the investigation 
of 7-ray sources. The applicability of this method is severely limited by the relatively 
large size of the EGRET error boxes. A typical 95% EGRET error box is ~ .5° 
in diameter ([Thompson et al. 1995| ). The number of potential counterparts to these 
7-ray sources that lie within these error boxes can be quite large. A query of the 
NASA Extragalactic Database (NED) yields an average of ~ 10 extra-galactic radio 
sources within the typical high latitude error box. A critical performance criterion 
for GLAST is the improvement in this area. 

Distributions of the 95% error box size were determined using the likelihood maps 
in the neighborhood of simulated sources. Likelihood maps in each independent 
energy range were constructed and summed to yield a joint likelihood map from which 
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the likelihood contours were measured. The relationship between error box size and 
source strength was explored for sources at high latitudes and on the galactic plane. 
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Figure 7.8: Scaling of the 95% error box radius with source strength. The crosses 
indicate the error box radius for the EGRET catalog sources and the dotted line ins 
the fit to S*"^/^ at the high fiuxes. The solid line represents the results of a simulation 
of the GLAST performance at high latitude. The dashed line shows the performance 
near the galactic center. 

The interpretation of the results trend is straight forward. In the limit of no 
background, the source location is determined by the width of the PSF, 9rms: and the 
source strength, S, which determines the total number of detected photons from the 
source, N. 

^ = ^ a OrmsS"^ . (7.6) 

As a result the confidence region should scale as S~^^'^ whenever the source dominates 
the background. This is clearly seen in the data for sources with source strengths 
greater than ~ 10~^. Weaker sources on the plane start to show the effects of back- 
ground fiuctuations and their error boxes are correspondingly enlarged. At high 
latitudes this effect begins at a lower fiux. 
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7.5.4 Point Source Confusion 

All the above analyses have been completed using the most ideahzed case of an isolated 
point source on a perfectly modelled background. On the actual sky this is never the 
case. Assuming a simple Euclidean power law extrapolation to low source fluxes 
from the EGRET data, we find that there will be ~ 8660 extra galactic sources 
with fluxes greater than the 5a flux threshold at high latitudes of 3 x 10~^. This 
implies that the mean separation between sources is ~ 1.2°. A certain fraction of 
these sources will inevitably be obscured by brighter objects. Furthermore, sources 
beneath the detection threshold will produce fluctuations which generate local areas 
of discrepancy between the diffuse model and the data. This will tend to obscure 
otherwise resolvable sources. 

This issue has been approached in two ways. The simplest way to estimate the 
number of confused sources is to use the above calculation of GLAST error box 
sizes to estimate the number of sources which would overlap other sources. For the 
purposes of this calculation a source is declared to be obscured if another source 
has an overlapping error box and is of equal or greater strength. If a logN — logS 
relationship is then assumed for the extra-galactic sources, an estimate can be made 
of the total number of obscured sources. 

Source A is confused if the separation between it and a second source {9) is less 



where Sa and Sb are the source fluxes of source A and B respectively. The average 
number of sources that confuse source A is given by. 



than 9, 



9., 



max 



^95('S'a) + ^95('S'b) , 



(7.7) 





and thus the probability that a given source will be confused is. 



P 



1-e 



NiSA) 



(7.9) 



In the simplest case where we assume that the 95% error box radius scales as S 
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and that there is a Euchdean logN — logS source distribution, this probabihty can 
be calculated analytically. 

N{Sa) = f^NM^r^/' , (7.10) 

where A^o represents the number of sources in the sky above the reference flux Sq and 
do is the 95% error box of a source with that strength. Using the best fit Euclidean 
source distribution to the EGRET data and the error box size determined above the 
number of confused sources at the detection threshold of 3 x 10~^ is ~ .03. Thus the 
probability of a source at the detection threshold being confused is ~ 3%. 

A second approach to this issue is to simply simulate an observation of a portion of 
the sky including a sample logN — logS. A Poisson sample is then simulated and the 
standard likelihood techniques are then used to reconstruct the sources in the field. 
The reconstructed logN — logS can then be compared to the input distribution in 
order to ascertain the effects of obscuration by brighter sources as well as the effects 
of unresolved sources creating local inaccuracies in the background model. 

Figure 7|9 shows the comparison between the input logN — logS curve and the 
corresponding reconstructed curve. The detection threshold which is determined by 
the smallest detected source strength is ~ 3 x 10^^. This is slightly higher than 
the corresponding isolated source detection threshold because the fluctuations in the 
unresolved sources tend to lower the likelihood values. Source confusion in this rep- 
resentation would appear as an inconsistency between the two distributions. There 
appears to be little confusion until S* < 5 x 10^^ at which point the reconstructed 
source distribution deviates toward the slightly higher number of sources. While this 
may at first seem counterintuitive, what is occurring is that unresolved sources in 
the vicinity of detected sources are contributing to the detected source's intensity. 
The result is to push that source to higher flux leading to an apparent over detec- 
tion of sources at fluxes slightly higher than the detection threshold. For a detailed 
discussion of this effect see QSchmitt 1986|) . 

The detection accuracy can be examined in more detail by looking at the spatial 
maps of the simulated region. Figure |7.10| shows the intensity map of a small region 
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Figure 7.9: The results of a simulation of a GLAST observation of the high latitude 
sky. The gray line shows the input distribution of sources. This distribution was 
chosen at random within a 16° x 16° field. The black hue shows the reconstructed 
source distribution after analyzing this field using the maximum likelihood techniques. 
Sources are detected down to a threshold of ~ 3 x 10~^. Near the detection threshold 
the source fluxes are biased toward slightly higher fluxes due to the contributions of 
unresolved sources. 




Longitude 



Figure 7.10: Simulated sky intensity for a high latitude region. The map shows the 
modelled intensities of point sources in this region convolved with the GLAST PSF. 
Also shown as cross hairs are the point source positions reconstructed using maximum 
likelihood from a Poisson sample of this map. The correlation in positions shows that 
the reconstructed sources do indeed represent true sources. In some cases multiple 
sources are shown to be reconstructed as a single source. 
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in detail. This map does not have Poisson statistical fluctuations added so as to 
present the actual model used in the simulation. The reconstructed source positions 
are shown as crosses. One can see that in some cases, a strong source has been 
grouped together with some weaker sources leading to the detection of a source with 
a slightly increased flux. 

One way to quantify this effect is to integrate the additional flux which leaks in 
in the form of source confusion. In the case of the 16° x 16° region simulated, the 
total integrated flux above S = 3 x 10^^ was 3.0 x 10^'' for the input sources and 
3.3 X 10^^ for the reconstructed sources. As a result the total flux source confusion 
flux is ~ 3 X 10^^ which amounts to ~ 10% of the total flux in the sky. This number 
is of course dependent on the exact form of the logN — logS relation chosen. For 
the purposes of this simulation the steepest logN — logS allowed by the fluctuation 
analysis of the EGRET data was used so as to get as conservative an estimate of the 
source confusion problem as possible. 



7.5.5 Spectral Performance 

The ability of EGRET to explore the highest energy 7-rays was limited by two factors. 
Firstly, its relatively small area limited the number of high energy photons incident 
on the telescope. Furthermore, those high energy photons that were incident on 
the telescope had a high probability of being self-vetoed which further limited the 
effective area. In GLAST there has been an attempt to overcome these problems 
with increased detector size and by elimination of the backsplash problem through 
modularization of the veto layer. 



The sensitive area of GLAST as a function of energy is shown in Figure |7.4 In 
contrast to the case of EGRET, the sensitive area reaches a maximum around 1 GeV 
and remains at that level as energy increases. What ultimately limits the sensitivity 



of GLAST is the falling spectra of the celestial sources of 7-rays. Figure [7.1 1| demon- 
strates the capabilities of GLAST in the energy range 1 GeV < E < 1 TeV. The 
brightest extra-galactic source in the EGRET sky is 3C279, a quasar that underwent 
a large flare during the first year of EGRET observations ( [Knitfen et al. 199^ ). This 
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blazar has a spectral index of ~ 2.0 and would be clearly visible to ground based 
air shower telescopes were its spectrum to continue unbroken to TeV energies. Its 
non detection leads to the conclusion that either this source has an intrinsic spectral 
break in the region between 5 GeV and 500 GeV or that there is sufficient extinction 
of the high energy 7-rays through 7 — 7 interactions with the interstellar infrared 
background to extinguish the source at high energies. 



10000.0 = 
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Figure 7.11: Measured spectrum of 3C279 after a one year GLAST exposure. The 
dashed line shows the extrapolated spectrum from the EGRET data. The solid line 
assumes no intrinsic source cutoff but allows for extinction by the interstellar infrared 
radiation field as calculated by Stecker et al. (1993). The lower graph shows the 
corresponding EGRET spectra. The GLAST spectra start to be limited by counting 
statistics at energies greater than ~ lOOGeV. It should nonetheless be possible to 
detect the spectral break which must occur at energies less than 1 TeV. 



Figure [7.11| shows the spectrum of 3C279 that would be produced by GLAST after 



a one year all-sky survey. The source is assumed to have undergone a flare similar 
to that observed by EGRET during the course of this year in order to produce an 
integrated time averag ed flux above 100 MeV of 1.0 x 10~^ A spectral index of 2.0 
is assumed consistent with the spectral index measured by EGRET. Two spectral 
models have been convolved with the GLAST response. The ffist (shown by a dashed 
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line) represents a simple extrapolation of the EGRET spectrum to TeV energies. This 
spectrum is inconsistent with the Whipple measurements at TeV energies and is shown 
for reference purposes. The second spectral model assumes that there is no intrinsic 
source cut-off but estimates the infrared extinction as calculated by Stecker (|Stecker 
et al. 199^ ). In order to demonstrate the statistics of such a spectrum directly, the 
values plotted are the number of photons in each spectral bin. It is worth pointing 
out that while the bin centered on 250 GeV contains only 2 photons, the chance that 
2 background photons in this energy range arrive within the 140 arc second error box 
of these high energy photons is 5 x 10"'^. The relative absence of background photons 
means that significant detections can be made with comparatively small number of 
photons if one has an a priori source position. Also shown for comparison is the same 
model convolved with the EGRET response. Notice that the EGRET self-vetoes on 
approximately half the converted photons at 8 GeV. 

The GLAST spectrum of 3C279 demonstrates the ability of GLAST to produce 
good quality spectra out to energies of ~ 100 GeV for the brightest EGRET sources. 
Beyond this the spectra will be limited by statistics but this instrument should 
nonetheless be capable of resolving the spectral break that must occur below 1 TeV. 
Of course an intrinsic spectral break in the source at lower energies would be easily 
resolvable. Detections of such breaks in the brightest 10 sources would provide crucial 
information as to the nature of the acceleration mechanism in blazars as well as hav- 
ing the potential to address fundamental issues of cosmology by gleaning information 
about the infrared extinction of distant sources. 

7.5.6 Diffuse Spectrum 

The extended energy reach of GLAST will allow the IDGRB to be probed out to 
higher energies. Furthermore, the detection of several thousand blazars will allow the 
blazar spectrum to be accurately measured. These two measurements will allow any 
remaining IDGRB to be probed for distinctive spectral features that might signal the 
presence of a non-AGN component of the IDGRB. 

In Chapter 1, the possibility of WIMP decays or annihilations contributing to the 
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IDGRB was discussed. The sensitivity of GLAST to the spectral signals predicted by 
these models can now be assessed. Figure [7.12| shows IDGRB spectrum as observed 
by GLAST (upper set of curves) and EGRET (lower set of curves). We have assumed 
that a large fraction of the sky will be unsuitable for the measurement of the diffuse 
background because of the presence of strong point sources of strong foreground 
galactic emission. As a result the data shown in the figure represents only 10% of 
the sky. Two assumptions are made in order to produce the various spectra shown. 
The dashed spectrum corresponds to an AGN-like spectrum scaled to account for 
~ 50% of the IDGRB. The remaining 50% is assumed to be found in resolved AGNs. 
The dotted line represents the predicted observed flux from a WIMP decay in the 
early universe. Following the model described in chapter 1, this curve is calculated 
for a WIMP which decays at some characteristic redshift zd with mass mx and 
whose density is given by VLxh"^- These parameters determine the spectrum shown 
using equation (1.3). The simulation shown in figure x makes the following choice of 
parameters, 

E, = =10GeV (7.11) 

2(1 + Zd) 

nxh^ = 10-1^ 



The solid curve in Figure |7.12| shows the sum of the AGN component and the 
WIMP decay component of the IDGRB. For this choice of parameters a bump at 
~ lOGeV would be clearly visible. Because, we have a shown that it is possible that 
~ 50% of the IDGRB will be resolved as point sources by GLAST, the dashed curve 
also represents the approximate summed spectrum over all detected AGNs. In this 
way the IDGRB and the AGN component can be measured to the accuracy shown in 
the figure, allowing the WIMP decay component to be reconstructed if present. 



The corresponding EGRET observation is also shown in Figure [7.12| . EGRET is 



not sensitive to such spectral bumps. This is due to its smaller collecting area and 
the reduced effective area caused by 'backsplash'. 

Another spectral signal postulated in chapter 1 was the appearance of a 7-ray line 
caused by the annihilation of two WIMP's into 7-rays. The sensitivity of GLAST 
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Figure 7.12: Diffuse extra-galactic spectra as measured by GLAST upper curves 
and EGRET (lower curves). The dashed lines indicate the IDGRB produced by 
unresolved AGNs. The dotted line shows a possible WIMP decay spectrum. The 
solid curve shows the summed spectrum. GLAST is sensitive to this WIMP decay 
spectrum whereas EGRET is not. 

and EGRET to these lines has been calculated in the following way. The 7-ray line 
is assumed to have negligible width compared to the instrumental energy resolution 
AE/Eq. The signal counts from a 7-ray line of integrated intensity I are then given 

rEo+A 

Ns = IS EDP{E,Eo)dE, (7.12) 

JEo-A 

where EDP{E, Eq) is the instrumental energy dispersion function and £ is the ex- 
posure. If we assume this line is measured against a power law background then the 
background in this energy range is 

rEo+A 

Nb = Bo£ / E-^dE , (7.13) 

JEo-A 



where we have ignored the higher order correction caused by the convolution of 
the background spectrum with the instrumental dispersion functeon. If we assume 
EDP{E, Eq) to be a Gaussian with the relative width, erg, shown in Figure |7^, the 
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signal to noise is maximized near A = cte- We can then evaluate the significance of 
a line detection as simply, 

aune = . (7.14) 



VNs + N, 
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Figure 7.13: 7-Ray line sensitivities of GLAST and EGRET. The dashed (solid) 
curve shows the 7-ray line intensity that would result in a 2a bump on the IDGRB 
when observed by EGRET (GLAST). The dot-dashed line indicates the estimated 
line intensity due to the annihilation of the lightest super-symmetric particle in the 
galactic halo. Although cross-section used is highly speculative, it is likely that the 
detection of such WIMP annihilation lines will remain beyond the scope of GLAST. 



Figure |7.13| shows curves of constant significance as a function of Eq = mx', 
the WIMP mass. The dashed (solid) curve is shows the 2a sensitivity of EGRET 



(GLAST). Also shown are the predictions by Kamionkowski ([Kamionkowski 19951) of 
the expected line flux from neutralino annihilation as a function of neutralino mass. 
These predictions indicate that the sensitivity of GLAST will remain insufficient for 
the detection of these 7-ray lines. However, the cross section for the coupling of 
WIMP's to 7-rays are highly speculative and contain large uncertainties. GLAST 
could provide interesting constraints on these processes. 



CHAPTER 7. THE GAMMA-RAY LARGE AREA SPACE TELESCOPE 



170 



7.6 Summary 

EGRET has provided a tantalizing look at some of the most violent processes in the 
universe while at the same time raised a number of questions which remain beyond 
its scope. The motivation for probing the 7-ray universe in more detail is compelling. 
The GLAST instrument concept described in this chapter represents an effort to 
develop a realistic instrument for the future of 7-ray astronomy. The simulations 
presented here indicate that GLAST will become a powerful tool in the quest find 
answers for many of the questions raised by EGRET. 

Table 7.1: Comparison of GLAST and EGRET Performance 
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Table 7.1 summarizes the improvements in point source detection achievable by 
GLAST over EGRET. The dramatic increase in the expected number of detectable 
sources represents a significant leap for the field of 7-ray astronomy. Up to 40% of 
the remaining isotropically distributed emission will be directly resolvable as point 
sources. A further 25% will be constrained by fluctuation analysis. In addition, a 
population of several thousand blazars will be available. If redshifts for a large fraction 
can be obtained, the luminosity function should be much better constrained. As 
mentioned in chapter 5, the luminosity of the X-ray Seyferts has now been sufficiently 
well constrained with a comparable number of sources as to be able to identify these 
sources as the overwhelming source of the X-ray background. 

In addition, the spectrum of the isotropic diffuse radiation should be very well 
measured, any discrepancy between the blazar spectrum and the diffuse spectrum 
could provide further evidence of a true diffuse component to the extra-galactic diffuse 
7-ray background. The ability of GLAST to detect spectral features of the sort that 
would be produced by WIMP decay has been established. 
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The intriguing question of the nature of the unidentified EGRET sources would 
be greatly elucidated by an instrument such as GLAST. Between the decrease in error 
box size which would reduce the number of counterparts in other wavelength bands, 
the ability to do long term monitoring of these sources which is afforded by the vast 
field of view of the telescope would allow multiwavelength campaigns to be easily 
mounted. Such observations have proven to be critical to unraveling the mysteries of 
unidentified sources (e.g. Geminga). 

The abihty of GLAST to extend the energy range of sateUite borne 7-ray telescopes 
has also been demonstrated. This development is particularly desirable in hght of 
the fact that the ground based 7-ray telescopes are pushing to extend their effective 
spectral range to lower energies. The ability to cross calibrate the two technologies will 
greatly enhance the usefulness of both telescopes. Furthermore, the complementarity 
of a large FOV telescope with comparatively smaller effective area (GLAST) and a 
pointed telescope with larger effective area (future air shower telescope) should prove 
to be a very effective combination. 

Perhaps the single greatest possibility of an instrument such as GLAST is of 
course the truly unexpected and unpredictable discoveries that it will undoubtedly 
make. In addition to opening a new spectral window onto the cosmos, the two decade 
increase in sensitivity will undoubtedly yield unexpected and exciting sources of 7- 
ray emission. As an example, AGN were detected by EGRET in numbers that were 
totally unpredicted prior to launch. It is the possibility of finding such staggering 
phenomena in the future which should always be the primary motivation for pushing 
outward the frontiers of the detectable universe. 



Appendix A 



Maximum Likelihood Analysis 

A.l Maximum Likelihood Analysis 

Because of low gamma ray counting rates as well as the relatively intense and highly 
structured background, EGRET data analysis requires the use of statistical tech- 
niques. An overview of the techniques used in the analysis of EGRET data is pro- 
vided in (Fichtel et al. 1994). The analysis undertaken by the EGRET team thus 
far includes time resolved studies of gamma ray pulsars (Nolan et al. 1993), diffuse 
studies (Bertsch et al. 1993), and point source studies (Mattox et al. 1994). The 
principal tool used in the analysis of EGRET point sources is that of maximum likeli- 
hood. The EGRET implementation of the maximum likelihood method is described 
extensively in Mattox et al. 1996 and will only be summarized here. 

Likelihood, introduced by Fischer (1925), expresses the joint probability that a 
data set are drawn from a hypothetical distribution. Likelihood maximization is a 
useful technique for parameter estimation (Cash 1979) and has been used extensively 
in gamma ray astrophysics (Pollock et al. 1981). 

In a case like that of photon counting, where noise in a given pixel of solid angle 
is assumed to be Poissonian, likelihood for a set of pixels ij with measured counts riij 
is defined as: 
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where 6ij = 6{ak)ij are the predictions of a model with parameters a^. In the context 
of a given model, the most likely value of a parameter is that value which maximizes 
the value of Lq. In practice InLo is more conveniently calculated: 

In Lo = XI ^ii 111 % - X! % - X! In ^ij' > (A.2) 

ij ij 

This last term is model independent and can be dropped. 

The theory of maximum likelihood fitting also predicts that in the limit of large N 
(where N is the number of relevant photons) , parameters will be normally distributed 
about the maximum value with a o" of 



<^--{H-^)K (A.3) 

where a is the parameter under estimation and H~p is the inverse of the Hessian 
matrix defined as, 

This allows one to estimate not only expectation values but also uncertainties. 

The likelihood ratio test (Neyman and Pearson 1928) extended the use of like- 
lihood to hypothesis testing and significance evaluation. This statistic is used for 
evaluating the extent to which one model better describes the data than a second 
such model. It is defined as: 

T5 = -21ni^^, (A.5) 

where Lj is the likelihood for a model with m free parameters and Lq is the likelihood 
for a model with p < m free parameters. 

The advantage of this test is that its distribution in the null hypothesis (model 
is correct) is known analytically. Specifically this quantity was shown by Wilks 



APPENDIX A. MAXIMUM LIKELIHOOD ANALYSIS 



174 



(1938) to be distributed as Xn where n = m — p is the difference in the number of 
free parameters between the two models. 

In the case of EGRET data, hkehhood is used to estimate the flux of a measured 
point source. The models used contain three free parameters: two parameters describe 
the diffuse emission described below and one describes the strength of a putative point 
source. 

Oij = gmGij + QbEij + CaPSF{a, 5)ij , (A.6) 

where Gij describes the structured galactic background, Eij the isotropic background, 
and PSF{a, 6)ij the point spread function for a point source at location a, 6. The 
free parameters in this model are the scaling factors gm,gb,Ca- This model is used to 
find the best fit points source fiux (c^) and is compared to the model with no point 
source present (ca = 0) in order to determine TS and thus the significance. 

A. 2 Monte Carlo Simulations 

At the basis of this investigation is the creation of simulated data distributed ac- 
cording to a model of sources and background convolved with the EGRET telescope 
parameters. This data can then be processed using EGRET data analysis techniques 
to test various responses. The creation of such simulated data requires two inputs: a 
spatially resolved background model and EGRET instrument parameters. 

A. 2.1 Instrumental Parameters 

The EGRET instrument is characterized for the purposes of this study by the point 
spread function and the exposure map. Comprehensive cahbration of the EGRET 
instrument at SLAG indicate a Gaussian point spread function with energy depen- 
dent width(P5'F(^, E)). All analysis described in this paper was done using photons 
with £■ > 100 MeV with an averaged point spread function generated by weighting 
monoenergetic point spread functions with an assumed power law spectrum of the 
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source and taking into account the energy dispersion of the EGRET instrument. 

1 rEmax roo 

PSF{9) = — dE PSF{9, E) / E'-'^dE' EDF{E', E) SAR{E') , (A.7) 

where a is the spectral index of the source, EDFE', E) is the energy dependent 
dispersion, SAR{E') is the sensitive area. The integration over E' represents an 
integral over true energies and the integration over E represents an integral over 
measured energies and N is the appropriate normalization. In most spectral 
index of 2.0 is used. The background model is convolved with this PSF in order 
to obtain a model of the observed background. All added point sources are also 
convolved with this PSF. 



600 F 




Figure A.l: The EGRET PSF for E > 100 MeV. 

The exposure maps are generated in order to transform flux maps (e.g. galactic 
flux map) into predicted counts maps. Sensitive area (SAR) as a function of photon 
energy, declination angle , and operational mode was cahbrated at SLAG (Thompson 
et al. 1992). Live time is recorded on board and these two quantities can be used to 
generate maps of exposure. 

Using the diffuse model with any arbitrary set of added point sources and exposure 
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maps, one can calculate the rate of arrival of photons into any pixel of an observation. 
These values form the smooth "infinite statistics" map. The simulated observations 
are then generated from this smooth map by adding Poisson distributed deviates 
to each pixel. The rejection method (Press et. al. 1992) is used for this purpose. 
Output counts maps together with the exposure maps can then be analyzed exactly 
as flight data using the likelihood software developed by the EGRET team ( i.e. LIKE 
(Mattox et al 1996)). This program optimizes the likelihood with respect to the three 
free parameters: gm the diffuse galactic intensity; gb the diffuse isotropic intensity; Ca 
the point source flux. 



A. 3 Point Source Analysis 

Point source parameter fitting was investigated using statistical ensembles of observa- 
tions containing test sources. Sources of varying flux were placed in different points 
on the sky in different background conditions. Both on and off axis sources were 
studied and the effect of lengthening the exposure was investigated. 

Distributions of fitted flux, positions, and significance were generated. Because 
likelihood analysis makes predictions not only about the value of these parameters 
but also their uncertainties, it was possible to check the first two moments of the 
distributions obtained against the likelihood results. 



A. 3.1 Flux Estimation 



Figure |A.2| shows the results obtained for sets of 1000 observations of the same source 
as one varies the flux. The median value of the fitted flux agrees well with the true 
value. The effect of disallowing negative fluxes shows up in these distributions as a 
piling up of the distribution at zero flux. 

The flux distributions are well approximated by Gaussian functions. The distri- 
bution widths were calculated and compared to the results obtained by averaging 
the individual likelihood estimates of the flux uncertainty. The agreement is good 



Figure |A.2| ). One can see the transition between the background limited flux region 
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population 1 best fH parameters: _ 
Peak:1G2 +/- 3 
Sigrria^62 +/- 5 

F'apuldtion II beat 'fit parameters; 
Peak'41 3.0+/- 7 _ 
Sigma^bS +/— 2 
Like estimate of 9igma:64 +/— 1 



ma 



Figure A. 2: Two sets of simulations siiowing tlie distribution of measured counts given 
a simulated point source flux. The distributions are peaked about the correct flux 
and the uncertainties in these measurements is consistent with the mean uncertainty 
derived from the likelihood function. 

(where the uncertainty is dominated by fluctuations in the background counts) to the 
source limited flux region (where the uncertainty is dominated by the fluctuations in 
the source counts). 

A. 3. 2 Upper Limits 

The most useful form of a point source upper limit is deflned as follows: it is the 
point source flux above which higher true source fluxes can be ruled out with specifled 
confldence based on the actual measurement of the flux. In other words one would 
like to construct the distribution of true point source flux given a measured flux 
{P{S \ m)). Such a distribution would typically be peaked about the measured flux 
and would have some characteristic width which would determine the na upper limit. 
Such distributions are not readily obtained from the data, however, and what one 
most often calculates is the point source flux that would result in a measured flux 
greater than or equal to the actual measured value with some specifled confldence. 
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This latter upper limit is calculated in the following manner. The measured flux 
(m)is obtained by maximizing the likelihood. In the limit of large counts this flux 
is normally distributed with uncertainty obtained from the equation (3). If this 
uncertainty is not strongly dependent on the point source strength, the upper limit 
of the source flux is then S = m + na. That is one will have found the point source 
flux {S) 95% of whose measured flux distribution lies above the actual measured flux 
if one chooses n = 2. 

Monte Carlo simulations allows one to determine P(m | S), the probability distri- 
bution of m given S, and thus to test the various hypotheses used above to calculate 
upper limits. The fact that measured flux is normally distributed has already been 
demonstrated (Figure |A.2[). Furthermore the likelihood estimate of the flux uncer- 



tainty has been confirmed within errors. The final assumption to be tested is whether 
the flux distribution changes as one increases the point source flux. As is shown in fig 
X. the width of the measured flux distribution remains roughly fixed in the background 
dominated region. Because upper limits are only calculated for the low significance 
sources, this assumption remains valid. 

With a Bayesian approach, one can use the knowledge of P{m \ S) to calculate 
{P{S \ m)). Using Bayes theorem: 

P(S I m)ds = Pirn I S) * P(S)ds/ / P(m I S) * P(S)ds , (A.8) 

Jo 

where P{S) is the prior distribution of S. The problem then becomes what to 
use as a prior distribution of S. If one uses a uniform prior P{S) = k, one finds 
P{S I m) = P(m I S) and the two definitions of an upper limit become equivalent. 
If, however, one uses a prior based on a more realistic distribution of sources(e.g. one 
estimated from a measured logA^(logS') relationship), one can significantly reduce 
the upper limits. However, it is probably better in most cases to assume the most 
conservative prior which in this case is that of uniformity. 
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A. 3. 3 Position Estimation 

The method used to evaluate the position contours for a given source is outhned 
in Mattox et al. (1994). If one tests the hypothesis that there is a source of some 
strength at position Iq, Bq against the null hypothesis that there is no source present, 
one obtains a likelihood result. One can now make this hypothesis the new null 
hypothesis and evaluate the likelihood of a model containing a source of the same 
strength located at a position This new likelihood ratio can be optimized to 

determine the best fit position lmax,bmax- 

One can determine position confidence intervals in much the same way as for flux. 
If the source were in fact located at a position Iq, Bq (the null hypothesis), the position 
likelihood ratio would be distributed as (two additional degrees of freedom li,bi). 
We can define : 

T Spositioni^^Oy bo) lflL(^l>iYiax : bj^iax )-lnL{lo,bo) , (A.9) 

Thus the locus of points /q, &o giving a TSposiUon value of x determines the y% confi- 
dence contour for containing the source. Where x and y are related by the formula: 

y= rxliTS)dTS, (A.IO) 
Jo 

The distribution of this test statistic in the null hypothesis was tested using Monte 
Carlo techniques. Maps of the position test statistic were generated for sample ob- 
servations of a source, and the values of TSposiUonikrue, hrue) were recorded. The 
resulting distribution was a good fit to xl (Figure [A. 3] ). The slight underprediction 
near TSposUion = results from the finite bin size used in calculating the likelihood 
maps. Differences in TSposUion less than some threshold were likely to fall in the same 
pixel and result be assigned a T Sposuion of zero. 

A. 3. 4 Inclination angle effects 

Analysis of off axis sources have the added complication of map edges and poor 
exposure. Source distributions for off axis observations were generated in order to 
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Figure A. 3: Cumulative distribution of the TS drop between the best fit position 
and the true simulated point source position. The distribution is consistent with xl 
except near zero where the finite bin size used in the analysis biased the distribution. 



test verify that this did not introduce systematic errors. 

Figure [A.4| shows the fitted flux of a source as one moves it away from the viewing 
axis. Flux estimation is shown to be unaffected. 



A. 4 Source Significance and Spurious Detections 
A. 4.1 Significance 

As mentioned previously, the likelihood ratio (TS) is used to determine point source 
significance. In the context of EGRET data, the significance of a point source can be 
characterized in several ways, each applicable in different scenarios. 

The most direct measure of significance is the probability that a likelihood ratio 
test at a given position with no source present would result in a TS value greater 
than or equal to the measured value. Significance defined in this way is determined 
by the distribution of TS in the null hypothesis. Application of Wilks' theorem in 
this case predicts that TS will be distributed as Xi ( free parameters in the null 
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Figure A. 4: The distributions of best fit:true flux as one travels out toward the edge 
of the FOV. The software does not introduce a bias near the edge of the FOV. 

hypothesis=2 {gm,gb), free parameters in the alternate hypothesis = 3 {gm, gb,Ca))- 
For this particular distribution the significance is well approximated by = VTS 
where is the number of standard deviations. 

Figure |A.5| shows the result of several thousand independent likelihood ratio tests. 
The discrepancy between the simulated distribution and Xi is a direct result of the 
physically motivated constraint that point source flux should be positive. Statistical 
fluctuations are equally likely to result in point source like deficits as positive excesses 
and such negative fluctuations result in positive TS values if left unconstrained. These 
negative fluctuations do not correspond to physically meaningful point sources so they 
are assigned the value of zero flux and TS. The result is that at every point except 
TS=0 the resulting distribution is a factor of two lower than the Xi prediction and 
all the remaining probability density is taken up at the origin. The resulting single 
trial significance is correspondingly shifted to give slightly higher significance for a 
given value of TS. This effect is small. For example the TS value corresponding 
to a significance of 3cr is 9.0 whereas the prediction was 10.3. In general the simple 
expression = ^/TS remains a reasonable approximation. 
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Figure A. 5: The distribution of TS in the null hypothesis. No point sources were 
used in this simulation and the resulting distribution of TS is shown compared to xl- 
The discrepancy is exactly a factor 2 which results from the fact that negative fluxes 
are assigned TS=0. 

In searching for potential point sources, it is common to make many point source 
tests throughout the sky. While the significance of any one source is determined as 
above, the significance of this set of trials is reduced through multiple trials. Moreover, 
the extended PSF couples nearby likelihood ratio trials making successive trials not 
completely independent. Two trials separated by .5 degrees yielding elevated TS 
values is not strong evidence for two separate sources. What is done is to generate 
likelihood maps in which likelihood ratio trials were made at the center of every .5deg 
X .5deg pixel in an observation. An "excess" is defined to be a set of contiguous pixels 
having TS greater than or equal to some threshold. The significance of such an excess 
is best defined as the chance probability of such an excess occurring per unit solid 
angle. 

Figure [A.6| shows the distribution of such excesses in several all sky likelihood 



maps generated with no point sources added. The distribution of excesses represents 
the chance probabihty of finding a spurious source with specified TS or greater in the 
entire sky. The distribution is no longer well approximated by xl but is significantly 
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flatter. The number of low significance excesses is less than what would be predicted 
from the simplistic argument about the number of trials because the significant ex- 
cesses obscure less significant fluctuations. The number of sources with TS > 16 in 
the whole sky is ~ 2. The first EGRET source catalog (Fichtel et al. 1994) includes 
all excesses with TS > 16 and as a result should be relatively free of spurious sources 
whose nature is purely statistical. 
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Figure A. 6: The distributions of spurious sources. The dashed line shows the distri- 
bution of TS within a given VP while the dot-dashed line shows the distribution that 
results from the simultaneous fitting of all excesses with TS > 5. Also shown are the 
phase 1+2 source distribution showing clear evidence for non-statistical fluctuations. 

The independence of TS on factors other than the number of free parameters can 
also be verified in this way. The spatial distribution of excesses in the sky reveals 
a constant detection rate per solid angle across the sky with no correlation with 
exposure or background intensity (Figure |A.7| ). 



Comparing these simulated TS distributions with that obtained from actual flight 
data reveals the strong departure of the flight data from the null hypothesis. Fig- 
A.6| shows a comparison between the differential distributions of TS for simulated 



ure 



data with no point sources and for Phase I and II flight data. It demonstrates the 
strong evidence for non-statistical excesses. Some of these excesses are unambiguous 
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Figure A. 7: A distribution of TS in longitude for a simulation of the full sky with no 
input point sources. There is no evidence for any nonuniformity in the distribution 
of TS. 

such as the enormous significance of sources such as the Vela, Crab and Geminga, 
but at lower significance levels, the excesses present could be due to inaccuracies in 
the diffuse model as well as statistical fluctuations and weak point sources. 



A. 5 Computational Approximations: Ranai 

The evaluation and maximization of the likelihood function is computationally inten- 
sive. This computation time becomes prohibitive as one increases the number of pixels 
under analysis. As a result, an expedient approximation is customarily made: only 
pixels within an angular distance Ranai of a test point are considered when evaluating 
the likelihood function at that point. An additional advantage of this truncation is 
that it makes the point source analysis somewhat independent of large scale inaccu- 
racies in the background model. The optimal value of Ranai is one which minimizes 
computation time without introducing systematic analysis errors. 

A study of the effect of decreasing Ranai on the best fit model parameters is shown 
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in Table 2. There is little sensitivity to changes in Ranai down to values of ~ 8 degrees. 
Beyond this point point source significance is reduced while flux increases. This study 
provides justification for the conventional use of Ranai = 15 degrees. 



Table 1.1: Effect of Ranai on Parameter Estimation 



Ranai 


Input Source Counts 


Ca 


Vts 


9 m 


9b 


2 


521.4 


509.5 


3.8 


0.9 


6.9 


4 


521.4 


522.4 


6.1 


1.0 


1.5 


6 


521.4 


521.9 


7.2 


1.0 


1.1 


8 


521.4 


521.7 


7.8 


1.0 


1.0 


10 


521.4 


521.4 


8.2 


1.0 


1.0 


12 


521.4 


521.4 


8.4 


1.0 


1.0 


14 


521.4 


520.8 


8.6 


1.0 


1.0 


16 


521.4 


521.0 


8.8 


1.0 


1.0 


18 


521.4 


521.1 


8.9 


1.0 


1.0 


20 


521.4 


521.5 


9.0 


1.0 


1.0 



A. 6 Multiple Source Analysis 

The principal shortcoming of the hkelihood technique as it has been thus far described 
is that it is only capable of fitting one point source and diffuse emission. If point 
sources are sufficiently spaced ( > Ranai) this does not pose any problem but of 
course actual sources can be quite closely spaced. The solution to this problem has 
been to develop an iterative technique to effectively do multi-source fitting. The 
hkelihood model is then defined as, 

% = QmGij + QbEij + CaPSF{a, S)ij + CkPSF{ak, 5k)ij , (A.ll) 

k 

where Ck,CKk, Sk represent the best fit point source parameters of all sources other 
than the active source currently under optimization. This technique amounts to 
altering the diffuse model to include features that account for the other sources present 

in an observation. In practice what is done is to find the strongest point sources in 
an observation and add them to the diffuse model which is then used to search for 
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further sources. When these subsequent sources are found, the iterative technique 
is used to simultaneously fit all of them. This process is continued until no excesses 
with TS > some threshold. 

This technique was tested to see if it accurately reproduced a set of sources inserted 
into the model. Tables 3,4 show a comparison between the source parameters input 
into a simulated observation of the Virgo region and what was found upon analysis 
of this data. It would have been preferable to do an analysis of a statistical ensemble 
of such observations but the amount of computing time required made this approach 
impossible. The analysis was carried down to a threshold of TS > 12. No spurious 
sources appear in this observation . The average number of spurious sources with 
TS > 12 in this amount of solid angle is ~ 2 so this is not unhkely. The fluxes and 
positions of the true sources are all within the respective error boxes of the measured 
data. The analysis demonstrates the power of the likelihood technique even when the 
sources must be dealt with simultaneously. 



Table 1.2: Sources Input into Virgo Region (from Phase 1 Analysis) 



EGRET Name 


Identification 


Position 


(J2000) 


Flux 






ra 


dec 


{IQ-^cts cm-^s-i) 


GRO J1256-05 


3C 279 


194.06 


-5.79 


170 


GRO J1232+02 


3C 273 


188.25 


2.82 


20 


GRO J1224+22 




186.08 


22.75 


13 


GRO J1230-02 




187.69 


-2.96 


10 



Table 1.3: Results of Multisource Fit for Sources with TS > 12 



Name 


Position 


Err. 


arcmin) 


Flux 


TS 


Ca 


9m 


9b 




ra dec 


68% 


95% 


val. err. 










GRO J1256-05 


194.12 -5.76 


5. 


7. 


167.6 5. 


2432.8 


1301.3 


1.37 


0.64 


GRO J1231+03 


187.94 3.12 


20. 


33. 


21.3 2.7 


102.7 


185.9 


1.01 


0.95 


GRO J1224+22 


186.08 22.80 


21. 


48. 


11.1 2.8 


24.2 


66.1 


0.42 


1.46 


GRO J1227-02 


186.97 -2.60 


35. 


53. 


9.4 2.3 


23.2 


75.7 


1.45 


0.64 



The distribution of spurious soiuxcs needs to be evaluated for this technique. It 
is easy to see that the distribution of excesses will differ from the simple case of 
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one source fitting. Whereas the distribution of excesses for the first hkelihood map 
is known in the null hypothesis, after altering the diffuse model to match this data 
and reevaluating the likelihood map, the new distribution will differ. Monte Carlo 
techniques have allowed the direct measurement of the probability per unit solid angle 
of finding an excess with TS greater than some threshold. 

What was done was to generate simulated observations with no point sources 
present and to analyze these data in exactly the same manner as flight data. Like- 
lihood maps are generated and strong sources identified and fit. These sources are 
then added to the diffuse map and new likelihood maps are generated. This process 
is repeated until no excesses above some threshold are present. The results are shown 
in Figure |A.(j| . At higher values of TS, the results are consistent with what was found 
for single source analysis. As one gets to lower values of TS the distribution rises 
more steeply as one would expect given the fact that significant fluctuations are no 
longer allowed to obscure smaller ones. 
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